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FOREWORD 


This  report  was  prepared  by  the  Systems  Engineering  Division  of  Daniel, 
Mann,  Johnson,  &  Mendenhall  under  Contract  AF  04(61 1)- 7540.  This 
contract  authorizes  studies  to  be  made  of  test  systems  for  toxic  propellant 
rockets.  This  volume  presents  the  meteorological  aspects  of  the  studies 
as  found  during  the  contractual  period  covering  1  August  1961  through 
31  October  1962,  and  is  part  of  the  final  report  of  the  contract. 

Conclusions  as  to  the  specific  use  of  the  results  of  the  "Air  Dispersion  of 
Rocket  Exhaust"  are  presented  in  Report  SSD- TDR- 62-  1 37,  entitled  "Design 
of  a  Toxic  Rocket  Test  System". 
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ABSTRACT 


An  analysis  of  the  ability  of  the  atmosphere  to  dissipate 
toxic  exhaust  products  from  a  rocket  lest  system  was  per¬ 
formed,  The  analysis  covered  dispersion  stack  height, 
total  integrated  dosage,  inversion  layer  penetration,  and 
sampling  and  weather  data  systems.  Methods  for  determin¬ 
ing  toxicity  levels  downwind  of  a  toxic  emission  were 
developed.  The  conditions  for  safely  releasing  various 
quantities  of  toxic  materials  were  defined  for  various  types 
of  rocket  firings. 
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INTRODUCTION 


Large  increases  in  payload  or  range  are  available  with  missiles  and  space¬ 
craft  from  relatively  small  increases  in  the  specific  impulse  of  the 
propellants.  Thus,  there  is  a  continuous  search  for  propellants  and  methods 
to  increase  specific  impulse  levels.  This  search  for  rocket  fuels  which 
contain  higher  energy  has  led  researchers  to  many  propellant  combinations 
which  are  extremely  toxic.  Testing  these  highly  toxic  rockets  at  the  Edwards 
Rocket  Site  presents  a  problem  in  toxicity  which  may  involve  the  health 
and  lives  of  people  living  in  surrounding  communities  as  well  as  personnel 
at  the  base. 

Since  it  is  imperative  to  our  national  interests  to  perform  further  propellant 
research,  it  is  of  the  highest  importance  that  toxic  rockets  are  tested  and 
developed.  It  follows  that  a  means  of  safely  testing  such  rockets  must 
be  evolved.  The  6593rd  Test  Group  at  the  Edwards  AFB  Rocket  Site  has 
been  assigned  Air  Force  responsibility  for  this  mission,  and  the  purpose 
of  this  report  is  to  present  the  theory  and  application  of  meteorology  to 
assist  in  the  accomplishment  of  this  goal. 

The  study  covered  in  this  report  presents  a  means  of  determining  the  extent 
of  ground  level  contamination  caused  by  the  rocket  effluent  discharged  into 
the  atmosphere.  The  primary  purpose  for  this  meteorological  study  is  to 
provide  a  basis  for  facility  design  for  testing  toxic  rocket  propellants. 

Trade  offs  in  cost  and  time  may  be  accomplished  between  the  two  extremes 
of  complete  filtration  or  decontamination  of  the  rocket  effluent  and  the 
direct  discharge  of  effluent  into  the  atmosphere.  An  elaborate  and 
extensive  meteorological  network  and  accurate  weather  prediction  may  be 
required  to  permit  safe  open-air  firing.  The  trade  -offs  in  design  must 
at  all  times  take  into  consideration  the  maintenance  of  a  high  degree  of 
safety. 

Section  I  of  this  report  covers  the  theoretical  and  em.pirical  derivations  of 
the  techniques  for  analysis  of  movement  of  clouds  of  toxic  products  in 
the  atmosphere.  The  study  includes  an  analysis  of  both  the  equations  and 
the  coefficients  according  to  present  state-of-the-art.  Atmospheric 
dispersion  phenomena,  effective  stack  height,  inversion  penetration,  and 
the  total  integrated  dosage  theories  were  evaluated  and  are  presented. 

Section  II  of  this  report  covers  the  application  of  the  theories  to  hardware 
and  support  system  designs  and  to  the  parameters  which  must  be  known 
for  these  designs. 
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SUMMARY 


The  theories  and  empirical  findings  concerning  the  capacity  of  the 
atmosphere  to  dilute  man-made  air  pollutants  were  evaluated  from  the 
standpoint  of  dispersing  rocket  exhaust.  The  primary  problem  found  in 
this  evaluation  centered  around  the  physical  size  of  the  rocket  engines 
under  consideration.  The  larger  size  (5,  000  to  100,  000  pounds  thrust) 
rockets  require  meteorological  data  of  a  magnitude  higher  than  micro¬ 
meteorology  but  a  magnitude  lower  than  found  in  atomic  explosions. 

The  general  considerations  covered  in  this  study  were  the  various  dispersion 
processes  encountered  in  diffusion,  effectii’e  stack  height,  inversion 
penetration  and  in  total  integrated  dosages  of  toxic  materials.  The  rate 
and  type  of  diffusion  are  influenced  by  wind  speed,  time  variation  in  wind 
direction,  gustiness,  vertical  stability  or  lapse  rate,  and  ground  surface 
influences.  Dispersion  theories  and  formulae  investigated  were  Sutton's, 
Bosanquet  and  Pearson's,  Holland's  Gaussian  form,  and  the  Ocean  Breeze 
derivation.  Continuous  point  source  dispersion  formulae  of  all  above 
types  are  essentially  similar  and  can  be  equated  to  each  other.  The 
reliability  and  the  procedure  for  application  of  these  formulae  are 
presented.  Any  of  the  formulae  may  be  used  with  equal  reliability, 
depending  on  the  particular  type  analysis  being  conducted. 

In  addition  to  the  dispersion  process,  the  effective  height  to  which  a 
rocket  exhaust  plume  will  rise  must  be  known  for  the  determination  of 
downwind  contamination.  The  factors  influencing  this  plume  rise  are: 
rocket  test  size,  propellant  characteristics,  rocket  firing  orientation, 
and  atmospheric  conditions.  A  number  of  equations  are  available  for  the 
determination  of  plume  rise.  Each  formula  has  been  discussed,  and 
optimum  relations  have  been  chosen. 

The  equation  which  will  be  used  for  stack  height  determination  depends 
on  the  type  of  emission.  A  hot  rocket  exhaust  pointed  upward  will 
require  the  use  ol  an  equation  for  momentum  rise  and  an  equation  for 
buoyancy  rise.  The  directing  of  the  exhaust  horizontally  will  omit  the 
need  of  the  momentum  rise  equation  and,  therefore,  reduce  the  overall 
effective  stack  height. 

The  above  relations  evaluated  for  stack  height  determination  are  for  the 
continuous  emission  of  an  effluent.  The  possibility  of  a  detonation 
occurring  during  a  rocket  test  firing  is  always  present  and  must  be 
provided  for.  Therefore,  the  buoyant  rise  of  explosion  clouds  was 
considered.  The  U.  S.  Atomic  Energy  Commission,  in  extensive  weapons 
tests,  has  accumulated  a  large  quantity  of  data  on  the  rise  of  nuclear 
clouds  as  a  function  of  weapon  energy  yield.  This  information  has  been 
compiled  into  graphical  form  for  the  rapid  estimation  of  cloud  rise  from 
explosion  sources. 
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As  a  result  of  the  effective  stack  height  studies,  an  empirical  expression 
was  formulated  to  cover  the  range  of  conditions  expected  in  rocket  test 
work.  In  addition,  an  equation  was  derived  by  dimensional  analysis 
which  includes  all  factors  affecting  cloud  rise.  The  empirical  expression 
may  be  utilized  (with  appropriate  safety  factors)  for  immediate  applications, 
and  the  broader  equation  derived  by  dimensional  analysis  will  provide  a 
higher  degree  of  reliability  when  directly  applicable  data  become  available. 

Inversions,  conditions  of  negative  lapse  rate  in  which  air  tends  to  mix  at 
much  lower  rates  than  normal,  have  been  studied.  The  indications  are 
that  the  energy  available  from  rocket  tests  is  sufficient  to  penetrate 
inversions  when  the  thrust  level  is  above  a  minimum  value. 

Applications  to  the  testing  of  toxic  rocket  engines  have  been  examined.  The 
thrust  levels  coi'ered  by  the  equations  presented  are  within  the  range  of 
6,  000  to  100,  000  pounds  thrust  for  solid  propellant  grains  of  100  to 
1  1,  000  pounds  weight.  The  formulae  presented  for  diffusion  are  applicable 
to  any  combination  of  pollutants,  whether  they  be  gaseous,  liquid  vapor, 
or  solid  phases.  This  means  that  the  presentation  covers  the  normal 
testing  of  liquid  or  solid  propellants  with  solid  additives,  test  stand 
explosions,  and  spillage  of  liquids.  In  addition,  the  diffusion  of  effluents 
fromi  nuclear  explosions  or  reactors  can  be  calculated  by  using  the 
equations  presented  in  this  report. 

In  many  cases,  the  cost  of  toxic  rocket  test  systems  may  be  appreciably 
reduced  by  taking  advantage  of  meteorological  phenomena.  Due  to  the 
absence  of  directly  available  experimental  data,  methods  of  improving  the 
dispersion  and  plume  rise  formulae  are  necessary  to  enable  the  maximum 
utilization  of  meteorology  in  test  systems  for  a  toxic,  but  otherwise 
desirable,  propellant.  This  improvement  can  be  accomplished  with  a 
suitable  micrometeorological  network  including  a  sampling  system  and 
data  reduction  equipment. 
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SECTION  I 


ANALYTICAL  CONSIDERATIONS 


This  section  of  the  report  covers  the  analytical  and  theoretical  treatments  of 
the  action  of  a  cloud  of  pollutants  in  the  atmosphere.  Extensive  literature 
research  was  conducted  on  the  metcorulogicai  aspect  of  toxic  rocket  testing. 
Properties  of  meteorology  which  have  been  coi’ered  include  the  general  con¬ 
siderations  affecting  atmospheric  dispersion.  To  this  end,  a  study  of  the 
reflection  of  a  cloud  as  it  reaches  ground  le\el,  the  effects  of  stack  heights, 
and  the  penetration  of  inversions  was  performed.  Dispersion  is  governed 
by  properties  which  include: 

A.  Wind 

B.  Daily  and  Periodic  Weather  Changes 

C.  Gaseous  Diffusion 

The  concentration  of  a  pollutant  downwind  of  its  source  is  dependent  upon 
these  influences,  as  well  as  the  interactions  between  them. 

Many  different  formulae  for  the  calculation  of  effecti\-e  stack  height  have  been 
covered,  and  definite  choices  ha\  e  been  made  as  to  the  optimum  relations. 
Inversion  penetration  is  basically  a  part  of  the  stack  height  analysis,  but  the 
atmospheric  conditions  are  sufficiently  different  to  cause  it  to  be  a  distinct 
problem  in  itself. 

The  variability  of  parameters  used  in  the  diffusion  equations  is  presented  in 
addition  to  the  above  analyses.  The  direction  in  which  the  toxic  concentration 
varies  is  a  complex  function  of  many  variables;  ail  of  which  should  be  known 
to  a  high  degree  of  accuracy  to  attain  the  optimum  design  of  an  actual  facility 
as  formulated  upon  meteorological  consideraLons. 


GENERAL  METEOROLOGICAL  CONSIDERATIONS 


The  concentration  of  a  pollutant  downwind  of  its  release  into  the 
atmosphere  is  dependent  upon  many  variables.  The  pollutant  is  trans¬ 
ported  with  the  air  mass  and  is  acted  upon  by  the  dynamic  forces  of  the 
air.  Instability  of  large  air  masses  is  caused  by  differential  surface 
heatingj  rotation  of  the  earth,  surface  friction,  mechanical  mixing  of 
moving  air  masses  by  surface  irregularity,  temperature  variation  with 
altitude,  and  a  host  of  other  influences  as  well  as  interactions  between 
those  various  influences. 

The  surface  air  behavior  can  only  be  determined  by  measurement  of  its 
performance  over  a  period  of  time.  By  application  of  the  data,  together 
with  an  under  standing  of  fluid  dynamics  as  it  is  applied  to  meteorology 
and  micrometeorological  phenomena,  reasonable  predictions  of  atmos¬ 
pheric  behavior  can  be  made.  With  this  background  of  accumulated  data, 
weather  knowledge,  and  current  weather  information,  much  can  be 
predicted  about  the  behavior  of  pollutants  which  might  be  released  into 
the  atmosphere 

Micrometeorology  is  normally  concerned  with  the  layer  of  air  within 
the  first  hundred  meters  above  the  earth's  surface  and  for  a  horizontal 
distance  of  a  thousand  or  so  meters  When  a  pollutant  is  very  toxic,  the 
areas  of  interest  become  greater  and  may  extend  to  the  stratosphere  and 
cover  the  hemisphere  as  is  the  case  with  atomic  explosions.  When  applied 
to  rocket  engine  test  and  missile  launch,  micrometeorology  may  involve 
large  areas  or  may  be  very  limited  depending  upon  the  toxicity  of  pro¬ 
pellants  or  exhaust  products  and  the  quantities  concerned.  As  rocket 
technology  is  advanced,  the  use  of  toxic  propellants  is  increased  since 
higher  reactive  fuels  and  oxidizers  tend  to  be  toxic.  Even  so-called 
nontoxic  fluids  handled  in  large  quantities  can  smother  human  life  by 
displacement  of  air  and  the  oxygen  needed  for  support  of  life.  Other 
propellants  can  be  lethal  at  levels  as  low  as  one  part  per  million  parts 
of  air  or  evep  one  part  per  billion  may  be  in  excess  of  tolerable  levels 
in  inhabited  areas.  Meteorological  information  for  extended  areas  must 
be  known  where  very  toxic  substances  are  released  into  the  atmospherCr 
whether  by  accident  or  of  necessity. 

Information  about  the  potential  source  of  atmospheric  contamination 
includes  toxicity  levels,  nature  of  the  toxic  reaction,  quantities  involved, 
relative  density  between  pollutant  and  air,  and  particle  size  of  pollutants 
where  solids  are  involved.  The  toxicity  level  of  the  pollutant  determines 
the  area  of  involvement  and  extent  of  meteorological  details  required  to 
predict  safe  atmospheric  dispersion. 
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Atmospheric  dispersion  processes  include  differential  movement  between 
parcels  of  air  which  results  in  mixing  of  adjoining  parcels.  This  causes 
dispersion  and  dilution  of  suspended  vapors  and  solids  contained  in  the 
air  parcels.  The  mixing  is  primarily  microscopic  in  that  intermingling 
of  air  masses  results  in  mixing  of  a  portion  of  one  air  mass  with  another 
where  individual  portions  slowly  lose  their  identity.  Air  masses  of 
different  density  retain  their  identity  with  the  more  dense  masses  re¬ 
maining  in  contact  with  the  surface  of  the  earth.  Resistance  to  inter¬ 
mingling  expresses  itself  in  the  formation  of  weather  fronts.  A  cold- 
front  is  formed  when  a  cold  air  mass  "under  runs"  a  warm  air  mass 
forcing  the  v.arm  air  upward  and  creating  very  turbulent  conditions. 

Rainfall  formed  from,  such  a  front  is  of  an  intermittent  "downpour"  type 
because  of  the  rapid  rise  of  warm  air  by  the  relatively  steep  cold  air 
wedge.  Warm-fronts  arc  produced  when  warm  air  blows  over  the  top 
of  the  trailing  edge  of  a  cold  air  mass.  The  trailing  edge  is  of  small 
slope  being  stretched  out  by  friction  between  the  cold  air  mass  and  the 
ground  surface.  Rainfall  produced  along  a  warm-front  is  continuous  but 
normally  light.  Smaller  parcels  of  air  mixed  together  tend  to  retain 
their  identity  for  a  period  of  time  but  always  eventually  lose  their  identity. 

Diffusion  processes  are  influenced  by  wind  speed,  variation  in  wind 
direction  and  velocity  igustiness’.,  and  vertical  stability  as  measured  by  the 
vertical  temperature  profile  Itapse  rate',  and  ground  surface  influences. 
Measurable  w-.nd  speeds  normally  result  in  turbulent  flow  where  there  is 
a  tumbling  action  rather  than  streamlined  flow.  Gustiness  results  in 
mixing  of  air  at  a  rate  significantly  greater  than  that  produced  by  molecular 
motion  and  resulting  Brownian  movement.  Because  of  Brownian  motion, 
diffusion  of  gases,  •■apors,  and  suspended  solids  will  always  be  greater 
than  zero.  Diffusion  by  Brownian  motion  result's  from  random  motion  of 
the  molecules.  The  net  transport  of  a  property  in  a  given  direction  is 
proportional  to  the  gradient  in  that  direction.  Eddy  currents  in  the  air 
result  in  an  apparent  random  motion  which  tends  to  distribute  a  property 
uniformly  by  reducing  its  gradient.  The  molecular  diffusion  equation  is 
extended  by  analogy  to  turbulent  diffusion,  but  the  rate,  of  diffusion  must 
be  determined  from  field  observations.  Exponential  or  Gaussian  dis¬ 
tribution  of  a  property  is  assumed  for  turbulent  diffusion  and  is  fairly 
well  confirmed  by  field  tests  where  known  quantities  of  tracer  substances 
are  released  in  the  atmosphere  and  samples  collected  using  relatively 
long  sample  periods. 

The  temperature  lapse  rate  is  defined  as  the  decrease  of  temperature 
with  increase  in  height.  When  the  lapse  rate  is  a  negative  value,  where 
temperature  increases  with  height,  the  condi':ion  is  called  an  inversion. 

The  atmospheric  lapse  rate  is  an  indication  of  atmospheric  stability  which 
affects  vertical  diffusion.  The  higher  the  lapse  rate,  the  greater  the 
vertical  diffusion.  A  negative  lapse  rate  'inversion  conditions)  results  in 
very  stable  air  layers  and  reduced  vertical  dispersion.  The  precise 
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relationship  between  dispersion  rates  and  atmospheric  lapse  rate  is  not 
established,  but  an  approximate  relationship  is  incorporated  in  dispersion 
equations,  One  difficulty  is  the  change  of  lapse  rate  with  altitude  as  well 
as  with  time,  humidity  conditions,  surface  heating,  wind  speeds,  and 
mechanical  mixing  by  surface  features.  Inversion  layers  are  often  formed 
near  the  earth's  surface  at  night  when  loss  of  heat  by  radiation  to  the  ground 
results  in  a  transfer  of  heat  from  the  lower  air  layer  to  the  ground  and  the 
temperature  of  the  lower  air  layer  drops  below  that  of  the  air  above,  A 
cloud  cover  reradiates  heat  back  to  the  ground  surface,  reduces  the  loss 
of  heat  by  the  surface,  and  minimizes  the  formation  of  an  inversion  at  the 
surface.  However,  the  cloud  cover  may  be  caused  by  a  weather  front 
producing  a  high  level  inversion. 

A  parcel  of  air  displaced  upward  will  expand  and  cool  adiabatically  as¬ 
suming  no  transfer  of  heat  between  the  parcel  and  the  surrounding  air. 

The  dry  adiabatic  lapse  rate  is  the  drop  in  temperature  with  change  in 
height  of  dry  air  when  displaced  upward.  The  dry  adiabatic  lapse  rate 
is  5.  5  F  per  1000  feet.  A  dry  parcel  of  air  displaced  vertically  in  at¬ 
mospheric  air  with  an  adiabatic  lapse  rate  will  have  neutral  buoyancy 
since  its  temperature,  pressure,  and,  therefore,  density  will  change  at 
the  same  rate  as  the  surrounding  air. 

Stability  of  the  atmosphere  affects  dispersion  of  a  pollutant  not  only  by 
reducing  vertical  dispersion,  but  during  stable  conditions,  diffusion 
transverse  to  the  plume  axis  is  usually  reduced.  Where  the  pollutant  is 
buoyant  due  to  lov  molecular  weight  or  reduced  density  because  of  high 
temperature ,  the  effluent  will  rise  above  the  actual  stack  and  produce  an 
effective  stack  height  greater  than  the  physical  stack.  When  the  atmos¬ 
phere  is  very  stable  due  to  isorhermal  conditions  or  a  negative  lapse,  the 
plume  buoyanc'  rapidly  decreases  wr.h  increased  heigh'  since  the  at¬ 
mospheric  density  decreases  fastei  than  the  pl...me  wiiic'-  cools  adiabatically. 

For  treatment  of  dispersion  problems,  atmospheric  stability  is  divided 
into  several  types  from  very  unstable  {superadiabatic  lapse  rate)  to  very 
stable  (strong  inversion),  Unstable  atmospheric  conditions  can  cause 
looping  of  an  exhaust  plume.  Looping  occurs  with  a  high  degree  of  turbu¬ 
lence,  especially  with  convective  turbulence,  and  is  typical  of  a  daytime 
condition  with  intense  solar  heating  of  the  earth's  surface.  Coning  plumes 
occur  under  more  nearly  neutral  thermal  conditions  when  mechanical 
mixing  of  small  scale  predominates,  Coning  is  more  likely  to  occur  when 
a  cloud  cover  reduces  thermal  effects  either  by  reducing  incoming  solar 
radiation  at  daytime  or  outgoing  terrestrial  radiation  at  night.  The  half¬ 
angle  of  the  coning  plume  is  on  the  order  of  10°,  from  canter  line  to  edge. 
Fanning  plumes  occur  under  stable  conditions  when  mechanical  turbulence 
is  suppressed.  The  vertical  component  is  suppressed  more  than  the 
horizontal  and  the  plume  is  wider  horizontally  than  vertically.  Fanning  is 
most  likely  to  occur  at  night  when  the  earth's  surface  is  cooled  by  outgoing 
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Fig 


ouliav-c 

Looping  -  Large  Lapse  With 
Stack-Low  Heat 


Coning  -  Smal]  Lapse  With 
S''2ck-Low  Heat 


2.  Various  Plume  Types 


radiation.  Fumigation  is  caused  by  unstable  air  at  and  below  the  plume 
with  stable  air  above.  Dispersion  of  the  plume  downwind  brings  the 
exhaust  to  ground  level  in  a  relatively  short  distance  before  dispersion 
can  become  effective  Lofting  is  caused  by  stable  air  below  the  plume 
and  turbulence  above.  The  turbulence  above  results  in  effective  dis¬ 
persion  yielding  very  low  ground  surface  concentrations.  Lofting  will 
occur  when  the  plume  is  released  above  or  penetrates  an  inversion, 

On  the  other  hand,  fumigation  will  result  if  the  inversion  is  not  penetrated 
and  ground  turbulence  is  present  due  to  thermal  effects  from  ground 
surface  warm-up  and  inversion  break-up.  The  described  plumes  are 
samples  of  conditions  which  may  be  found  during  actual  operation  of  a 
facility,  but  the  capacity  to  test  must  be  determined  at  the  given  test 
site  and  at  the  time  of  test. 

Stable  conditions  will  keep  a  plume  compact  and  prevent  contact  with 
the  ground  for  a  considerable  downwind  distance  if  its  density  if  not 
greater  than  that  of  the  surrounding  atmosphere.  Contact  of  the  plume 
with  the  ground  and  the  point  of  maximum  concentration  occur  closer 
to  the  point  of  release  with  increasing  atmospheric  instability.  For 
extreme  instability  at  low  wind  speeds,  portions  of  a  looping  plume  may 
reach  the  ground  within  a  distance  of  about  one  stack  height.  High  wind 
speeds  increase  mechanical  mixing  caused-by  ground  features  and  reduce 
thermal  effects.  Plume  features  which  arise  from  thermal  effects  are, 
consequently,  more  likely  to  occur  at  low  wind  speeds. 

The  following  discussion  on  diffusion,  dispersion,  effective  stack  heights, 
and  inversion  layer  penetration  includes  the  mathematical  techniques 
utilized  in  explaining  these  phenomena. 
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ATMOSPHERIC  DIFFUSION  FORMULAE 


Diffusion  is  the  motion  of  fumes  or  dust  particles  in  the  atmosphere.  The 
principles  of  diffusion  or  dispersion  are  well  established  Many  theoretical 
analyses  on  this  phenomena  have  been  prepared.  In  addition,  a  vast  amount 
of  experimental  data  has  been  accumulated. 

Diffusion  is  sometimes  limited  to  fluids  in  suspension  and  dispersion  to 
the  activity  of  particulate  matter  in  suspension.  Diffusion  due  to  molecular 
activity  is  small  compared  with  mixing  and  other  effects  due  to  atmospheric 
turbulence.  A  smoke  plume  released  in  the  atmosphere  will  spread  outward 
by  the  energy  within  the  plume  and  the  activity  of  the  atmosphere.  The  con¬ 
centration  of  smoke  will  be  greatest  near  the  center,  decrease  toward  the 
edges,  and  diminish  along  the  axis  of  the  plume  with  time  or  with  distance 
from  the  source.  The  degree  of  stability md  velocity,  and  gustiness 
affect  the  rate  of  dispersion  of  the  plume  in  the  air.  The  buoyancy  of  the 
plume  also  affects  the  rate  of  dispersion,  but  the  most  important  effect  of 
buoyancy  is  that  of  raising  the  plum.e  above  the  lower  air  'ayers.  By  the 
time  the  edge  of  the  plume  touches  the  ground,  diffusion  will  have  been 
effective  in  considerably  reducing  the  concentration  of  the  cloud, 

This  discussion  is  intended  to  present  the  formulae  required  for  use  with 
toxic  rocket  testing.  Emphasis  will  be  placed  on  problems  involved  with 
rockets  using  toxic  propellants  and  those  yielding  toxic  exhaust  products 
in  order  to  provide  maximum  usefulness  for  rocket  test  system  design 
and  operation.  Spills  tests  have  been  made  on  liquid  propellants  in  which 
toxic  liquids  have  been  allowed  to  evaporate  with  and  without  combustion 
either  by  air  or  in  conjunction  with  other  oxidizers  for  example, U.DMH+NTO), 
For  liquid  propellant  spills  accompanied  by  detonation  the  source  strength 
(Q)  will  be  instantaneous  and,  therefore,  the  dow'nwind  toxic  concentration 
(X)  can  be  easily  found.  This  also  holds  true  of  solid  pi  opellant  detonations. 
When  a  liquid  spill  simply  evaporates  or  burns  for  a  period  of  time,  variables 
of  spill  volume,  surface  area  exposed  to  the  atmosphere,  liquid  temperature 
and  vapor  pressure  or  flame  temperature  and  heat  of  combustion,  and  wind 
velocity  must  be  known  to  determine  the  continuous  source  strength. 

Atmospheric  diffusion  has  been  expressed  as  a  differential  equation  which 
takes  into  account  the  rate  of  diffusion  in  the  three  dimensions  of  space. 

The  major  problem  is  the  determination  of  dispersion  rates  in  the  air. 

Fick's  law  of  molecular  diffusion  (Reference  1)  may  be  stated  as  follows; 
"Diffusion  of  material  is  in  the  direction  of  decreasing  concentration  and 
is  proportional  to  the  concentration  gradient.  "  Atmospheric  diffusion  is 
primarily  caused  by  turbulence  or  gustiness  of  the  atmosphere  and  is 
influenced  by  vertical  stability  as  indicated  by  the  lapse  rate.  The  influence 
of  atmospheric  turbulence  in  changing  the  properties  of  a  smoke  cloud  or 
plume  can  be  compared  to  that  of  molecular  activity  where  the  time  interval 
is  sufficient  to  average  out  random  variations.  By  fitting  atmospheric 
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observ^it  ions  to  busic  differential  equations  the  resulting  diffusion 
v-ilues  will  vary  from  about  0„  2  cm^  /sec.  for  molecular  diffustion  to  10^° 
cm^  /sec  for  diffusion  due  to  large-sc-ilc  cyclonic  storms  in  the 
atmosphere  The  pronouiiced  influence  of  conditions  of  the  atmosphere 
diffusion  has  resulted  in  the  development  by  Sutton  and  Frenkiel 
iK.d'ercnc  c  2)  of  tim  statis'ic  il  theories  of  turbulent  diffusion,  following 
Tfivlor  s  st.itist  ici  1  tlieory  of  turbulence  'Reference  R.  Sutton’s 
statistic-!  theory  seems  to  apply  to  intermediate  range  atmospheric 
dispersion  problems  wliere  hundreds  of  meters  to  kilometers  are 
inv  ul ved 

A  pronounced  sliear  effect  on  the  mem  wind  is  exhibited  for  layers  of 
iir  near  tht-  ground  Tins  is  due  to  tlie  wind  s  increasing  from,  zero  at 
ground  level  to  relatively  high  velocities  .it  several  thousand  feet  in 
elev-^tion.  Tlie  shearing  -jction  affects  dispersion  by  elongating  the 
upper  portion  ol  a  smoke  cloud  or  plume  The  ground  and  structural 
fe  Oures  of  tlie  ground  influence  dispersion  by  mechanical  mixing  of 
‘f.e  air  Diffusicjii  ('quations  iiave  bi'cn  cb'rived  based  on  this  shearing 
action  bijf  are  limi'ed  in  -qiplic.-t  ion  to  the  lower  air  levels.  Statistical 
dot  ennin  I'lon  of  disjiersion  pa  r-jmi‘tcrs  appears  to  be  more  reliable  than 
the  theoretic.:!  i‘xp>ressions  -md  is  easily  determined  at  or  near  the 
earth  s  surface  There  is  cun fi rm.a t ion  that  concentrations  in  a  diffusing 
cloud  a  re  distributed  according  to  a  three-dimensional  Gaussian  law 
tReference  1.  Frenkiel  -idvanccd  the  fol'ouing  equation  for  the  concen¬ 
tration  distribution  from  an  instan'aneoiis  point  source  in  an  infinite 
region 


y ' 


X ,  N  v  z  t 2  tt  y  ) 


3/2 


The  coordin-.te  axes  arc  imagined  to  be  affixed  to  the  center  of  the  cloud 
and  nre  translated  wi'h  the  mean  wind  value  For  nonisotropic  diffusion 
the  expression  may  be  written  as  follows 
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4.  Downwind  Concent  rati 


Sutton's  instantaneous  point  source  equation  for  the  nonisotropic 
case  is: 


(^.y,  C  C  (ut)^(^"^)/^ 

X  y  z 


-(ut) 


Sutton  considered  th.il  the  ground  could  act  as  a  perfect  reflector  of  diffusing 
particles  and  accounted  for  it  by  using  the  "method  of  images.  "  By 
doubling  the  right  hand  side  of  the  expressions,  this  reflection  influence 
can  be  accounted  for  when  determining  ground  level  concentrations.  This 
reflection  effect  can  be  expected  for  gases  that  do  not  react  chemically 
with  the  ground  or  vegetation  and  small  paniculate  matter.  The  effect 
of  thermal  stability  on  atmospheric  dispersion  in  Sutton'  s  equation  is 
represented  by  a  stability  factor  "n"  which  ranges  from  almost  zero  for 
extremely  turbulent  conditions  to  a  value  of  1  as  a  limit  for  extremely 
stable  conditions.  Sutton's  formula  has  been  successfully  utilized  to 
predict  diffusion  over  distances  up  to  several  kilometers  under  certain 
meteorological  conditions,  and  this  formula  and  modifications  thereof 
have  been  fairly  widely  accepted. 


Bosanquet  and  Pearson  (Reference  3)  derived  the  following  formula  for 
ground  level  concentration  distribution  from  a  continuous  elevated  point 
source  near  the  ground: 


L— 

p  X 


Where  p  and  q  are  vertical  and  lateral  diffusion  coefficients.  The 
equation  is  very  similar  to  Sutton's  continuous  point  source  formula 
where  extremely  turbulent  conditions  exist  (  n  =  0  ). 

Various  point  source  formulae  may  be  integrated  with  respect  to  time 
to  give  equations  for  concentration  distribution  downwind  from  a  con¬ 
tinuous  point  source.  Sutton's  equation  for  continuous  elevated  point 
source  yields  the  following  expression  when  integrated; 


X 
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Continuous  line  source  formulae  can  be  derived  from  point  source 
formulae.  Sutton's  continuous  infinite  elevated  crosswind  line  source 
for  nonisotropic  diffusion  is  as  follows: 

/  h'*  \ 


This  formula  may  be  applied  to  a  row  of  stacks  or  facility  roof-level 
emission.  For  rocket  testing,  where  intermittent  firings  and  a  limited 
number  of  test  cells  are  involved,  the  crosswind  line  source  would 
seldom  apply. 

All  of  the  above  equations  have  been  developed  based  on  gaseous 
diffusion.  If  a  smoke  plume  has  any  solid  or  liquid  particles  within 
it  the  values  calculated  for  elevated  sources  will  be  in  error  by  a 
function  of  the  settling  rate  of  the  particles.  The  settling  rate  of  particles 
varies  with  the  particle  diameter  and  density  and  may  be  determined 
by  use  of  Stoke's  law,  by  other  suitable  formulae,  or  by  observation. 

In  general,  the  simplest  solution  is  to  define  an  average  particle  size 
and  determine  the  settling  velocity  by  Stoke's  law.  This  particle 
settling  will  effectively  tilt  the  axis  of  the  plume  (from  a  continuous 
point  source)  downward  so  that  the  ground  will  be  intersected  at  an 
angle  (6).  An  expression  for  this  angle  can  be  written  as  a  function  of 
settling  rate  (v)  and  wind  velocity  (u) . 

tan  0  =  -^ 

u 


This  angle  will  generally  be  rather  small.  The  elevation  of  the  plume 
(z')  may  then  be  approximated  at  any  distance  (x)  downwind  from  the 
maximum  effective  stack  height. 

z'  =  h-4^ 
u 
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The  continuous  elevated  point  source  formula,  by  the  substitution  of 
z'  for  h  in  the  power  of  e,  becomes:  2 
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The  factor  of  2  has  been  dropped  since  complete  deposition  is  now 
assumed,  rather  than  partial  reflection. 

Holland  derived  the  following  formula  for  an  instantaneous  volume 
source  using  a  Gaussian  concentration  distribution,  in  terms  of  Sutton's 
theory, 


_ 2_Q _ 

tt"/"  C^(xo  +u 


The  term  is  the  distance  upwind  from  the  stack  to  the  virtual  source. 
For  distances  of  many  miles  it  may  be  neglected  in  computing  X  . 


The  Holland  derivation  for  an  instantaneous  volume  source  can  also  be 
modified  for  particle  diffusion  from  an  elevated  source.  By  analogy  with 
the  continuous  elevated  point  source  formula  it  will  become: 

I—  /  V  X'2  -I 
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The  factor  of  2  has  been  dropped  since  complete  deposition  is  now 
assumed,  rather  than  any  reflection. 

In  the  present  use  of  these  equations  for  predicting  the  fallout  from 
rocket  testing,  the  particle  size  is  in  the  range  of  less  than  one  micron 
up  to  about  20  microns  in  diameter.  Generally  the  particles  will  be 
about  1  micron  for  normal  tests  and  less  than  twenty  for  malfunctions. 

If  the  assumption  is  made  that  their  diameter  is  10  microns,  the  settling 
velocity  in  air  will  be  about  0.  02  feet  per  second  at  a  particle  density  of 
2.  0.  This  will  mean  that  the  particles  will  drop  about  seven  feet  per 
mile  traveled  downwind.  Any  error  introduced  into  a  calculation  of  X 
by  neglecting  the  settling  velocity  of  particles  less  than  10  microns  will 
be  much  less  than  the  probable  error  of  the  calculation  induced  by  other 
variables.  Therefore,  in  the  calculation  of  toxic  concentrations  where 
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the  particles  are  sufficiently  small  the  gaseous  diffusion  equations  may 
be  used.  This  can  be  shown  by  the  following  example. 

If  a  detonation  occurs  in  which  1650  it  of  toxic  material  is  released  up  to 
an  effective  altitude  of  12,400  ft.  and  a  calculation  is  made  of  X  versus 
distance,  for  the  assumptions  of  a  gaseous  release  and  also  for  particles 
of  ten  microns  diameter  the  resulting  curves  are  found  (Figure  5)  .  The 
maximum  difference  between  the  two  curves  is  about  forty  percent,  -'nd 
this  error  will  decrease  on  either  side  of  the  maxima  for  X.  If  the 
particles  are  less  than  10  p.  the  particle  curve  will  approach  the  gas  curve. 
The  significance  of  this  is  that  the  toxic  concentration  can  be  calculated 
based  on  the  diffusion  of  gases  instead  of  the  sedimentation  of  solids. 

For  small  particle  sizes,  the  error  introduced  is  small  using  tl'is 
method. 


Maximum  ground  concentration  resulting  from  various  sources  of  toxic 
gas  or  particulates  in  the  atmosphere  is  usually  of  primary  concern.  By 
differentiating  the  above  expressions  and  solving  for  maximum  concen¬ 
tration  at  ground  level,  the  following  expressions  result: 


Instantaneous  point  source  equation  from  Sutton's  equation  for  the 
nonisotropic  case: 


max 
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Continuous  point  source  equation  from  Sutton's  equation  for  an  elevated 
point  source: 
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For  isotropic  conditions: 


Continuous,  infinite  crosswind  line  source  equation  from  Sutton's  equation 
for  the  nonisotropic  case: 
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These  equations  may  be  utilized  for  calculating  the  distance  from  the 
source  to  the  point  of  maximum  concentration  and  for  the  calcu¬ 

lation  of  maximum  concentration  ( 


If  the  Gaussian  concentration  distribution  is  integrated  with  respect  to 
time,  the  resulting  expression  (Reference  4)  describes  the  total  inte¬ 
grated  dosage  (TID): 


TID 


-  h' 
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(  u  t) 


2  -n  i 


TT  u  (tr  t)^ 


TID 


max 


2  Q 


max  dosage 


-A 

IX  e  -  — 

Vc 


TT  e  u  h"^ 

:/(2  -n) 


The  maximum  total  integrated  dosage  (TID^^g^,^)  and  the  distance  from 

the  source  to  the  point  of  maximum  dosage  dosage^  described 

by  these  equations  are  useful  tools  for  determining  toxicity  problems 
downwind  of  a  test  facility. 


So-called  fumigation  conditions  may  develop  when  the  accumulation  of 
material  in  very  stable  air  layers  is  brought  to  ground  levels  by  surface 
turbulence  during  the  initial  stages  of  stable  layer  break-up  or  by  strong 
turbulence  below  relatively  stable  layers.  When  effluent  is  limited  in 
vertical  dispersion  by  a  stable  atmosphere  or  inversion  condition,  the 
center  of  the  plume  which  contains  the  highest  concentration  may  be  forced 
down  near  ground  level.  The  following  formula  is  used  to  express  the  fumi¬ 
gation  concentration  when  the  material  is  confined  within  a  height,  H. 


X,  ,  _  w 

(downwind) - ;;; - 

Cy  3  HX 

The  equations  presented  provide  fairly  accurate  results  provided  the 
proper  diffusion  coefficient  can  be  determined  for  the  particular  situation. 
Additional  experimentation  and  observation  is  needed  to  define  the 
conditions  possible  for  an  infinite  number  of  surface  features  and  climatic 
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conditions.  Precise  predictions  can  probably  never  be  attained,  but 
a  confidence  level  of  predictability  may  be  built  up  through  years 
of  experimentation.  Dispersion  coefficients  are  discussed  in  a  separate 
section  of  this  report  along  with  a  discussion  of  the  reliability  of  the 
coefficients  and  their  effect  on  the  over-all  reliability  of  predicted 
concentrations. 


Continuous  point  source  dispersion  formulae  of  the  Gaussian  distribution 
type  and  the  Sutton  form  are  essentially  equivalent.  Sutton's  dispersion 
coefficient  can  be  expressed  as  a  simple  function  of  distribution  proba¬ 
bilities.  Recent  work  accomplished  for  study  of  toxic  propellant  vapor 
diffusion  under  the  name  of  "Ocean  Breeze"  (Reference  5)  has  resulted 
in  statistical  determination  of  a  purely  empirical  formula  for  continuous 
ground  source  diffusion  with  essentially  zero  buoyancy.  The  "Ocean 
Breeze"  composite  formula  is  expressed  as  follows: 


Q  / 

log  -^  =  0.  1309  -  1.  855  log  x  +  2.  273  log  (ATT  5)-0.  7261  log  a(A) 


Or: 


S' 


1.  352  Q  (ATT  5) 


2,  273 


1.  855  ,  .  ,  0.  7261 

X  a(A) 


Where: 


Cp /=  Maximum  concentration  along  plume  axis,  gm/m^  at 
distance  x,  meters 


Q  =  Source  strength,gm/ sec 


0  (A)  =  Standard  deviation  of  wind  direction,  degrees 


Sutton's  equation  for  continuous  point  source  axial  concentration  with 
reflection  from  ground  surface  is: 


X 


2  Q 

IT  Cy  U  X^ 


n 


C  /  and  X  have  the  same  meaning  and  can  be  expressed  in  the  same  terms. 
P 

Therefore,  the  Sutton  and  "Ocean  Breeze"  equations  can  be  equated  as 
follows: 

2  Q  1.  352  (ATT  5)  ’  Q 


^  C  C  u  x^"'^ 
y  z 


X 


1.  855 


a(A) 


0.  7261 
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The  value  of  n  in  Sutton's  equation  can  vary  between  0  and  1  and 
normally  falls  within  a  range  of  0.  10  to  0.  50.  Where  n  =  0.  145, 

X  =  855^  values  4  a  (A)®'  ^^61  representative 

of  horizontal  diffusion  perpendicular  to  the  plume  axis  in  the  "Ocean 
Breeze"  and  Sutton  equation  respectively. 


Then  the  above  equation  reduces  to  expressions  of  vertical  dispersion. 

- - -  ^  1,  352  (AT'+  5)^‘ 

2  ir  u 

Or:  P  _  _ 1 _ 

2TTu(i.  352)  (AT'+  5)^- 


To  convert  from  the  "Ocean  Breeze"  diffusion  formula  to  Sutton's 
diffusion  formula,  continuous  point  source,  the  following  values  may 
be  equated; 


X 


n  ^1.855 


Cy  -  4  a(A) 


0.  7261 


C. 


2,  704  ^  u  (ATA  5)^' 


This  analysis  should  help  to  illustrate  the  compatibility  of  the  two 
equations  and  suggestsaway  to  apply  diffusion  parameters  and  variations 
of  one  expression  with  respect  to  the  other.  By  analogy,  the  "Ocean 
Breeze"  equation  can  be  modified  for  continuous  elevated  point  source 
applications.  The  following  expression  results: 


C  , 
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1.  352  Q  (AT't  5) 


2.  273 


^.1.  855 


a(A) 


0.  7261 


But: 
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=  C  .  C 

y  z 


2o  (A) 
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1.  352  TT  u  (AT+  5) 
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x2  -  n  _  j^l.  855 
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Therefore-  the  "Ocean  Breeze"  formula  with  stack  height  effect  added 

becomes  /  o  97 

I-  h^)  1-  352^  u(AT'+  5r- 


C‘  1  352  {AT'+  5) 

P  - - — 


2,  273, 


2  a  (A) 


0.  7261 


/ 


Q  ,,1,855 


x*'  a(A) 


0.  726  1 


For  maximum  ground  level  concentration  for  instantaneous  elevated  point 
source  using  terms  available  under  the  "Ocean  Breeze"  experiments: 
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For  total  integrated  dosage  at  the  point  of  maximum  concentration  downwind 
of  an  instantaneous  elevated  point  source,  using  terms  available  under  the 
"Ocean  breeze"  experiments 


TID 


2  Q 


max 


d 


max  dosage 


TT 

e  u  h^ 

1 

/hf^ 

1=4=-') 

/  2  27  3^ 

1,  352  TT  G  h^(ATT  5) 

1.  855 

[c^l 

i  7  ^261 

y  2  a(A)  j 

These  formulae  should  provide  a  suitable  method  for  predicting  downwind 
maximum  concentrations  and  distances  for  relatively  low  stack  height, 
where  the  diffusion  constants  are  not  greatly  different  from  those  determined 
by  instruments  on  network  meteorological  towers.  The  maximum  point  of 
total  integrated  dosage  and  the  distance  to  this  point  are  expressed  in 
standard  diffusion  and  "Ocean  Breeze"  parameters. 

The  foregoing  analytical  treatments  of  diffusion  theory  are  equivalent  to  the 
extent  of  availability  of  experimental  data.  The  chart  on  the  following  page 
(Figure  6)  illustrates  the  differences  in  results  utilizing  two  of  the  equations. 
As  more  appropriate  experimental  data  become  available,  the  curves  will 
more  nearly  coincide.  These  data  are  for  a  continuous  emission  source. 


Due  to  the  extremely  short  duration  of  rocket  tests,  continuous  emission 
equations  are  not  reliable.  Rocket  dest  durations  are  more  nearly 
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Comparison  of  Sutton's  and 
Ocean  Breeze  Equations  for 
Continuous  Ground  Level  Emission 

Assumed  Parameters 

Sutton  Ocean  Breeze 

n  =  0.  Z5  a(A)  =  13° 
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Figure  6.  Atmospheric  Diffusion 
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comparable  to  instantaneous  point  source  emissions.  The  chart  on  the 
following  page  illustrates  this  difference  (Figure  7). 

Also,  on  the  same  chart,  a  curve  of  concentration  at  various  downwind 
distances  of  an  assumed  source  with  a  stack  height  of  100  meters  is  shown. 
It  may  be  seen  that  stacks  will  prove  highly  beneficial  in  the  analysis 
of  meteorology  problems. 
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1.  10  mph  Wind 

2.  =  C  =  C  =  0,  22 

3.  n  =  0.^25 

4.  Q  =0.  5#toxic/sec  for  10  sec  (Rocket  Engine) 


The  continuous  source  curve  assumes  0,  5^/sec 
continuously,  while  the  instantaneous  source  curves 
assume  a  source  strength  of  5  pounds  released 
instantaneously. 


DISTANCE  FROM  SOURCE,  METERS 


Figure  7.  Atmospheric  Diffusion  Sutton's  Equation 


VARIABLE  PARAMETERS  IN  ATMOSPHERIC  DIFFUSION  THEORIES 


The  frictional  drag  of  the  earth's  surface  on  moving  air  results  in  a  decrease 
of  the  wind  velocity  at  the  surface.  Energy  from  the  moving  wind  aloft  is 
rapidly  brought  to  the  earth  by  eddy  currents.  This  process  results  in  a 
change  of  location  of  material  within  the  air  mass  or  in  a  mixing  of  clean 
air  with  effluent  material  or  both.  The  eddy  motion  produces  diffusion  of 
the  material  in  the  air  mass  and  is  caused  by  mechanical  disturbance, 
convection  or  thermal  currents,  or  both.  Diffusion  in  the  polluting  cloud 
will  begin  almost  immediately  after  material  is  released  into  the 
atmosphere  due  to  these  eddy  currents. 

The  diffusive  capacity  of  an  air  mass  is  influenced  also  by  its  vertical 
thermal  structure.  The  temperature  profile,  together  with  atmospheric 
buoyancy  effects,  is  a  necessary  item  for  consideration  when  calculating 
downwind  pollution  concentrations.  (The  rate  of  decrease  of  temperature 
with  altitude  of  an  atmospheric  condition  is  referred  to  as  its  lapse  rate,  ) 
Also  to  be  considered  is  the  degree  of  windiness,  since  high  winds  favor 
rapid  diffusion. 

In  the  most  widely  used  calculations,  Sutton  (Reference  1)  and  modified 
Sutton  type  equations,  the  diffusion  values  are  broken  into  three  distinct 
groupings  which  are  dependent  upon  a  thermal  stability  parameter,  n. 

The  diffusion  coefficients  are  taken  as  C^,  Cy,  and  C^.,  which  are 
measured  in  the  horizontal  downwind,  horizontal  crosswind,  and  vertical 
directions  respectively.  The  dimensionless  number,  n,  is  described  as 
approaching  zero  under  extreme  turbulence  and  tends  to  approach  unity 
as  the  turbulence  decreases.  For  low  altitudes,  less  than  100  meters, 
the  empirical  values  of  the  diffusion  coefficients  and  stability  parameters 
have  been  used  under  certain  meterological  conditions  to  give  successful 
predictions  of  downwind  concentrations  for  distances  of  a  few  kilometers. 

Any  diffusion  due  to  varying  wind  velocities  in  the  atmosphere  giving  rise 
to  shear  forces  is  ignored.  Sutton  states  that  the  diffusion  coefficients, 

Cy  and  C^,  depend  on  the  value  of  n  and  the  magnitude  of  surface  gustiness. 
The  diffusion  coefficients  decrease  with  increased  height  due  to  a  normal 
steady  decrease  of  turbulence  and  tend  to  approach  each  other  as  the 
mechanical  turbulence  decreases  with  increased  height. 

The  Sutton  form  for  ground  level  concentration  for  a  continuous  point 
source  emission  is  given  as: 
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Where : 


X  =  Downwind  concentration 

u  =  Mean  wind  velocity 

X  =  Distance  downwind 

Q  =  Weight  emission  rate  per  unit  time 

y  =  Lateral  distance  from  vertical  plane  through  source 

and  parallel  to  mean  wind  direction 

h  =  Stack  or  cloud  height,  meters 

The  evaluation  of  shell  bursts  in  the  air  has  been  used  by  Sutton  to  obtain 
the  decrease  in  diffusion  coefficients  with  gain  in  height.  The  decrease  in 
value  has  been  found  to  follow  the  empirical  form  (Reference  1): 

C  =  C  (o)  -  0.  075  log^^  Z 

where  C(o)  is  the  ground  level  diffusion  coefficient  and  Z  is  the  height 
measured  in  meters. 

The  values  of  C  and  are  assumed  to  be  equal  for  heights  above  Z5 
meters  and  Cy  =C^  =  C.  The  value  of  n  is  obtained  through  its  relationship 
to  the  mean  wind  velocity  profile  in  the  equation  (Reference  1). 


n 


Whe  re  : 


u, 


Average  wind  velocity  at  elevation  Z 
Average  wind  velocity  at  reference  level  Z^ 


The  diffusion  coefficients  can  be  obtained  from  the  relationship  (Reference  1). 

/  \  i-n 
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(l-n)(2-n)U^  \ 


Where; 


V  =  Kinematic  viscosity 

V  =  Average  wind  variation 
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n  -  Stability  parameter 
u  -  Average  wind  velocity 

In  the  case  of  rough  surfaces,  the  kinematic  viscosity  is  replaced  by 
n;  icroviscosity,  N. 

Experimental  determinations  by  Holland  have  shown  the  value  of  Cy 
for  various  stability  parameters  to  be  greater  than  postulated  by 
Sutton  ,  and  the  values  of  n  to  generally  be  larger  than  Sutton's 
value  for  approximately  the  same  defined  stability  conditions.  Holland 
also  obtained  data  showing  that  the  diffusion  coefficients  decrease 
with  increasing  wind  speed.  The  discrepancies  betv'^een  measured 
and  theoretical  values  are  attributed  to  time  differences  in  sampling, 
wind  speeds,  and  terrain  The  values  of  Sutton's  parameters  and 
the  experimentally  obtained  values  by  Holland  are  shown  in  the 
following  table  : 

Values  from  Sutton  for  C  and  C-,  (na  ^  )  follow'; 

y  ' 


Stability  Condition _ Elevation  (Meters) 


Large  Lapse 

25 

SO 

75 

100 

n  r  20 

0.  21 

0.  17 

0.  16 

0.  12 

Small  Temperature  Gradient 

n  -■  ,25 

0.  12 

0  10 

0.  09 

0.  07 

Moderate  Inversion 

n  ^  .  33 

0  08 

0  06 

0.  05 

0  04 

Large  Inversion 

n  -  .50 

0.  06 

0.  05 

0.  04 

0.  04 
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Oak  Ridge  Experimental  Diffusion  Parameters  by  Holland  are  presented 
below: 


Stability  Condition 

Paramete  r 

Elevation 

(Meters) 

5  5 

16.  5 

47 

Z4.  4 

Daytime  (Superad- 
iabatic  Lapse  Rate) 

C  hiij^  ) 

V  2 

0.  34 

0  Z8 

0  ZZ 

0,  Zl 

n  -  1 5  1 0  20 

y  ^ 

(.  'mil) 

^  s 

0  Z3 

Li !  m  /  s  e  c ) 

1  5 

1  8 

Z  4 

Z  7 

Average  'Neutral 

Stabil  ity ) 

C,.  (m^l 

0  zo 

0  17 

0  1  5 

0  13 

n  -  Z  5 

C^(mf ) 

0  Zl 

u  (m/ sec) 

0,  O'J 

1  Z 

1,  8 

Z,  7 

Night  (Moderate 
ln\-ersion) 

Cy  'mf) 

o'l:  ' 

0^  14 

0  11 

0.  08 

0.  08 

n  --  .  30  to  40 

u  14 

u(m/ sec) 

0.  3 

0  6 

1  5 

Z  7 

On  the  basis  of  the  Round  Hill  (Reference  Zl)  and  Project  Prairie  Grass 
'References  ZZ  and  Z3l  experiments,  a  different  type  diffusion  equation 
has  been  developed.  The  values  of  Cy  and  are  replaced  by  standard 
deviations  of  gas  concentration  distributions  at  distance  x.  The  form  of 
the  equation  is- 
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Where  X,  Q,U)  and  y  are  the  same  quantities  defined  in  Sutton's 
equation,  and  z  is  the  distance  measured  from  a  horizontal  plane 
through  the  source  and  the  various  diffusion  parameters  are: 


V  Z  Oy  - 
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it  IS  well  to  note  that  in  the  Gaussian  dispersion  equation  the  value  of 
n  IS  replaced  by  n  and  n  respectively,  since  recent  work  indicates  a 
\’ariation  does  exist  for  lateral  and  vertical  diffusion  at  low  altitudes.  At 
hiyh  .ititudes  the  assumption  that  n^  :  n^  "  n  seems  reasonable  but  cannot 
be  supported  due  to  lack  of  experimental  data. 

From  experimental  da’:a  obtained  at  Edwards  Air  Force  Base,  the 
Stanford  Research  .institute  computed  a  series  of  diffusion  coefficients 
for  tha'  particular  hiyh  desert  area  (Ref.  24).  Using  the  instantaneous 
poin’  source  equanon  postulated  by  Sutton,  the  Stanford  Research 
j’nsti'ute  back  calcinated  the  coefficients  by  obtaining  concentrations  of 
tr.'C';r  materials  downwind  of  a  tracer  source.  The  Sutton  equation  for 
jioinr  source  and  ground  level  emission  is: 


X 


Q 


3  [i-^)/7: 


and  llie  i-aiues  computed  by  Stanford  Research  Institute  are: 


Stabilit y  C-ondirion 

- 

Cx 

Cy 

Midday  .  Turbulent  .• 

0.  20 

0.  23 

0.  51 

0.  23 

Hea'uJy  Oi'ercast 

0  25 

0.  16 

0.  38 

0.  16 

inversion 

0.  33 

0.  07 

0.  10 

0.  07 

Since  the  tracer  was  released  at  a  relatively  low  altitude,  the  diffusion 
coefficients  approximately  coincide  with  Sutton's  values. 


The  value  of  the  diffusion  coefficient,  Cy,  as  obtained  by  the  Stanford 
Research  .institute  a^  Edwards  Rocket  Site,  is  gi\'en  as  0.  51  for  midday 
lapse  rates.  The  value  was  obtained  by  substituting  into  the  Sutton 
formula  for  plume  width: 

2 


When : 


2  |ln 


lod'i 

Pj 


C, 


,3-n/2 


Vq  I  Plume  width  from  centerline,  meters 


p  t  Per  cent  of  peak  concentration  at  distance  from  centerline 
Cy  Lateral  diffusion  coefficient 
X  -  Distance  to  point  of  determination,  meters 
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n  =  Stability  parameter 


Usi/ig  the  information  collected  from  a  series  of  diffusion  tests,  an 
average  plume  width  was  obtained  and  the  above  equation  solved  for  C 


y’ 


As  may  be  surmised,  the  conditions  under  which  the  average  plume 
distribution  curve  concentrations  were  obtain .  d  varied  for  each  test. 
Accurate  measurements  of  the  distance  to  the  arc  of  sampling  stations 
were  not  taken  but  were  between  six  and  seven  miles 


From  the  graphical  data  presented  in  the  Stanford  Research  institute 
report,  a  more  accurate  determination  of  versus  p  can  be  obtained, 
and  this  will  give  a  revised  value  of  C,. 

The  new  value  of  is  equal  to  2900  meters  rather  than  3220  meters;  a 
new  value  of  Cy  is  found  to  be  0,4b  This  value  is  approximately  what  is 
found  by  interpolation  of  Sutton's  values  to  one  foot  above  ground  level. 
Interpolation  to  tin-  height  of  3  feet  by  Sutton's  values  gives  0.  395  for  Cy 
which  is  the  approximate  height  at  which  the  samples  were  collected. 

The  difference  between  the  Cy  values  can  be  attributed  to  terrain  roughness 
lack  of  accurate  measurement  of  the  sampling  arc,  and  plume  width 
measurements 

The  values  for  Cy  recalculated  on  the  basis  of  six  and  seven  miles  in 
lieu  of  6,  5  miles  give  0  49  and  0.  43  respectively,  or  a  possible  error  of 
6  per  cent  on  this  variable  alone.,  It  would  then  be  assumed  that  Sutton's 
values  of  Cy  are  in  good  agreement  for  the  Rocket  Site. 

Logarithmic  plots  of  Sutton's  values  have  been  extrapolated  to  heights  of 
over  1,  000  feet,  as  shown  in  Figure  8,  in  an  attempt  to  obtain  approximate 
diffusion  coefficients  at  high  altitudes,  but  the  values  obtained  have  not 
been  checked  with  actual  experimental  operating  conditions. 

The  discrepancies  between  Sutton  s  diffusion  coefficients  and  Holland  s 
experimentally  derived  values  are  pointed  out  in  Figures  9  and  10  for 
low  and  medium  altitudes.  Extrapolation  of  the  values  given  on  Page  29 
to  great  heights  was  used  to  illustrate  the  converging  characteristics  of 
the  data.  The  Sutton  values  were  modified  where  applicable  by  the 
recommended  correction  factor  given  by  Sutton  for  n=  0.  20  and  corrected 
by  the  Wanta  equation  for  n  -  0  2  5 

Figure  11  shows  dispersion  coefficients  and  stability  parameters  advanced 
by  Sutton  and  determined  by  Holland  'Ref,  2).  By  taking  the  average 
values,  fairly  accurate  predictions  can  be  made  Data  taken  at  the  test 
site  will  be  helpful  in  further  improvement  of  the  values  and  will  lead 
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HORIZONTAL  DISPERSION  COEFFICIENT  C^ 

Figure  11.  Recommended  Diffusion  Coefficients 

for  Sutton  s  Equations 


to  values  applicable  to  Edwards,  Surface  diffusion  coefficients  determined 
for  the  Edwards  Rocket  Site  (Ref.  24)  agree  reasonably  well  with  surface 
values  indicated  by  this  graph.  The  decrease  in  dispersion  coefficients 
with  altitude  has  been  applied  *‘o  elevated  source  and  stack  height  effects 
throughout  this  report.  The  reduced  values  result  in  increased  distances 
to  maximum  points  of  concentration  and  T.lD,  Use  of  these  dispersion 
values  is  recommended  until  such  time  as  improvement  can  be  made 
through  observations  and  interpretation  of  experimental  releases  of  non¬ 
toxic  plumes  and  those  from  actual  toxic  tests. 


Variances  in  relating  Ocean  Breeze-type  equations  to  peak  concentrations 

downwind  of  an  emitting  source  are  a  definite  factor  to  consider.  A  major 

factor  which  must  be  taken  into  consideration  when  applying  Ocean 

Breeze-type  equations  for  obtaining  downwind  pollutant  concentration  is 

plume  meandering,  A  correction  factor  of  this  type  does  not  seem  to 

appear  in  the  original  report.  From  wind  trace  photos  and  actual  plume 

photos  of  meandering  smoke  plumes,  a  pattern  strongly  suggesting  a 

sinusoidal  wave  appears.  The  Ocean  Breeze  repoit  mentions  normalized 

peak  concentration,  ( ^ ^ ^ ^ .  which  was  obtained  by 
^  source  strength  rate  ’ 

experimentation  as  the  concentration  downwind  for  varying  parameters. 

It  is  assumed  that  the  concentrations  obtained  during  the  experiment  are 

actually  average  concentrations  for  the  period  of  tracer  emission;  hence 

the  actual  peak  concentration  has  a  somewhat  greater  value  than  that 

obtained  by  Ocean  Breeze  experiments. 


If  the  assumption  is  correct  that  the  pattern  of  the  plume  cross  section  is 
sinusoidal,  the  peak  value  wou.d  be  1/0,  636  times  the  average  peak 
concentration,  i.  e.  ,  1.  57  times  average. 

Comparison  of  the  values  obtained  using  the  Ocean  Breeze  equation  and 
Sutton's  equation  for  a  continuous  emitting  source  and  multiplying  a 
correction  factor,  1.57,  shows  the  results  to  be  quite  close  for  a  point 
5.  000  meters  downwind.  Of  course,  through  judicious  choice  of  Sutton's 
diffusion  coefficients,  it  should  be  possible  to  obtain  coinciding  values  for 
the  two  equations,  A  mathematical  model,  which  would  be  based  on 
probability  curves,  could  be  set  up  to  take  this  meandering  effect  into 
account.  A  probability  function  is  necessary  in  view  of  the  nearly 
infinite  combination  of  variables  which  affect  the  plume  dispersion,  such 
as  diffusion  coefficients,  wind  velocity,  rate  of  direction  change,  and 
height  of  emitting  source.  Based  on  these  assumptions,  the  peak  con¬ 
centration  of  contaminant  downwind  of  a  long  time  emitting  source,  such 
as  a  rocket  test  lasting  for  several  minutes,  will  be  greater  than  cal¬ 
culated  by  Ocean  Breeze  equations. 
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TOTAL  INTEGRATED  DOSAGE 


Integration  of  ground  level  concentration  for  instantaneous  point  source 
releases,  with  respect  to  time,  produces  an  equation  that  can  be  used  to 
calculate  the  concentration-time  exposure  for  atmospheric  release  of  toxic 
products  downwind  of  the  source.  Differentiation  with  respect  to  distance 
from  the  source  yields  an  expression  from  which  the  maximum  dosage  can 
be  found  and  distance  determined.  The  equations  (Ref.  Z)  have  been  pre¬ 
sented  under  the  discussion  of  diffusion  formula.e  and  are  repeated  here 
for  convenience. 


TID  = 


2  Q 


2  /  -  4.\2  -n 


u  t) 


TID 


u(u  t)^‘  ^ 
2  Q 


max 


"^ax  dosage  2 


7T  e  u  h^ 

V(2  -  n) 
h^' 


C‘ 


Using  these  equations,  the  downwind  dosage  can  be  plotted  with  respect  to 
distance.  The  TID^^g^^^  expression  defines  the  maximum  point  on  the  TID 

curve  at  distance  d^i^x  dosage. 

The  two  most  significant  toxic  emission  parameters  are  the  maximum  con¬ 
centration  and  the  total  time  of  exposure.  With  some  toxic  products,  the 
limits  are  expressed  only  as  maximum  allowable,  but  with  others  the 
integrated  dosage  is  apparently  of  prime  significance. 

Integrated  dosage  will  be  important  when  the  maximum  allowable  concen¬ 
tration  (MAC)  varies  with  the  interval  of  exposure.  This  occurs  when  the 
allowable  exposure  for  personnel  is  greater  per  eight-hour  work  day  than  a 
continuous  exposure  for  the  general  public.  The  computation  of  TID  is  done 
by  multiplying  the  period  allowable  concentration  by  the  time  period  to 
obtain  the  dimensions  of  mass  time  per  volume  for  TID.  The  allowable  TID 
varies  between  inhabited  areas  and  in-plant  short  time,  daily,  or  weekly 
doses.  Therefore,  TID  would  have  to  be  applied  to  the  specific  occurrence 
as  required.  The  allowable  TID  for  a  definite  period  would  be  found  by 
muliplying  the  MAC  for  the  period  by  the  duration  of  the  period.  The  units 
obtained  could  be  gram  seconds  per  cubic  meter. 

The  standard  methods  of  reducing  TID  are  followed  for  toxic  rockets  ;  that 
is,  reduction  of  toxic  effluent  rate  by  scrubbing  or  filtration,  increasing  the 
stack  height  via  afterburning  and/or  flaring,  and  scheduling  for  optimum 
atmospheric  conditions. 


38 


effective  stack  height  formulae 


A  key  factor  in  considering  the  fallout  from  a  toxic  rocket  test  is  the 
maximum  altitude  reached  by  the  exhaust  plume.  This  is  because  the 
dispersion  of  the  toxic  material  effectively  begins  at  the  exhaust  plume 
maximum  altitude,  and  the  toxic  material  disperses  in  three  dimensions. 
Therefore,  a  relatively  small  portion  reaches  the  ground. 

The  factors  governing  the  maximum  altitude  or  effective  stack  height  are 
a  combination  of  the  following  parameters: 

A.  Test  rocket  size 

1.  Engine  thrust 

2.  Duration 

B.  Propellant  characteristics 

1.  Exhaust  temperature 

2.  Mixture  ratio 

3.  Average  molecular  weight  of  exhaust 

C.  Rocket  orientation 

1.  Vertical  upward 

2.  Horizontal 

D.  Atmospheric  conditions 

1.  Wind 

a.  Velocity 

b.  Direction 

c.  Time  variations 

d.  Altitude  variations 

2.  Diurnal  changes 

a.  Ambient  temperature 
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b.  Temperature  variation  with  altitude  (lapse  rate) 

(1)  Adiabatic  lapse  condition 

(2)  Inversion  (negative  lapse  rate) 

c.  Convective  circulation 

From  these  many  parameters,  a  large  number  of  distinct  formulae  have 
been  derived,  utilizing  data  from  wind  tunnel  tests  of  smoke  models,  full 
scale  power  plant  observations,  mathematical  derivations,  and  actual 
rocket  firings.  Very  few  have  a  common  resemblance,  as  they  differ  in 
specific  variables  treated,  general  form,  and  method  of  approach.  This 
will  not  produce  any  contradictory  information  in  the  ultimate  results  if 
the  proper  equation  is  used  for  the  particular  application  .  Therefore,  the 
best  of  the  available  equations  for  the  present  situation  may  be  selected 
on  its  particular  .applicability  or  by  a  process  of  elimination.  Analternate 
method  would  be  the  derivation  of  a  completely  new  formula.  The  new 
formula  would  be  designed  to  cover  the  range  of  stack  heights  required  for 
rocket  usage.  Derivations  of  two  new  formulae  are  presented.  One  formula 
is  intended  for  use  in  the  near  future.  The  other  formula  depends  on  an 
accurate  determination. of  constants  and  coefficients;  however,  the  use  of 
this  equation  will  eventually  provide  highly  accurate  results. 


The  following  survey  of  the  existing  stack  height  formulae  is  presented. 
The  sources,  accuracy,  and  general  applications  are  presented  in  a  later 
section  of  the  report. 


This  formula  (Ref.  2,  3)  was  derived  from  wind  tunnel  experiments  and 


gives  the  height  to  where  the  centerline  of  the  plume  is  near  horizontal. 

This  is  not  to  be  misconstrued  as  the  maximum  stack  height.  Factors 
which  are  not  considered  are  the  stability  of  the  atmosphere  (i.  e.  ,  lapse 

rate)  and  the  velocity  scaling  ratiq  ^  ,  which  varies  due  to  compressibility 

u 

effects  as  the  stack  gas  ejection  speed  (Vg)  becomes  large.  These  factors 
would  indicate  optimistic  plume  rise  and  a  high  degree  of  plume  dilution  as  the 
velocity  ratio  increases.  As  the  above  equation  is  presented,  there  is  no 
limit  to  the  effect  that  a  large  velocity  ratio  gives  to  the  calculated  height. 
This  equation  also  does  not  include  the  effect  of  energy  dissipation  at  high 
rocket  gas  exhaust  temperatures  which  would  decrease  the  temperature 


ratio,  AT 


Therefore,  this  equation  is  not  recommended  for  use  in 


analysis  of  rocket  test  systems. 
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2,  Bosanquet,  Carey,  and  HaJton  equation 


The  equation  iRef,  6)  gi\-en  by  Bosanquet,  et  al  ,  for  the  maximum  plume 
rise  in  n  stable  atmosphere  is; 

f 


_ 4,  7_7 _ iQv,  VJ 
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The  Bosanquet,  ei  a!  equation  is  based  on  Theoretical  developments  and 
employs  some  fundamental  experimental  constants  of  diffusion. 


This  equation  ln:S  been  tested  using  chimneys  and  similar  low  heat  sources. 
A  report  by  Moses  and  .Strom  'discussed  below  under  the  Scorer  formula), 
stated  that  a  \  alue  of  stack  height  calculated  from  this  equation  gave  a 
minimum  a'  crage  d'fference  between  the  caiciuated  and  observed  heights. 

It  would  seem  that  tins  equation  would  co'er  most  of  the  parameters,  but 
certain  important  ones  are  left  out  These  are  entrainment  rate,  which 
would  modify  'he  '  elocity  rate  at  something  other  than  a  constant  rate, 
and  the  effect  of  extreme  exhaust  temperatures  on  the  temperature  ratio 
and  buoyancy  of  the  gases,  in  order  to  be  useful  for  the  purpose  of  rocket 
firings,  this  equation  would  require  modification  and  testing. 
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Oak  Ridge  formula 


Ah 


1  5  Vg  d  +  0  000134  Qh 
ri 


Where  ■ 

is  heat  emission  rate  of  s'ack  gas  relative  to  ambient 
air,  calories  per  second 

and  u  are  in  meters  per  second 

Ah  and  d  are  in  feet 

The  Oak  Ridge  formula  'Ref  3)  was  denied  empirically  from  observations 
of  hot  plume  rise  from  steam  generating  plants,  assuming  the  velocity 
effect  as  denied  b\  Rupp,  et  al  The  purpose  of  the  Oak  Ridge  equation 
was  to  impro'.'e  agreement  between  predicted  and  observed  results.  This 
formula  does  no<  mve  maximum  plume  Iieighi,  since  it  was  derived  from 
data  taken  only  a  few  siack  heights  downwind. 
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The  experimental  data  used  for  deriving  this  empirical  equation  did  not  take 
into  consideration  variations  in  lapse  rate.  Holland  suggests  from  10  to 
Z5  %  should  be  added  to  Ah  for  large  lapse  rates  and  a  like  amount  sub¬ 
tracted  for  small  lapse  conditions  Various  investigators  found  that  the 
Oak  Ridge  formula  gives  consistently  lovv  values  for  stack  height.  There¬ 
fore,  a  safe  assumption  may  be  made  that  any  equation  which  gives  a 
lower  numerical  answer  that  the  above  equation  should  be  disregarded 
for  the  application  at  EAFB  rocket  site 

The  effluent  from  the  X-10  reactor  cooling  air  exhaust  at  Oak  Ridge  and 
two  other  hot  exhaust  systems  were  used  to  empirically  derive  the  buoyancy 
term  of  this  equation  The  disappointing  factor  is  that  the  observed  value 
of  the  stack  height  is  not  the  maximum  since  observations  were  made 
only  a  short  distance  dovnv-ind  from  the  stack  This  equation  is,  therefore, 
not  recommended  since  it  docs  not  give  the  maximum  stack  height. 

4.  Thomas  formula 


Ah 


1  5  V's  d  r  0  000268  Qj^ 


The  units  in  tin?  equation  are  the  same  as  the  Oak  Ridge  formula.  The 
Thomas  equation  (Ref  7)  is  uscti  to  derive  the  maximum  plume  rise. 

The  possible  variation  due  to  changes  in  stability  is  ±10  to  25%as  mentioned 
in  the  Oak  Ridge  formula  This  method  is  considered  by  Stern  (Ref.  3)  to 
give  a  fairly  accurate  value  for  stack  height  Whether  it  can  be  extrapolated 
to  the  degree  required  for  large  rocket  firings  must  be  determined  by 
experiment . 

5.  Sutton  for mula 


Ah 

Where: 


And:  Cp :  -  Cpg 

The  distance  along  the  plume  path  (s)  is; 


When  ifis  less  than  10°,  the  effective  stack  height  approaches  the 
maximum,  and  when  the  wind  speed  is  high,  s  equals  the  downvind 
distance  from  the  stack  (x) 


[  cot  t  cosec  v 


log(,  (  cot  'A-  cosec  i!/  )  ] 
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Sutton's  formula  only  takes  buoyancy  into  consideration  and  does  not 
include  momentum  rise  (Ref.  1).  One  should  refer  to  his  book  prior 
to  using  his  equations,  since  the  value  of  Cj^  is  not  clear  as  to  the  units 
which  should  be  used.  However,  for  dimensional  continuity,  must 
be  dimensionless.  In  the  determination  of  the  effective  stack  height, 
the  value  of  t  chosen  is  quite  important  since,  as  the  angle  between  the 
plume  axis  and  horizontal  {\|f)  approaches  zero,  the  function  of  iji  goes 
through  a  maximum  and  then  becomes  discontinuous  as  a  limit  of  zero 
degrees  is  reached. 

The  results  from  this  method  are  generally  too  high  to  be  of  use,  and 
are  considered  to  be  too  high  in  most  of  the  referenced  material. 

6.  Bryant  and  Cowdrey  experiment 

This  procedure  (Ref.  8)  is  based  on  wind  tunnel  tests  and  on  observations 
of  stack  effluents.  A  series  of  graphs  has  been  formulated  but  has  to 
be  greatly  extrapolated  and/or  extended  to  obtain  a  result  for  the 
purpose  of  this  application. 

One  factor  apparent  from  the  graphs  is  the  decrease  in  effective  stack 
height  with  increase  in  turbulence.  Another  important  factor  is  that  of 
condensation  from  the  cooling  plume.  As  the  gas  rises  and  cools,  any 
water  formed  from  the  combustion  process  will  condense  when  the 
temperature  decreases  below  the  dew  point.  Therefore,  if  the  rising 
plume  does  not  mix  with  the  surrounding  air,  the  water  will  tend  to 
rain  out  and  carry  any  solids  or  water  soluble  gases  with  it.  However, 
under  actual  field  tests,  Bryant  and  Cowdrey  state  that  effluent  dilutions 
of  the  order  of  magnitude  of  several  hundreds  to  one  occur.  This  then 
would  be  sufficient  to  unsaturate  any  rocket  exhaust  plume. 

In  conclusion,  it  is  believed  that  this  procedure  is  of  the  correct  general 
type,  but  the  conversion  to  rocket  testing  would  involve  extrapolations  of 
such  magnitudes  that  the  accuracy  of  information  would  be  invalidated. 

7.  Priestley  procedure 

Priestley  gives  a  two-phase  solution  to  stack  height  (Ref.  9).  The  first  is 
that  of  a  vertical  jet  moving  through  a  medium.  The  motion  of  the  jet 
causes  resistance,  spreading,  and  dilution  due  to  eddy  turbulence.  As 
the  jet  is  decelerated  by  the  action  of  the  air,  the  intensity  of  this  jet 
energy  becomes  less  than  that  of  normal  air  turbulence.  At  this  point, 
the  second  phase  of  dilution  begins;  one  that  is  dependent  on  atmospheric 
properties. 
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By  this  method,  Priestley  derived  the  following  relations; 


Maximum  plume  rise 


Ah  =  Ahi  +  g  A  0i' 


G  ' 


In  these  equations,  subscripts  of  "o"  refer  to  the  chimney  top,  and 
subscripts  of  "1"  refer  to  conditions'  at  transition  from  phase  1  to 
phase  2.  In  phase  1,  the  axial  speed  of  ’  the  plume  is  given  as  w^  . 


2  K  ■  Ah/y  Ah/ 

Values  of  w^^  between  "o"  and  "  may  be  found,  w^^  being  the 
necessary  one.  Other  relations  are  given  below; 


TT  0  r 
S  PS 


This  is  a  buoyancy  term.  Ah^  is  measured  from  a  point  located  a 
distance  Ah^^  below  the  chimney  top  where  the  lines  of  plume  expansion 

form  a  virtual  point  of  origin.  Radius  R  of  the  plume  cross  section  is 
related  to  Ah^. 

R  =  c  Ah,, 


Also  at  point  "o"  this  relation  hold  true. 


R  -  c  Ah 

0  VO 


Priestly  suggested  that  an  alternate  form  for  the  determination  of  R^ 
should  be  used. 


w  A0 
ao  ^°o 


The  value  of  the  spreading  coefficient  (c)  in  the  above  equations  is 
based  on  smoke  plume  photographs  in  wind  tunnel  experiments. 

(  .  \  0.  5 


=  0.  75 


2.  44 
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Where  u  is  wind  speed  in  meters  per  second  and  c  is  dimensionless. 
The  value  of  Ahy^  is  then  defined  as  : 

.6  3  Ai  g  , 

vi  -  —  1  "ao  ‘-‘“vo - -  1=  0 


I^Ah 


VI 


Ah 


2  e'c® 


I 
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2  e' 


The  mixing  rate  (k)  is  given  as  proportional  to  wind  speed  and  is  in  units 
of  inverse  seconds  when  u  is  in  meters  per  second. 


k  =  0.  0197  u 


The  value  of  Ah^  may  be  replaced  with  (Ah^  ^  -  Ah^^)  when  solving  for 
Ahj,  .  In  the  procedure,  the  time  to  reach  the  transition  point  can  also 
be  found.  This  is  given  below: 


t  = 


^  P  A.  g 

2  Aig  V  2 


2 


(Ahv  -  Ah^o  )  +  ]  .  [w^Q  Ah^o  ] 


The  value  of  Wjj_  is  now  found,  and  then  the  potential  temperature  is 
found  at  "1". 


A0/=  - - 

c^  Ahvi^w^  1 


Finally,  the  value  of  Ah^  can  be  found.  This  method  gives  the  lowest 
calculated  stack  height  of  any  of  the  formulae.  It  is  of  a  form  that 
should  give  "order  of  magnitude"  results  but  apparently  does  not.  A 
simpler  empirical  relation  is  more  to  be  desired. 


8.  Bosanquet  formulae 


Ah  =  Ag  u 


Where: 
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= 


a  = 


t  +  t. 


Aa 


+ 


t  -  — 


The  above  procedure  (Ref.  3,  10,  11)  is  evaluated  using  consistent  units, 
such  as  the  English  ones  used  in  previous  equations  ;  t  is  time  and  (a) 
and  (a^)  are  graphically  solved  on  page  142  (Ref.  3).  For  values  out¬ 
side  the  range  of  the  graphs,  the  following  apply: 

Where  (a)  is  very  large: 

fi  (a)  =  logg  a  -  0.  12 

When  (a)  is  very  small; 

fi  (a)  =  1.  054  a^ 


When  (a^)  is  very  large: 

:  1,  311  a^t  - 


^11  (^0^ 

When  (a^)  is  very  small; 

fu  (ao)  =  -  0.  527  a^‘ 


-  1 


The  value  of  Cg  is  0.  13  and  t  +  t^  <200  sec.  For  momentum  rise  only: 


Ah 
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2  Q 
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Cg  TT  U 
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When  Vg  /u  >0.  5 


For  small  values  of  V,  /u 
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The  two  equations  for  Ah^^^^^^  give  the  same  value  at  Vg/u  =  0.  48. 
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The  maximum  rise  of  a  buoyant  plume  in  a  stable  atmosphere  is  found 
for  t  +  tp  of  200  sec.  ,  or  the  following,  whichever  is  smaller; 


t  "I"  t 


o 


1.  527 


i 


The  accuracy  of  this  method  has  been  stated  by  Bosanquet  to  be  about 

±21  %. 


Bosanquet  derived  this  procedure  by  assuming  that  when  a  cloud  of  hot 
gas  is  rising,  the  total  heat  content  and  total  upward  momentum  are 
unaffected  by  dilution  with  atmospheric  air.  It  is  also  assumed  that  the 
upward  momentum  increases  at  a  rate  proportional  to  the  heat  content. 

The  rate  of  dilution  is  assumed  to  be  proportional  to  the  surface  area  of 
the  cloud  multiplied  by  a  function  of  the  wind  velocity  and  the  velocity  of 
the  cloud  relative  to  the  surrounding  atmosphere.  This  procedure  gives 
a  low  value  for  stack  height  when  com.pared  with  other  methods. 

It  can  be  noted  that  the  calculation  of  the  time  or  rise  is  a  function  of 
lapse  rate.  Since  the  lajise  rate  varies  with  altitude,  any  calculation  of 
time  must  weigh  this  factor. 

The  statement  by  Bosanquet  that  the  time  of  plume  rise  is  equal  to  or 
less  than  200  seconds  seems  to  imply  low  powered  systems  since  time 
periods  greater  than  this  have  been  observed  for  large  explosions.  For 
example,  up  to  ten  minutes  for  atomic  bombs  has  been  reported  by  the 
AEC  (Reference  17). 

The  velocity  ratio  (Vg/u)  and  the  buoyancy  or  temperature  ratio  (AT^/Ti) 
have  been  considered  by  Bosanquet,  but  perhaps  not  enough  for  the  purpose 
of  this  report.  This  is  because  the  compressibility  effect  of  high  velocity 
exhaust  streams  becomes  noticeable  and  the  large  buoyancies  due  to 
temperature  differences  also  become  important. 

The  dilution  coefficient  (C^)  has  been  evaluated  at  0,  13,  which  provides 
an  adequate  fit  to  the  calculated  stack  heights.  Other  values  have  been 
offered,  0.  09  and  0.  17  for  example,  which  give  less  accurate  answers. 

The  Bosanquet  formula  is  not  recommended  due  to  the  consistent  low  values 
of  stack  height  obtained  by  this  method.  The  equation  is  not  recommended 
until  further  trials  at  stack  heights  more  pertinent  to  rocket  tests  have  been 
made. 
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9  Scorer  formu'lc'i 
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Scorer  assigns  three  phases  to  the  rise  of  a  plume.  Phase  one  has  the 
plume  nearly  vertical  and  the  upward  velocity  exceeds  the  horizontal 
wind  Velocity.  This  phase  ends  when  the  vertical  velocity  approximately 
equals  the  horizontal  wind  velocity.  Phase  two  is  a  transition  stage  in 
which  the  plume  bends  over  and  become  horizontal.  This  phase  ends 
when  the  vertical  velocity  equals  the  turbulent  edrly  velocity  of  the 
atmosphere.  This  occurs  when  the  upward  velocity  is  about  10%  of  the 
wind  velocity,  according  to  Scorer,  Phase  three  is  the  regime  of 
diffusion  at  the  effective  stack  height.  This  is  accomplished  by  the  Sutton 
and  other  diffusion  formulae.  Scorer's  method  is  derived  from  the 
observations  of  chimney  plumes.  Upon  application,  it  appears  as  accurate 
as  other  methods.  Moses  and  Strom  (Ref.  12)  have  compared  observed 
plume  rises  with  values  obtained  from  various  formulae. 

They  calculated  the  stack  height  from  (1)  Holland  (Oak  Ridge),  (2)  Bryant 
and  Davidson,  (3)  Sutton,  (4)  Scorer,  (5)  Bosanquet,  Carey,  and  Halton, 
and  (6)  Bosanquet,  and  compared  the  results.  Their  conclusions  state, 
"There  is  no  one  formula  which  is  outstanding  in  all  respects.  "  They  do 
not  recommend  the  Sutton  formula,  but  consider  the  others  about  equally 
valid  for  the  low  energy  level  conditions  tested.  Whether  the  Scorer 
method  can  be  extrapolated  to  rocket  testing  is  not  known,  but  it  does 
seem  to  give  far  too  optimistic  results.  Scorer's  assumptions  for  the 
values  of  and  (i  are  not  clearly  explained,  and  these  constants  may 
vary  with  the  system  being  considered. 


10.  Machta  formula 


Ah 


M 


- - -  •  ^og 

M  ■  a  z 
mass  of  cloud. 


aM 
a  z 
ae' 


a  z 

kilograms 


ae' 

\ 

(A0')^  +  a  z 

1  a  M 

M  ’  a  z 

/ 

z  =  upward  pointing  coordinate,  meters 

9^  =  temperature,  degrees  Kelvin 

(A  6^)^=  potential  temperature  difference  between  cloud  and 
environment,  degrees  Kelvin 


—  .  =  0.  5  •  10  cm  (cumulus  clouds) 

M  a  z 
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This  equation  (Ref.  2,  13)  is  for  an  instantaneous  source,  an  atomic 
bomb,  and  gives  the  height  at  which  a  cloud  is  no  longer  warmer  than 
its  surroundings.  Close  agreement  was  obtained  upon  comparison 
with  actual  bombs. 

The  factor  which  tends  to  give  optimistic  heights  from  this  formula 
is  the  entrainment  rate.  Cloud  size  is  an  important  factor,  the  larger 
the  cloud  the  smaller  the  rate,  with  the  given  value  being  for  cumulus 
clouds.  In  the  case  of  rocket  testing,  an  order  of  magnitude  increase 
is  given  as  a  function  of  cloud  radius  to  rocket  cloud  radius.  Any 
further  improvement  in  this  method  would  be  subject  to  actual  tests. 
This  equation  could  possibly  be  used  after  an  explosive  malfunction 
of  a  rocket  engine. 


1 1.  Empirical  formula  for  small  solid  rockets 


Ah  =  3.  4  (  1.  2  -  n)^ 


0.  5 


Where: 

n  is  a  stability  parameter  equal  to  Sutton's  n  so  that  n  =  0.  25 
for  adiabatic  lapse  rate. 

=  propellant  weight,  pounds 

Tg  '  -  gas  effluent  temperature,  degrees  Centigrade 
u  =  wind  velocity,  miles  per  hour 
Ah  =  cloud  height,  feet 


This  equation  (Ref  14)  is  for  small  solid  rockets  fired  horizontally, 
and  does  not  contain  upward  momentum.  It  was  developed  empirically 
during  actual  testing  and  is  not  necessarily  valid  upon  extrapolation  to 
high  thrust  levels.  It  is  not  known  if  the  calculated  height  from  this 
equation  is  a  maximum. 
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12.  Jet  Penetration  Equation 


A  wind  tunnel  investigation  was  made  under  NACA  direction  into  the 
relationship  the  momentum  of  gases  from  an  orifice  or  stack  has  to 
the  altitude  achieved  by  the  gases.  (Ref.  26) 

Based  on  experimental  data,  an  empirical  formula  was  derived  far  the 
penetration  of  an  air  stream  by  a  jet  discharging  through  an  orifice 
perpendicular  to  the  air  stream.  The  equation  for  the  penetration  was 
obtained  in  terms  of  the  jet  diameter,  distance  downstream  from  the 
jet,  the  ratios  of  the  jet,  and  air  stream  velocities  and  densities: 

5  _ 

=  2.  91 

Pa  Va  V  Dg 

depth  of  jet  penetration  (stack  height  due  to  momentum) 

distance  downstream  from  orifice  centerline 

orifice  diameter  (stack  diameter) 

mass  density  of  jet  air  at  vena  contracta 

mass  density  of  free  air 

velocity  of  jet  at  vena  contracta 

velocity  of  free  stream  air 

Use  of  the  above  equation,  with  proper  substitution  of  terms  related  to 
the  exhaust  gas  stream,  gives  a  stack  height  of  120  meters, for  a  333  #  /sec 
engine  of  40  %  gaseous  exhaust  at  4000  F  and  13.  2  psia  atmospheric 
pressure  with  a  molecular  weight  of  29  discharging  through  a  2  meter  stack 
at  150  meters  per  second  initial  velocity  ,  at  a  distance  of  6.  2  miles.  As 
the  distance  downstream  is  increased  to  62  miles  the  height  due  to  velocity 
is  240  meters  and  at  400  miles  is  only  425  meters.  Although  no  term 
appears  for  different  lapse  rates,  it  is  believed  that  only  a  very  small  error 
is  introduced  by  the  omission  since  the  height  due  strictly  to  exhaust  velocity 
is  small  in  comparison  to  the  buoyancy  effect  of  the  hot  gas. 

This  formula  has  merit  for  obtaining  stack  heights,  since  it  is  based  on 
actual,  physical,  controlled  experimental  data. 
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A  further  simplification  of  the  formula  is: 


Ah., 


1.  92  D 


0.  70 


0.'606 


0.  606 

X 


0.  30 


A  major  drawback  to  this  formula  is  that  the  information  was  gathered  at 
extremely  large  free  stream  air  velocities,  260  and  360  feet  per  second. 
However,  smooth  curves  are  shown  in  the  original  report  for  the  Penetration 

Coefficient,  il  ,  versus  the  density-velocity  ratio,  Ps 

D„  - TT"  ;  from  these 

Pa  ^a 

data,  the  assumption  can  be  made  that  the  ratio  should  hold  for  small  free 
stream  air  velocities.  The  variation  in  penetration  coefficient  was  found 
to  be  negligible  for  Reynolds  numbers  between  60,  000  and  500,  000  for  the 
jet  and  for  viscosity  ratios  of  1.5  to  1.9. 


13.  Buoyancy  of  Explosion  Clouds 


The  buoy-incy  term  for  stack  height  has  been  shown  to  be  a  function  of  stack 
diameters  rat io  of  the  kinematic  viscosities  lapse  rate,  temperature 
■  fference  flow  rate,  wind  speed,  and  stack  velocity.  This  term  is 
I  "esented  below  as  the  buoyancy  equation  derived  for  stack  height  due 
to  buoyancy  from  a  continuous  source  The  buoyancy  equation  was 
dei'ived  for  norm.il  runs  For  explosions  a  sornewh'it  different  view 
must  be  taken  Instead  of  stack  diameter  flow  rate,  and  stack  velocity, 

'i  I  miss  of  explosive  must  be  considered  s  well  as  the  other  terms  men- 
t).jiied  'ibove  Some  of  these  tf^rms  will  generatp  d  relatively  large  change 
m  licigh.t  of  th,e  plume  with  a  small  change  in  tlieir  magnitude.  The  exact 
degree  to  which,  this  ch.ange  will  be  import-mi  is  not  presently  known,  but 
’!,e  tot.:l  'unction  has  an  experimentally  determined  maximum  from  nuclear 
explos  ions 

The  '.".■entual  height  reached  by  a  nuclear  explosion  cloud  depends  on  the 
I  e-it  energy  of  the  bomb  -ind  on  the  temperature  gradient  (lapse  rate)  and 
density  of  the  surrounding  air.  The  greater  the  explosive  yield,  and 
therefore  tlm-  iieal  reledsed  the  greater  will  be  the  upward  thrust  due  to 
buoy.Tncy  and  therefore  a  higher  st.i  ck  h  eight ,  Fac'ors  which  influence 
tlie  cioucl  heiglit  other  than  those  given  above  are  of  importance  in  con¬ 
sidering  rocket  test  explosions  as  well  as  nuclcdr  explosions.  These 
f.jctors  include  the  height  of  the  explosive  device  at  the  instant  of  detona¬ 
tion  since  any  large  distance  between  tlie  fireball  and  the  ground  will  be 
<dded  to  tl'.e  total  height  readied  by  the  cloud  wlicreas  if  the  fireball 
touches  the  ground  a  large  amount  of  debris  will  be  vaporized  and  picked 
up  Witt  the  oth.er  buoyant  material  and  cause  a  net  loss  in  cloud  height. 

Tlie  m-jximum  height  achieved  by  the  nuclear  cloud  is  strongly  influenced 
by  the  tropopause  the  boundary  between  the  troposphere  and  the  strato¬ 
sphere  assuming  the  cloud  reaches  the  height  of  the  tropopause.  The 
tropop-ause  is  between  the  tropospliere  and  the  rcl-atively  stable  air  of  the 
isothermal  stratosphere  It  varies  in  altitude  with  seasons  and  latitude, 
ranging  from  23  000  feet  ncnr  the  poles  to  about  55,  000  feet  in  equatorial 
regions  When  -i  nuclear  cloud  reaches  the  tropopause  there  is  a  tendency 
to  spread  lateially.  This  effect  is  caused  by  the  relative  magnitude  of  the 
diffusion  coefficients  changing  at  the  altitude-  the  horizontal  ones  being 
greater  fhr,n  the  vertical  ones,  A  portion  of  the  cloud  will  penetrate  the 
tropopmse,  and  it  will  ascend  into  the  more  stable  air  of  the  stratosphere 
provided  sufficient  energy  remains  in  the  cloud  at  this  height. 

The  nuclear  detonations  carried  out  by  the  United  Stales  h.ave  been  heavily 
instrumented  for  scientific  purposes.  The  properties  of  interest  for  this 
report  -ire  the  explosive  yield  of  a  test  tlic  height  of  the  burst,  type  of 
burst  (air  drop  tower,  surface  or  underground),  the  mean  height  of  the 
cloud  riiid  1 1- e  height  of  the  tropopause  From  this  information  for  many 
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different  magnitude  explosions,  a  graph  has  been  constructed  of  the  thermal 
yield  versus  the  effective  height  of  the  cloud  in  feet.  The  US  Atomic  Energy 
Commission  has  stated  that  the  equivalent  of  one  pound  of  TNT  is  900  B  T  U, 
Therefore,  if  the  heat  release  of  an  explosion  is  known,  the  effective  height 
of  the  resulting  cloud  can  be  determined. 

From  the  following  graph  (Figure  12)  the  height  that  the  explosion  cloud 
will  reach  is  obtained.  A  complete  detonation  of  a  600  ff  solid  propellant 
grain  will  give  an  effective  height  of  4,  700  feet  (assuming  50%  of  the  total 
heat  is  availaljlc). 

Also,  a  complete  detonation  of  a  1  1,300  #  grain  will  give  an  effective  height 
of  8,  500  feet  (assuming  50%  of  the  total  heat  is  available).  The  basis  for 
this  height  calculation  is  that  the  total  enthalpy  of  the  solid  propellant  is 
4346,  calories  per  gram. 

The  probability  of  a  detonation  of  such  large  solid  propellant  grains  is 
lov/  but  it  definitely  is  present.  A  test  stand  designed  for  toxic  rockets 
must  provide  for  such  an  eventuality.  This  can  be  done  by  completely 
containing  any  explosion,  or  by  creating  an  open  structure  which  will  give 
the  maximum  upward  impetus  to  the  toxic  materials.  The  open  structure 
would  be  constructed  with  sloping  sides  (upward  and  outward)  so  that  the 
maximum  fraction  of  the  explosive  forces  will  be  directed  upward.  Any 
vaporizing  or  powdering  of  the  structure  should  be  avoided  since  the 
explosive  energy  would  be  attenuated  by  this  action. 

A  final  recommendation  for  the  use  of  the  graph  would  be  to  calculate  the 
explosive  yield  and  the  minimum  safe  effective  stack  height  for  the  rocket 
to  be  tested.  If  the  height  given  by  the  graph  is  too  low,  then  either  the 
test  should  be  forbidden  or  a  means  of  augmenting  the  explosive  stack 
height  should  be  found.  One  possible,  but  not  too  feasible,  way  to  accom¬ 
plish  this  would  be  to  place  a  sufficient  quantity  of  shock  sensitive  non- 
toxic  explosive  adjacent  to  the  toxic  rocket. 
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uoyancy  Rise  of  Explosion  Clouds 


14.  Dimensional  Analysis  Formula  for  Effective  Stack  Heights 


The  effective  stack  height  is  a  function  of  many  variables  which  are 
both  dependent  and  independent.  By  the  use  of  dimensional  analysis 
these  variables  may  be  related  to  each  other  and  to  the  effective  stack 
height.  Dimensional  analysis  is  useful  since  it  assists  in  the  determi¬ 
nation  of  convenient  arrangements  of  variables  in  a  physical  relation 
and  helps  plan  systematic  experiments.  The  first  step  is  to  list  all 
the  variables  involved.  This  list  is  a  result  of  experience  and  judgment 
as  to  the  variables  of  importance.  After  the  list  has  been  made,  the 
formal  procedure  of  analysis  is  used  to  relate  the  variables. 

The  stack  height  equation  can  be  written  symbollically  as; 

Ah  =  f(G,  Vs,  u,  d_e  ^  ^  ^  T,  AT, 

d  z  d  z 

S>  .  Cp  ,  d,  p,  M,  u) 

The  value  of  Ah  will  be  computed  in  two  distinct  parts;  those  of  upward 
momentum  and  of  buoyancy.  That  is: 


Ah  -  Ah  +  Ah. 

m  D 

Ah  f  (Momentum) 

m 

Ah|^  -  f  (Buoyancy) 

The  assumption  is  made  that  all  of  the  momentum  is  dissipated  before 
buoyancy  is  added. 


By  use  of  the  Buckingham  tt  Theorem,  a  versatile  mathematical  tool, 
the  two  expressions  for  the  plume  rise  were  obtained.  The  momentum 
equation  is; 


Ah 

m 


{  K  '  d^  q 


(Frf 


This  expression  is  applicable  to  jet  flow  systems  and  also  possibly  to 
instantaneous  upward  releases  or  puffs.  For  a  continuous  flow  the 
mass  (M)  would  equal  G  . 
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where: 


1  ,,  2 
-  P 

2  s 

=  Dynamic  pressure 

of  stream 

Re 

d  VgP 

=  Reynolds'  Number 

of  stream 

— 

M- 

o 

Fr 

V 

s 

=  Froude's  Number  : 

_  Inertia  Force 

d  g 

Gravity  Force 

The  plume  is  assumed  to  lose  no  buoyancy  during  the  time  .  Therefore, 

an  expression  for  Ah,  would  start  from  Ah  and  continues  up  to  Ah,^  ,  , 

^  b  m  ^  Total 


Therefore; 


Ah^^ 


where 


V 


r  gf  d^  Av"^ 


„  -c'  -d' 

a,  “ 


P 


kinematic  viscosity 


Therefore,  it  is  shown  that  the  buoyancy  rise  is  a  function  of  the  force  of 
gravity,  stack  diameter,  gas  viscosity,  lapse  rate,  gas  temperature, 
volume  flow  rate,  wind  speed,  and  gas  velocity  Finally,  from  the  above 
derivations,  the  equation  for  effective  stack  height  is  written  by  combining 
the  two  equations,  one  of  momentum  and  one  of  buoyancy.  The  equation  for 
effective  stack  height  must  be  evaluated  in  terms  of  real  numbers  for  K^, 

K"',  and  the  various  exponents  before  it  can  be  of  any  practical  use.  A 
prog'"am  to  accomplish  this  is  possible  The  parameters  to  be  measured 
are  part  of  the  equation.  A  series  of  tests  in  which  the  various  known 
quantities  are  changed  would  allow  an  evaluation  of  the  constants  and 
exponents.  Such  a  series,  at  low  energy  levels,  has  been  conducted. 


Using  the  given  data,  a  series  of  equations  were  set  up  to  find  the  unknown 
constants  and  exponents.  As  is  shown  in  the  Appendix,  the  evaluation 
of  the  problem  leaves  an  equation  with  twelve  unknowns. 


For  a  complete  evaluation  of  the  exponents  and  constants,  this  equation 
must  be  evaluated  a  total  of  twelve  (12)  times  using  twelve  different  runs. 
One  factor  which  will  assist  in  evaluating  the  buoyancy  term  is  the  necessity 
of  maintaining  dimensional  consistency.  The  unit  equation  is: 


(m) 


(m  ) 


Therefore  these  relations  must  hold  true; 


1  =  f  +  k  -  (a')  +  3  (-c')  +  (  -d')  +  (  -  e') 
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and 


0 


2f +  (  -  c')  +  (  -d')  +  (  -e') 


The  above  analysis  indicates  the  degree  of  complexity  involved  in  the 
calculation  of  effective  stack  heights.  This  holds  true  even  when  all 
constants  and  exponents  are  known,  and  when  there  are  twelve  unknowns 
the  problem  becomes  one  for  a  computer  program. 
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15.  Buoyancy  Equation  for  Continuous  Source 


An  expression  for  the  rise  due  to  buoyancy  has  been  derived,  based  on 
the  basic  properties  of  the  atmosphere  and  the  laws  of  thermodynamics. 
The  final  form  consists  of  a  relation  between  buoyancy  rise,  heat 
emission  rate,  atmospheric  stability  parameters,  and  wind  velocity. 

The  equation  is  presented  below;  the  entire  analysis  is  in  the  Appendix. 


Ah,,  =  0.  050  Qh 

0.463 

(  u) 


for  standard  lapse  (n  =0.25) 


The  above  expression  assumes  100  per  cent  thermodynamic  efficiency 
and  is  based  on  standard  lapse  conditions.  The  lapse  rate  will  affect 
the  work  done  against  gravity  in  lifting  the  diffused  air.  For  dry  adiabatic 
lapse  rates,  the  work  done  would  be  zero.  However,  a  strong  lapse  will 
be  limited  in  height, or  infinite  buoyancy  would  result.  The  strong  lapse 
will  result  in  an  appreciable  increase  in  buoyancy  rise  over  isothermal 
or  inversion  conditions.  The  overall  lapse  rate  effect  from  limited 
plume  rise  data  is  shown  to  be  the  following  power  of  to  provide 
lapse  rate  adjustment.  Assuming  80  per  cent  thermal  efficiency,  the 
equation  takes  the  following  form: 

2  -  n 

Ah^=0.04(Qh)  ^  (u) 


Where  Ah^  =  maximum  rise  due  to  buoyancy,  meters 

=  heat  release  above  atmospheric,  BTU/sec 
H  =  mean  wind  velocity  to  maximum  rise,  meters/sec 
n  =  stability  parameter  of  Sutton  where: 

n  =  0. 20,  strong  lapse 
n  =  0.  25,  standard  lapse 
n  =  0.  33,  inversion 
n  =  0. 50,  strong  inversion 
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The  above  formula  is  the  buoyancy  effect  rise  of  a  hot  plume  from  a 
continuous  point  source.  This  formula  will  at  least  give  a  theoretical 
maximum  buoyancy  height.  When  any  other  factors  which  might  influence 
the  buoyancy  rise,  are  included  in  the  expression,  they  might  lower  the 
calculated  height.  All  indications  are  that  the  height  will  not  be  greater 
than  the  value  obtained  from  this  expression  for  standard  lapse  conditions. 
When  strong  lapse  or  inversion  conditions  are  present,  the  equation  will 
only  give  approximations,  the  answers  being  directed  in  the  proper 
direction.  As  is  recommended  with  any  theoretical  expression,  this 
equation  should  be  tested  under  a  wide  range  of  conditions  before  complete 
confidence  is  placed  in  it. 

Considering  only  buoyancy  due  to  heat,  a  comparison  of  several  formulae 
has  been  made  (Figure  13).  The  formulae  shown  are  those  considered  to 
be  the  most  useful  for  rocket  engine  testing.  The  graph  plots  plume  rise 
in  meters  versus  heat  release  rate.  The  Thomas  formula  is  considered 
conservative  but  is  based  primarily  on  a  heat  release  rate  of  3500  BTU. 
Only  the  buoyancy  portion  of  the  Thomas  formula  is  plotted  on  the  graph. 

The  thermally  derived  equation  based  on  published  data  provides  a  reason¬ 
able  fit  for  much  of  the  published  data  on  observed  plumes  where  the 
molecular  weight  of  the  cloud  is  essentially  that  of  air.  It  is  recommended 
that  the  buoyancy  equation  for  continuous  source  be  used  until  proven  or 
modified  by  test  observation.  The  wind  velocities  assumed  on  Figure  13 
are  estimated  averages  for  the  various  lapse  conditions.  Also  shown  is 
the  altitude  attained  by  the  horizontal  firing  of  small  rockets  from  the 
empirical  equation.  It  is  well  to  note  that  no  consideration  is  gi\'en  to  the 
rate  at  which  the  propellant  burns.  For  this  reason.,  the  use  of  this 
formula  is  not  recommended  for  application  to  larger  engines  where  high 
rates  of  heat  release  are  obtained,  Also  since  some  energy  must  be 
expended  in  diffusion  in  the  horizontal  plume  with  further  loss  due  to  the 
momentum  energy  expended  in  the  horizontal  direction,  much  lower  plume 
heights  would  be  expected  than  with  conventional  formulae. 

A  presentation  of  the  most  useful  inversion  penetration  equation  is  also 
shown  to  indicate  how  easily  formulae  break  down  at  low  heat  release 
rates . 
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igure  ^  Buoyancy  Rise  frorn  Continuous  Sovirces  of  Hot  Gas 
from  Recommended  Equation 


General  Discussion  of  Stack  Formulae 


From  the  above  descriptions  of  the  available  formulae,  the  conclusion 
is  that  there  is  no  perfect'  equation  for  the  prediction  of  stack  height  over 
an  infinite  range.  This  is  true  whether  the  reference  is  to  rocket  firing, 
steam  plants,  atomic  reactors,  or  industrial  waste  discharges.  There  is 
no  equation  that  will  give  complete  accuracy  outside  the  region  from  whence 
it  was  developed.  Various  investigators  have  attempted  to  show  the  regions 
of  application  for  certain  formulae.  Their  findings  have  been  that  each 
equation  will  provide  results  with  a  certain  magnitude  of  error.  Part  of 
this  error  is  due  to  defects  in  the  formulae  and  part  is  due  to  inherent 
properties  of  the  atmosphere. 

The  equations  of  Sutton,  Bosanquet,  et  al.  ,  and  the  Oak  Ridge  formula 
were  studied  by  A.  C.  Best  (Ref.  15)  who  calculated  the  stack  heights 
obtained  under  varioas  wind  velocities,  and  strength  of  heat  source.  His 
conclusions  were  ".  .  .  .  that  if  the  stack  designer  is  concerned  with  the 
maximum  ground  level  gas  (or  toxic)  concentrations  ever  likely  to  be 
experienced,  he  will  arrive  at  similar  answers,  probably  within  the  limits 
of  the  variability  of  human  reaction  to  specified  concentrations,  whichever 
formula  he  chooses.  "  Best  states  chat  Sutton's  equation  gives  only 
approximate  answers.  Bosanquet,  et  al.  ,  has  been  shown  to  give 
reasonable  agreement  with  some  full-scale  observations,  and  the  Oak 
Ridge  gives  "a  reasonable  approximation  to  the  results  given  by  the  other 
two,  "  How  Best  can  come  to  these  conclusions  when  he  has  calculated 
differences  of  o\'er  ten  times  for  some  of  the  examples  of  stack  heights  is 
not  known. 

Upon  discussion  of  the  article  by  Best  before  the  institute  of  Fuel  in 
London,  various  dissident  voices  were  raised.  Dr.  R,  S.  Scorer 
(Equation  9)  stated  that  many  phenomena  of  great  importance  were  com¬ 
pletely  ignored  in  Best’s  presentation.  The  most  obvious  were  the 
results  of  variations  of  height  with  changes  in  the  stability  of  the  atmosphere. 
Also,  a  Mr.  Lucas  stated  that  he  thought  Sutton's  formula  underesti¬ 
mated  the  plume  rise,  and  that  the  Oak  Ridge  formula  underestimated  it 
by  such  a  large  factor  that  he  considered  it  to  be  practically  worthless. 

J.  E,  Hawkins  and  G.  Nonhebel,  in  their  article  on  "Chimneys  and  the 
Dispersal  of  Smoke"  (Ref.  16),  state  that  the  Oak  Ridge  formula  cannot  be 
extrapolated  over  large  ranges,  but  the  Thomas  equation  is  much  better. 

Upon  comparison  of  the  Bosanquet  formula  with  actual  plume  rise,  an 
accuracy  of  i  25%  was  found  in  the  chimneys  observed. 

As  was  mentioned  above,  under  the  Scorer  formula,  the  work  of  Moses 
and  Strom  indicates  that  there  is  no  one  formula  which  is  outstanding 
in  all  respects. 
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Some  of  the  formulae  are  of  an  extremely  simple  form  (i.  e.  ,  The  Oak 
Ridge  formula)  and  probably  do  not  cover  all  the  necessary  parameters 
to  insure  accuracy.  In  order  to  insure  a  reasonable  approximation, 
they  must  include  these  parameters  in  empirically  determined  constants. 
These  constants  will  apply  more  clcseJv  in  regions  near  which  they  were 
determined. 

The  decision  as  to  which  stack  height  equat  on  to  use  depends  upon  the 
type  of  effluent  system  present.  In  the  case  of  explosive  detonations  the 
graph  of  explosive  heat  release  versus  effective  height  represents  the 
most  advanced  information  available.  This  graph  is  therefore  recom¬ 
mended  for  use  with  detonations- 

Where  the  rocket  is  oriented  to  fire  hcrizontallv,  and  a  flame  deflector 
is  not  used  to  direct  the  exhaust  vertically,  the  equation  for  buoyant  gases 
applies.  This  equation  was  derived  from  laws  of  thermodynamics  and  is 
estimated  to  give  answers  to  1  f  or  2  places  of  accuracy- 

The  directing  of  the  exhaust  out  of  a  vert'cal  stick  will  nrovide  a  momentum 
rise  component  to  the  gases.  This  rise  i  >  on -- -ch  r  e  H  r-dH u ive  lo '.he  buoyanoy 
rise  and  thus  a  summation  of  the  two  is  made  to  ob’'.am  the  total  effective 
stack  height.  The  momentum  rise  equation  chosen  is  based  cn  jet  penetration 
studies  in  wind  tunnels  and  is  believed  to  give  reasonable  results. 

The  effects  of  atmospheric  variation--  are  'mportant.  They  will  introduce 
scatter  in  observed  versus  calculated  heights  and  scatter  v/ill  increase  with 
distance  downwind  from  the  source  of  emission.  Atmospheric  parameters 
of  most  consequence  are  changes  in  lapse  rate  with  altitude,  wind  velocity 
with  altitude,  and  wind  direction.  These  changes  will  cause  eddies  in  any 
direct  path  that  a  plume  might  follow  when  unhindered.  A  forceful  example 
of  the  effects  produced  by  atmospheric  variation  cn  man-made  clouds  would 
be  in  the  differences  in  height  of  nuclear  weapons  tests.  A  graph  of  mean 
cloud  height  versus  explosive  yield  will  show  a  scatter  of  points  for  any  one 
nominal  yield  (Ref  17).  For  example,  in  the  one  hundred-ton  range,  the 
observed  heights  vary  between  seven  and  fifteen  thousand  feet  above  mean 
sea  level.  Also,  in  the  twenty  thousand  ion  r  ange,  the  heights  vary  between 
25  and  35  thousand  feet.  Factors  causing  these  variations  include  the 
atmospheric  stability,  moisture  content.-  and  the  height  of  the  tropopause 
(30,  000  -  55,  000  ft. ) 

All  of  the  stack  height  equations  contain  a  wind  velocity  parameter.  This 
parameter  is  considered  fairly  constant  during  the  calculation  but,  in 
actuality,  has  considerable  variation  with  altitude.  At  ground  level  the 
wind  velocity  may  be  about  ten  to  twenty  feet  per  second  but  as  altitude 
increases,  the  velocity  increases.  The  me?m  wind  velocity  at  the  altitude 
of  maximum  winds  is  generally  between  80  -i-.d  180  fps  (Ref.  18)  and  these 
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upper  winds  may  have  veJocvties  of  ZOO  to  300  fpa  one  per  cent  of  the 
time.  The  altitude  at  which  this  occurs  is  between  30.  000  and  40,  000 
feet,  or  approximately  lha,t  of  the  trcpc-priuse. 

The  mean  wind  direction  of  the  upper  a;r  is  fromi  the  v-esl  by  southwest, 
at  the  latitude  of  Edwards  Air  Force  Base,  at  an  elevation  of  about  35,  000 
feet  The  steadiness  parameter  .'mean  vector  wind/mean  scalar  wind) 
is  approximately  0.8  for  the  mean  wind  at  35,  000  feet.  This  will  enable  a 
prediction  cf  fallout  tc  be  made  with  a  greater  average  confidence  level. 

From  ’he  above  discuss. on,  i>  is  eh''wn  'h-rt  the  exact  height  attained  by 
a  plume  cannot  at  present  be  predicted  wTh  a  high  degree  of  accuracy. 
However,  a  minimum  heght  and  a  p-obable  maximum  height  are  possible 
to  calculate. 

The  maximum  h^^ights  would  be  based  ca  exp'  .s.on  heights  for  an  instan¬ 
taneous  source,  and  the  buoyancy  plus  momentum  equations  for  continuous 
sources. 

The  maximum  concentration  of  contsm.nants  at  ground  level  can  be  found 
if  the  minimum  height  is  known.  The  energy  of  a  plume  may  be  augmented 
by  flaring  if  the  con'.ami:.-*  on  it  excels. ve  but  i.n  order  to  reduce  the  cost 
per  test  caused  by  flaring,  the  stack  height  should  be  known.  A  definite 
program  for  >he  determination  cf  emp.r'cai  formulae  might  be  implemented. 
This  program  would  nc'ude  the  .a’atic  tes'ing  of  rocket  engines  cf  different 
thrust  levels. 

During  test  firings,  information  should  be  recorded  as  *:o  lapse  rates,  wind 
vectors,  exhaust  temperature  and  ccmpcsiticn,  engine  thrust,  test  duration, 
and  plume  height  as  a  function  of  d' stance  downwind  .;-"'m  the  test  stand. 
From  these  data,  the  dimensional  analysis  equation  can  be  further  developed 
to  reduce  the  probable  error  of  calcv.iated  •■■‘•irk  heights  from  other  methods, 
The  probable  error  of  the  refined  data  would  si-c  be  available  for  use,  and 
from  'his  combination  a  method  of  stack  heigh’  prediction  would  be  possible. 

In  the  absence  of  precise  data  and  experience  w  ’h  s’^ack  heights  for  large 
rocke’s  at  the  Edwards  Rocket  S.te,  it  is  'ecommended  that  the  following 
equations  be  used  ■ 

1.  For  detonation  the  relation  between  energy  release  versus 
stack  height  -curve  Fig  12,i. 

2.  For  horizontal  fir'ng  the  buoyancy  equation  equation  15. 

3.  For  vertical  firing  the  buoyancy  equat  on  plus  the  jet 
penetration  equation,  equation  12- 


4.  After  sufficient  data  are  accumulated,  a  final  solution 
of  the  dimensional  analysis  equation  may  be  used. 

5.  During  inversion  conditions,  a  combination  of  equations 
could  be  used  for  cross-checking  of  predicted  results. 
That  IS,  the  buoyancy  and  momentum  equation's  (15  and 
13)  answers  would  be  compared  to  the  answer  from  the 
inversion  penetration  formula  presented  below, 


INVERSION  PENETRATION 


The  atmospheric  temperature  has  a  daily  variation  with  height.  This 
variation  is  caused  by  solar  radiation,  absorption  and  emission  of  heat  by 
clouds  and  the  ground,  and  mixing  of  layers  of  air  by  wind  forces.  When 
the  temperature  of  the  air  decreases  with  an  increase  in  altitude,  this  is 
called  a  lapse  condition  and  the  rate  of  change  is  the  lapse  rate.  The  U.  S. 
Standard  Atmosphere  has  a  lapse  rate  of  3.  5  F/IOOO  feet.  If  the  lapse 
rate  is  zero,  an  isothermal  condition  exists.  When  the  air  temperature 
varies  directly  with  altitude,  a  "negative"  lapse  rate  or  an  inversion 
condition  exists.  An  inversion  is  caused  by  the  rapid  cooling  of  surface  air. 
The  nocturnal  radiation  of  the  ground  without  interference  by  clouds  often 
results  in  inversion  conditions 

An  inversion  is  a  critical  condition  during  a  firing  of  a  toxic  rocket.  The 
hot  exhaust  gases  wall  rise  until  they  reach  an  equilibrium  height  within 
the  inversion  layer  and  will  remain  at  this  level  until  the  inversion  is  dis¬ 
sipated  by  outside  forces  such  as  those  produced  by  sunrise.  This  dissi¬ 
pation  of  the  inversion  will  distribute  the  pollutant  vertically  due  to  the 
turbulent  mixing  which  occurs  during  the  formation  of  normal  lapse 
conditions. 

If  the  rocket  exhaust  gases  have  sufficient  potential  energy,  they  will 
penetrate  an  inversion  This  energy  is  supplied  by  upward  gas  velocity, 
gas  buoyancy,  and  the  large  amount  of  heat  produced  in  the  rocket  combustion 
chamber.  The  inver  sion  height  which  can  be  penetrated  is  mathematically 
definable  as; 

Z  f  (g,  Ta,  p,  Q  ,  D 

Where: 

Z  -  Vertical  distance  to  top  of  inversion 
g  =  Gravitation  constant 
Q  =  Heat  potential  of  exhaust  gases 
p  =  Air  density 
Cp  =  Specific  heat  of  air 

r  =  Average  lapse  rate  from  ground  to  top  of  inversion 
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The  expression  becomes.’ 


Z  = 

or:  Z  = 
where  K  = 
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This  expression  has  previously  been  used  in  rhe  prediction  of  power  plant 
exhaust  plume  paths  (Reference  19  and  20).  In  that  function,  its  form  was: 

3  1 

Z  =  784  (  D  ®  Q. 

n 

r  is  in  F  per  1000  feet, 

Z  is  in  feet. 


Q  is  the  megawatt  rating  of  power  piarr  based  on  nhe  assumption  that 
one  third  of  the  power  goe-s  out  the  stack. 


Converting  Q  to  BTU/sec  available  from  a  rocket  exhaust,  the  expression 
becomes; 


And; 


Z  ^ 

K' 

K'  = 


K'  (D 
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3  times  the  power 


y(948  BTU/sec)  per  megawatt 


1 

T 


186. 


Then; 


_3_ 

Z  =  186  (  D®  Qh  ■ 


This  will  provide  a  formula  that  will  allow  the  scheduling  of  a  toxic  rocket 
test  during  an  inversion.  Given  the  parameters  of  inversion  intensity  and 
height,  the  minimum  thrust  level  can  be  calculated.  Conversely,  the 
strongest  inversion  that  a  given  thrust  level  may  penetrate  can  be  found. 
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As  an  example: 


Given:  lOOK  toxic  solid  rocket  engine 


Total  heat  available  with  afterburner  (approx)  =  1,  5  x  10^ 
BTU/sec 

Strong  inversion  where  F/1000  feet 


Solving  for  Z: 


(186)  (20)  (1,  5  X  10®  )’^  Feet 


2120  Feet, 


Since  diurnal  inversions  of  both  this  strength  and  height  are  extremely 
rare,  the  lOOK  engine  could  pierce  alm.ost  any  inversion. 

There  are  other  available  means  of  overcoming  inversions.  The  simplest 
is  to  avoid  them  by  proper  scheduling  of  tests  for  periods  of  minimum 
inversion  probability.  This  would  mean  firing  during  late  morning  to  early 
afternoon  and  avoiding  the  early  evening  to  early  morning  inversion.  Due 
to  normal  working  hours  and  check  out  procedures,  this  should  not  be  an^ 
inconvenience  to  test  personnel. 

Additional  methods  of  overcoming  inversions  include  the  use  of  tall  stacks, 
heat  implementation  of  the  rocket  source,  and  proper  upward  velocity.  A 
tall  stack  is  the  device  generally  used  at  power  generating  plants  to 
achieve  the  dispersion  of  hot  exhaust  gases.  These  stacks  are  quite  high, 
ranging  to  700  feet  (Reference  3)  and,  therefore,  expensive.  It  is,  of 
course,  evident  that  any  physical  stack  will  help  in  normal  testing  to 
increase  the  effective  stack  height.  The  most  likely  stack  height  is  one 
higher  than  all  other  adjacent  structures,  perhaps  one  of  50  to  100  feet. 

Normally,  the  exhaust  from  the  rocket  will  by  itself  provide  an  after¬ 
burning  effect  due  to  its  fuel  rich  nature,  since  a  large  portion  could  consist 
of  gaseous  hydrogen  and  carbon  monoxide.  These  gases  will  be  emitted 
above  their  ignition  temperature  and  will  spontaneously  ignite.  Heat  imple¬ 
mentation  would  consist  of  injection  of  some  inexpensive  fuel  into  the  high 
temperature  fuel  rich  exhaust  in  the  stack.  After  mixing  with  the  toxic 
exhaust,  the  combined  stream  would  be  ignited  by  a  pilot  flame  and  the 
resulting  torch  would  penetrate  the  inversion. 

The  velocity  of  the  rocket  exhaust  gases  emitting  from  the  stack  has  an 
effect  on  the  entrainment  rate  of  the  plume  in  the  atmosphere.  Excessive 
air  turbulence  could  be  generated  if  the  gas  velocity  is  too  great,  which 
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would  cause  dissipation  of  energy,  preventing  penetration  of  the  inversion 
layer.  The  same  applies  to  a  gas  velocity  which  is  too  low.  A  velocity 
of  less  than  half  that  of  sound  and  perhaps  closer  to  a  few  hundred  miles 
per  hour  would  be  o;  timum. 

As  has  been  stated  above,  inversions  present  a  resistance  to  penetration. 
This  holds  true  in  both  directions  That  iS,  once  the  inversion  has  been 
penetrated,  any  subsequent  fallout  of  toxic  material  will  be  inhibited  by 
the  top  of  the  inversion,  This  condition,  called  lofting,  will  cause  the 
plume  to  diffuse  upward  and  sideways  and  very  little  downward.  This  is 
a  most  desirable  situation  and  perhaps  one  that  might  be  planned  for, 
since  the  inverion  impedes  the  toxic  materia!  from  reaching  the  ground. 

At  the  same  time,  the  toxic  materia!  is  rap.dly  diluted  by  the  normal 
lapse  conditions  above  the  inversion.  It  shoald  be  noted  that  an  inversion 
need  not  disappear  w.th  the  sunrise,  but  couM  persist  for  several  days, 
Atmospheric  conditions  are  not  constant,  and  transient  conditions  may 
occur  at  any  location  and  at  any  time.  The  probabii.ty  of  these  changes 
taking  place  at  any  particular  time  can  be  predicted  for  the  test  site,  but 
for  large  distances  downwind,  no  rigid  prediction  can  be  made  unless 
data  are  taken  at  these  points.  General  weather  forecasts,  location  of 
fronts  and  storms,  and  the  time  of  day  will  assist  in  prediction  of  diffusion 
changes.  When  changes  occur,  averag.ng  of  diffusion  calculations  may  be 
made  to  estimate  the  hazard  of  a  change  from  normal  lapse  to  inversion  or 
washout  conditions. 

For  proper  test  scheduling,  it  would  be  advisable  to  determine  the  strength 
and  height  of  inversions  at  the  test  s.te  The  standard  method  is  the  use 
of  a  radiosonde,  a  balloon-borne  devic-a.-which  trar.rmits  temperature  and 
pressure  to  a  ground  station.  The  data  transmitted  to  the  ground  are  then 
reduced  to  a  plot  of  temperature  versus  altitude,  since  altitude  is  a 
function  of  air  pressure.  The  lapse  rate  is  then  directly  obtainable,  and 
the  existence  and  strength  of  any  inversion  wi’l  be  known.  Whether  the 
inversion  can  be  penetrated  can  then  be  calculated  and  a  test  firing 
scheduled  accordingly.  The  existence  and  strength  of  the  inversion  layers 
need  to  be  monitored  before,  during  and  after  tests  to  predict  hazards.  A 
computer  is  the  only  means  of  accomplishing  thi-,  The  conTidence  in  the 
data  IS  dependent  directly  on  past  experience  and  the  immediate  availability 
of  applicable  data. 
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SECTION  n 


APPLICATION  OF  METEOROLOGICAL  CRITERIA 


The  value  of  meteorology  in  the  efficient  design  and  operation  of  a  toxic 
rocket  test  stand  will  be  found  upon  consideration  of  the  influencing  factors 
pertinent  to  rocket  propulsion  testing.  These  factors  involve  dispersion, 
attainment  of  stack  heights,  the  penetration  of  inversions,  and  the  effects 
of  total  integrated  dosage  upon  test  system  design.  To  attain  the  degree 
of  reliability  required  for  an  operational  system,  a  knowledge  of  specific 
toxicity  levels,  improvement  of  predictabilities,  meteorological  data 
systems  and  sampling  systems  is  required. 

The  following  section  of  this  report  presents  discussions  and  analyses  on 
the  various  toxicity  considerations  and  the  effect  of  a  finite  emission  time 
versus  a  continuous  emission  as  it  would  pertain  to  limitations  in  test 
durations.  Also  covered  are  the  variability  and  reliability  of  toxic  con¬ 
centrations  as  predicted  by  equations  and  as  actually  determined  by 
instrumentation.  The  degree  of  prediction  error  is  calculated  in  order  to 
establish  a  basis  for  improving  the  ability  to  predict  toxicity  levels. 

The  instrumentation  requirements  for  the  determination  of  meteorological 
and  toxicity  parameters  are  presented.  Weather  stations  and  the  data 
reduction  system  required  for  a  complete  dispersion  analysis  are  intro¬ 
duced.  The  toxic  sampling  network  for  determination  of  hazardous  con¬ 
ditions  is  also  related  to  data  reduction. 

It  has  been  found  that  specific  toxic  rockets  may  be  safely  tested  in  the 
proper  environment.  A  complete  analysis  of  the  various  test  stand 
configurations  and  operational  descriptions  will  be  found  in  the  report, 
"Design  of  a  Toxic  Rocket  Test  System",  SSD-TDR-62- 1 37,  (Reference 
30)  which  was  published  simultaneously. 
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TOXICITY  CONSIDERATIONS 


The  necessity  for  any  study  of  the  dispersion  of  rocket  exhaust  products  is 
based  t.-ri  the  observation  that  certain  substances  are  injurious  to  human 
'ife.  TJ.eir  effects  are  observed  in  bolJi  chronic  and  acute  diseases  as  well 
as  systemic  poisons  or  annoyances  That  the  use  of  these  toxic  chemicals 
in  rockel  technology  is  becoming  more  common  is  apparent  The  justifi¬ 
cation  for  this  increased  use  of  toxic  chemicals  is  found  in  their  higher 
energy  content,  improved  storabiiity.  or  other  properties  These  improve- 
nients  in  propellants  have  the  inherent  advantage  to  improve  vehicle  per- 
formrince  either  in  the  increase  in  rani;e  or  payload. 

Propellants  are  genera  llv  divided  into  two  tvpes:  liquid  and  solid,  although 
a  tln'rd  type  i.s  also  pos.'^ible  (hybrid  or  iriprope  llanl)  in  which  both  liquid 
and  solid  phases  occur  The  presence  of  two  varieties  of  propellants  and 
three  possible  forms,  solid,  liquid,  and  gas  for  rocket  exhaust  products 
increases  the  comple.Mties  of  the  toxicity  problem 

The  effects  of  atoxic  substance  upon  anima’.  or  human  life  depend,  in 
general,  on  its  physical  state.  A  liquid  may  be  absorbed  through  the  skin, 
sw'allowed,  entrained  in  air  and  breathed,  or  the  vapors  may  be  breathed. 
Typical  eifects  of  a  toxic  fuel  spill  such,  as  UDMll  or  hydrazine,  are 
chocking,  difficult  breathing,  nausea,  and  upon  e.xtreine  exposure,  con¬ 
vulsive  seizures  and  death. 

Of  course,  these  symptoms  would  not  occur  unless  the  concentration  of 
vapors  are  of  dangerous  proportions,  and  any  concentrations  of  that 
magnitude  would  exceed  allowable  dosages 

There  are  generally  three  limits  of  concentration  specified  for  a  chemical; 
an  environmental  one  for  the  general  public,  an  eight-hour  or  working  day 
concentration,  and  an  evacuation  limit  for  everyone  who  might  be  exposed 
to  the  chemical.  Depending  on  the  chemical,  the  eight-hour  limit  might  vary 
from  5,  000  ppm.  for  a  mildly  toxic  substance  like  carbon  diokide  to  0  01 
ppm  for  an  extremely  toxic  substance  like  pentaborane 

The  methofl  of  poisoning  can  come  from  a  number  of  ways  but,  in  general, 
is  limited  to  inhalation  of  vapors  or  cutaneous  action  for  corrosive  propel¬ 
lants.  The  chances  of  personnel  dosages  by  oral,  intramuscular,  intra¬ 
venous,  or  subcutaneous  means  are  low  with  regard  to  rocket  propellants 
at  test  stands  or  launch  sites.  Before  any  usage  of  toxic  chemical  at  a 
facility,  the  noxious  properties  must  be  tlcfined  and  suitable  safeguards 
provided  for  personnel  These  safeguards  include  safety  clothing  and 
breathing  apparatu.s,  toxic  detection  ecpiiptnont,  and  decontamination 
provisions . 


The  liquid  propellants  can  be  introduced  into  the  body  by  many  different 
routes,  but  solid  toxic  substances  have  only  three  ways  to  enter  the  body; 
through  cuts,  ingestion  and  by  inhalation. 

A  brief  present -ition  t)f  the  respiratory  system  and  its  action  on  solid 
phase  ’-articulate  matter  is  included  to  aid  in  the  understanding  of  deposi¬ 
tion  of  toxic  matter  in  The  lungs.  With  th.e  point  of  initial  intake  being 
ei’iijr  the  nose  or  the  mouth,  dir  passes  in  succession  through  the  trachea, 
bro.ichi  bronchioles  alveolar  ducts  and  finally  into  the  alveolar  sacs, 

The  uvcr-oll  structure  consists  of  a  series  of  branching  passages  decreas¬ 
ing  in  size  and  increasing  in  number  as  the  distance  from  the  point  of  intake 
.ncreases.  Th.e  lidirs  and  small  bones  of  llie  nose  provide  an  efficient 
fihering  svstem  for  large  particles,  while  the  hair-like  cilia,  lining  the 
respirTtorv  ducts  above  the  bronchioles  trap  and  transport  to  the  mouth 
anv  insoluble  particles  tripped  in  this  region.  The  very  fine  particles 
penetrate  to  the  deeper  recesses  of  the  lungs  and  are  deposited  there  to 
varying  degrees  depending  on  their  size.  Particles  deposited  in  the  alve¬ 
oli  andalveo'ar  sacs  arc  ingested  by  migrcitorv  cells  and  transported,  with 
d  feic  passing  into  the  bloodstream,  through  tlie  lung  tissue  to  collect  in 
the  lympliatic  sy.stem  Eventuctlly  the  particles  are  carried  to  the  bron¬ 
chioles  fnen  to  the  ciliated  regions  and  1  rom  here  eliminated  Analysis 
of  particles  size  distributions  in  the  inhaled  and  expired  air  of  human  sub¬ 
jects  and  examination  of  lung  deposits  of  deceased  workers  has  been  used 
to  study  the  mechanism  of  deposition  The  result  of  this  analysis  is  repro¬ 
duced  and  shown  in  Figure  14.  afti-r  Green  and  Lane  (Reference  25).  If  the 
characteristics  of  the  toxic  rocket  motor  are  similar  to  an  aluminum  motor, 
one  would  expect  the  solid  exhaust  products  to  be  in  the  range  of  maximum 
retention  ^ipp roximar cly  0.8  to  1  6  microns  in  diameter. 

The  possibility  cf  poisoning  from  toxic  solid  propellants  is  presently  less 
Than  from  liquid  propellants  because  they  are  not  yet  in  general  usage. 

This  condition  promises  to  be  only  a  temporary  one  since  extensive 
throught  and  research  is  being  expended  for  various  toxic  solid  compounds. 
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PARTICLE  DIMAMETER, 
Figure  14.  Particle  Deposition  Rate 
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LIMITED  TEST  DURATION 


During  rocket  motor  test,  the  emission  is  normally  short  time  firing 
rather  than  instantaneous  or  continuous.  A  study  of  meteorology  will 
reveal  that  short  time  firings  or  finite  emissions  will  never  exceed  the 
downwind  concentration  found  with  a  continuous  point  source  of  the 
same  emission  rate.  When  the  test  duration  is  long,  the  downwind 
concentration  from  a  finite  emission  approaches  the  continuous  point 
source  concentration  as  a  limit.  When  the  plume  reaches  the  point  of 
interest  downwind  before  the  emission  is  discontinued,  the  concentration 
will  be  very  nearly  that  for  a  continuous  source.  When  the  plume  reaches 
about  a  quarter  of  the  way  toward  the  point  of  interest  before  the  emission 
is  stopped,  the  peak  concentration  will  be  about  three  quarters  of  that 
from  a  continuous  point  source. 

The  downwind  concentration  from  a  finite  source  will  never  exceed  the 
concentration  from  an  instantaneous  point  source  if  the  quantities  of  toxic 
products  emitted  are  identical.  The  total  amount  of  substance  emitted  is 
compared  in  this  instance,  while  the  rate  of  emission  was  used  in  com¬ 
paring  finite  with  continuous  emissions.  For  short  duration  and  finite 
release,  the  instantaneous  point  source  plume  is  approximated  and  so  are 
the  downwind  concentrations.  In  fact,  with  increased  distance,  the  clouds 
are  practically  identical  in  shape  and  concentration  distribution.  For  a 
release  time  as  long  as  500  seconds  in  a  turbulent  atmosphere  with  a  ten 
mile  per  hour  wind,  the  maximum  concentration  for  the  finite  release'  is 
about  90%  of  the  concentration  from  an  instantaneous  release  ten  miles 
downwind.  At  greater  distances  the  concentrations  become  essentially 
identical. 

For  short  duration  testing  of  rocket  motors  up  to  several  minutes  running 
time,  the  instantaneous  point  source  emission  equations  are  appropriate 
for  estimating  downwind  concentrations.  The  total  integrated  dosage  for 
instantaneous  point  source  release  can  be  applied  to  finite  release  to  de¬ 
termine  dosages  downwind. 

The  relationship  between  center  of  plume  concentrations  for  continuous 
point  source  emission  and  instantaneous  emission  can  be  expressed  as 
follows; 

X  i  ^  (2-  n)/2 

^c  _  ttE  C  X  _ 

X.  u  0 

1 

Wherej  “  concentration  downwind  of  continuous  point  source 
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-  concentration  downwind  of  instantaneous  point  source 
C  -  diffusion  coefficient 

n  -  Sutton's  stability  parameter 

X  -  distance  downwind  of  point  source 

u  -  wind  velocity,  meters/sec 

9  -  time  in  seconds  required  to  emit  from  the  continuous 

point  source  the  amount  of  material  issued  instantaneously 
from  the  instantaneous  point  source 

The  duration  of  integration  for  determining  concentrations  along  the 
axis  of  the  plume  is  given  as  follows:  (References  4  and  27) 

Duration  of  integration; 

n  5.  48  C 


Where  d' -  distance  downwind  or  is  equal  to  x  and  other  symbols  have 
the  usual  meaning 

By  dividing  the  duration  time  for  a  finite  emission  by  the  duration  of 
integration,  a  fraction  is  obtained  from  which  the  ratio  of  concentration 
for  a  finite  release  compared  with  a  continuous  emission  can  be  de¬ 
termined.  A  chart  in  Reference  27  can  be  used  to  determine  the  relation¬ 
ship,  The  curve  is  a  straight  line  on  log-log  plot  to  the  point  of  60  % 
finite  release  concentration.  The  straight  line  portion  of  the  curve  can 
be  expressed  as  follows; 

log  E.  0.  9813  log  +  0.  4362 
V  01 

Where  =  concentration  downwind  of  a  finite  point  source 

release  for  time  0, 

From  these  expressions,  a  table  can  be  developed  to  show  the  relation¬ 
ship  between  concentrations  downwind  of  continuous,  finite,  and  instan¬ 
taneous  point  source  releases.  This  table  is  based  on  unstable  atmospheric 
conditions,  but  the  relationships  hold  true  for. other  atmospheric  conditions 
with  increased  downwind  distances.  The  first  column  of  the  chart :gives 
distance.  The  second  column  Xi )  presents  the  ratio  of  downwind 

concentrations  for  continuous  as  compared  with  instantaneous  releases, 
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The  period  of  time  required  for  the  continuous  point  source  to  evolve 
the  same  amount  of  effluent  establishes  the  relationship  between  the 
source  strengths  and  are  shown  for  periods  of  5,  30,  and  500  seconds. 

The  third  column  of  the  chart  shov;s  the  relationship  between  Xf  /  X^,  for 

various  finite  release  times  where  the  rates  of  release  are  identical. 

The  ratios  are  shown  to  increase  with  increased  release  time  and; to 
decrease  with  downwind  distance.  The  last  column  shows  the  product  of 
column  (2)  times  column  (3),  which  shows  Xf/X^  .  For  example,  where 
a  five-second  finite  firing  is  compared  with  an  instantaneous  source 
releasing  the  same  amount  of  effluent  in  a  very  short  period  of  time, 
the  downwind  concentrations  are  the  same  within  3%  at  one  mile  distance. 

RELATIVE  DOWNWIND  CONCENTRATIONS  FROM 
CONTINUOUS,  FINITE,  AND 
INSTANTANEOUS  POINT  SOURCES 


X  finite 

X  continuous 

X  finite 

X  instantaneous 

Distance 

y  instantaneous 

X  continuous 

Finite  Firing 

Time  for  Equal 
Total  Emissions 

Miles 

5  sec*  30  sec*  500  sec* 

5sec$30sec^  500  sec| 

5  sec 

30  sec 

500  sec 

1 

24.2 

4.03  0.242 

0.04  0.23 

1.  0 

0.968 

0.927 

0.242 

5 

104.0 

17.33  1.04 

0.0095  0.056 

0.77 

0.988 

0.972 

0.802 

10 

194.0 

3.23  1.94 

0.0051  0.0305 

0.46 

0.990 

0.985 

0.893 

Where  ,  continuous  emission  rate,  is  such  that  an  amount  equal  to 
is  liberated  every  5,  30,  and  500  seconds,  respectively. 

Duration  of  finite  emission  at  rate 

In  summary,  instantaneous  point  source  concentration  equations  may  be 
utilized  for  finite  releases  for  rocket  test  systems. 


76 


CONCENTRATION  VARIABILITY  AND  RELIABILITY 


Dispersion  coefficients  for  diffusion  formulae  have  been  determined  for 
a  number  of  locations.  The  process  of  diffusion  is  influenced  by  many 
conditions  of  the  atmosphere  and  topographic  features.  However,  as 
shown  by  Ocean  Breeze  experiments  run  in  the  Atlantic  and  Pacific 
Missile  Ranges,  prediction  of  ground  level  diffusion  can  be  made  within 
a  factor  of  two  with  a  99%  confidence  level.  Ground  level  diffusion  varies 
more  with  atmospheric  conditions  and  surface  influence  than  does  elevated 
diffusion,  but  data  are  more  accessible  near  the  ground  than  at  points 
above.  With  sufficient  upper  air  measurements,  diffusion  in  the  upper  air 
layers  will  be  as  predictable  as  that  at  the  surface.  A  factor  of  between 
2  to  1/2  the  predicted  concentration  should  be  possible  given  the  height  of 
the  diffusion  process  and  atmospheric  condition  along  the  predicted  path. 

The  height  at  which  the  supposed  diffusion  process  takes- place  is  not 
easily  predicted,  For  ground  level  concentrations,  the  actual  vertical 
profile  of  the  smoke  plume  is  less  important  than  the  ultimate  height  to 
which  the  plume  rises  and  from  where  the  greater  part  of  diffusion  can 
be  assumed  to  take  place.  Various  effective  stack  height  formulae  have 
been  derived  to  predict  the  plume  path  and/or  maximum  height  of  plume 
rise.  The  effective  stack  height  includes  not  only  the  physical  stack  height 
when  used  as  the  source  of  emission,  but  also  exit  velocity  effect  and 
cloud  buoyancy, 

Atmospheric  influence  on  effective  stack  height  is  very  complex,  and 
most  stack  height  formulae  are  derived  for  limited  conditions  and  are 
overly  simplified.  Some  formulae  do  not  consider  effluent  exit  velocity, 
others  consider  plume  rise  within  a  short  distance  of  the  stack,  and  still 
others  neglect  the  effect  of  atmospheric  lapse  rate  on  plume  path.  The 
Thomas  effective  stack  height  formula  has  been  derived  empirically  to 
predict  actual  smoke  plumes  issuing  from  steam  generating  plants.  It 
is  intended  to  yield  maximum  rise  of  effluent  rather  than  the  short  distance 
downwind  used  in  the  Oak  Ridge  formula,  which  is  admittedly  conservative. 
The  energy  release  of  power  generating  plants  is  more  nearly  equivalent 
to  third  stage  ICBM  motors  and  is  probably  most  reliable  for  this  size 
motor.  Smoke  plumes  from  smaller  engines  tend  to  lose  buoyancy  by 
rapid  diffusion  and  larger  plumes  are  rendered  less  buoyant  by  stable 
atmospheric  lapse  rates.  Detonation  of  explosives  equivalent  to  second 
stage  Minutemen  have  produced  smoke  clouds  with  center-line  elevations 
in  excess  of  2000  meters.  Explosion  of  400-pound  charges  have  produced 
smoke  clouds  with  top-of-plume  elevations  of  1500  meters.  The  Thomas 
stack  height  formula  appears  to  be  fairly  accurate  and  possibly  conserva¬ 
tive  for  small  upper  stage  motors,  but  possibly  optimistic  by  a  factor  of 
2  for  very  small  propellant  charges  and  for  second  stage  ICBM  solid 
motors. 
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Assuming  dispersion  and  stack  height  errors  are  cumulative,  the 
resulting  error  in  predicted  concentration  can  be  off  as  much  as  eight 
fold.  Until  data  have  been  collected  from  dispersion  and  stack  height 
measurements,  a  safety  factor  of  10  is  recommended  for  initial  rocket 
testing. 

The  maximum  downwind  concentration  from  an  elevated,  instantaneous 
point  source  varies  inversely  as  the  cube  of  the  stack  height.  The 
distance  to  this  point  is  influenced  by  the  ratio  of  stack  height  over 
dispersion  coefficient  raised  to  a  power  of  1.  1  to  1.  3,  depending  on 
atmospheric  stability.  However,  the  total  integrated  dosage,  which 
may  be  considered  of  more  significance  than  maximum  concentration, 
varies  inversely  as  the  square  of  the  effective  stack  height.  This  means 
an  estimated  stack  height  twice  actual  will  result  in  a  total  integrated 
dosage  of  4  times  predicted.  The  maximum  total  integrated  dosage  is 
not  affected  by  dispersion  coefficients  assuming  constant  atmospheric 
conditions  and  the  absence  of  weather  fronts,  rainfall  and  other  major 
disturbances  between  the  source  and  point  of  maximum  contamination. 

The  difficulty  of  monitoring  for  toxic  products  at  various  locations  and 
at  great  distances  downwind  suggests  a  safety  factor  of  up  to  10-fold  for 
larger  engines.  Small  engines,  which  are  more  easily  monitored  and 
where  area  control  can  be  maintained  might  soon  be  fired  using  snialler 
safety  factors  than  10-fold  as  the  confidence  level  is  improved.  Even 
the  dispersion  and  stability  parameters  for  major  types  of  atmospheric 
conditions  may  vary  as  much  as  50^  from  one  type  to  the  next.  Without  an 
accurate  statement  of  atmospheric  conditions,  prediction  of  concen¬ 
tration  of  toxic  products  cannot  be  precise. 

Actual  prediction  of  concentration  of  toxic  products  is  not  necessary. 

The  intent  is  to  establish  safe  limits  and  avoid  overexposure  of  any 
inhabited  areas. 

In  general,  toxic  limits  as  established  by  the  Manufacturing  Chemists' 
Association  and  the  Atomic  Energy  Commission  are  apparently  conserva¬ 
tive.  Another  safety  factor  of  10  appears  overly  conservative.  However, 
human  lives  are  at  stake  and  the  results  of  overexposure  are  conclusive. 
Once  confidence  levels  are  established  within  a  narrower  range  for  three 
dimensional  dispersion,  a  safety  factor  of  10  may  not  be  required  below 
allowable  concentrations  recommended. 
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JMPROVEMEN  I  OP'  PREDlCTABlLiT  Y 


The  amount  of  error  inherent  in  any  one  case  of  dispersion  can  be  reduced 
to  a  minimum  when  more  complete  data,  recorded  from  previous  tests, 
have  been  added  to  the  prediction  equations.  These  increased  data  will 
include  topographical  features  of  the  test  site,  characteristics  of  the 
tested  pr opeilant  systems,  and  atmospheric  conditions.  The  topographi¬ 
cal  features  comprise  the  mechanical  eddies  generated  by  terrain  and 
structures  in  the  immediate  vicinity  of  the  stack,  and  by  the  stack 
Itself.  Also  included  in  this  category  are  the  terrain  features  downwind 
of  the  stack,  since  these  features  influence  the  surface  boundary  layer 
at  the  distance  of  iriaximum  concentration  'when  dispersed  from  a  stack). 

.Yecessary,  therefore,  to  a  comprehensive  test  program  is  a  thorough 
survey  of  any  proposed  test  sites  to  determine  the  optimum  location.  A 
trade-off  wilt  be  im.  oived  in  which  access  to  the  site,  load  bearing  ability 
of  soil,  and  topographv  related  to  meteorology  are  all  considered. 

The  characteristics  of  the  propellant  system  consist  of  the  toxic 
ingredient's'  and  their  quantities,  the  thrust  level,  and  orientation  of  the 
exhaust.  These  are  easily  determined  parameters  which  are  of  great 
importance  in  predicting  s’ack  heights. 

Observations  of  maximum  plume  rise  can  be  equated  to  rocket  parameters. 
Standard  formujae  for  plume  rise  include  the  parameters  of  gas  buoyancy,  ' 
upward  velocity,  and  atmospheric  lapse  rates. 

Atmospheric  conditions  affecting  dispersion  are:  wind  velocity  and  gusti¬ 
ness,  wind  variation  with  altitude,  lapse  rate  and  variation  in  lapse  rate 
•with  altitude,  vertical  air  density  structure,  atmospheric  turbulence, 
and  wind  shear. 

.if  any  of  the  abo\  e  sets  of  parameters  is  not  properly  considered  in  the 
prediction  of  dispersion  of  toxic  effluents  the  resulting  calculations  will 
not  agree  with  the  actual  observed  concentration  of  contaminants. 

Assuming  ail  variables  of  importance  are  known  and  determined,  the 
accuracy  to  which  they  are  known  is  of  next  degree  of  importance  to  the 
determination  of  the  predictability  of  the  final  result. 

The  probable  error  of  a  calculated  function  can  be  derived.  The  simplest 
example  in  the  present  study  is  the  calculation  of  the  maximum  toxic 
concentration,  X  max  3,  continuous  point  source.  It  is  written  as 


e  TT  u  h^  C^. 
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where  C^/  =  1  for  isotropic  conditions.  To  determine  the  probable 

error  of  Xmax.  the  probable  errors  of  Q  (emission  rate  ),  u  (wind 
velocity),  and  h  (effective  stack  height)  need  to  be  known.  Error  theory 
can  be  applied  to  functions  of  several  variables,  the  best  value  of  the 
derived  quantity  being  obtained  by  substituting  the  mean  values  of  the 
several  variables  in  the  function,  Defining  the  probable  error,  ,  of 

the  best  value,  p,  of  the  function: 

p  =  f  (  S;L  ,  Sg  ,  S3  ,  .  ,  .  ,  Sm  ) 

and  relating  it  to  the  probable  errors,  ,  Qg  ,  Q3  ,  .  .  .  .  Q^,  ,  of  the  mean 
values  m^,  rog  ,  m3  ,  .  .  .  .  m^  ,  of  the  several  independently  measured 
quantities,  ,  Sg  ,  S3  ,  .  ,  S^j  by  the  following  equation: 


^  al.)  q.2 

j  V  j; 

If  this  expression  is  expanded  for  the  present  case,  it  becomes: 


(Ref.  29) 


^  /  (l 


9X 

dQl 


Qq^  + 

^  Idu 


Qu"  * 


An  example  problem  is  presented: 

^max  is  desired, 

and  Q  =  1000,  ±10.  gm/sec 

u  =10  ±  0.  5  meter s /  sec 


Then:  X 


max 


X 


max 


50  ±5.  meters 

2  (1000) 

e  7T  (10)  (50)^ 


gm/m'" 


=  0.  009368  gm/m^ 

=  9.  368  X  10  '  gm/m" 


or 


'max 


=  9.  368  milligram/meter" 
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The  probable  erroJ  then  computed; 


9_A 

2 

9  368  X  1 0”'’ 

9  Q 

e  TT  u  h^ 

/  .  J\J  U  .A.  Jk  V/ 

9  X 

2  Q 

-  9.  368  X  lO""^ 

9  u  ‘ 

e  7T  h^  ( -  u  ^  1 

9X 

2  Q 

‘o'h 

e  TT  u  (-2  h^) 

-  9  368  X  10"^' 

Then: 

/  ~  ^  i  ~  2 

/  q  368  lO'"''  (10)  +  (-9  368  lO"’^/  (0.  5)  +  (-9.  368  •  lO""  )  (5) 

V 

r  — 

/ 

■  /  87  76  10'"  +  87  76  10'®  10.  25.)  +  87.  76  •  1 0"'°  (25) 

V 

6.69  X  10’'*gm/m^ 

Therefore.  Xj-^ax  ^  ®  ^^9  mjlhgram/meter" 

The  number  of  signif.icant  figures  present  in  the  above  values  of 

too  great  A  value  more  m  keeping  with  the  given  values  of  Q,  u,  and  h 

Xmax  9  4  ±0  7  milhgram/meter^ 

and  the  vcilu'=-  is  correct  to  about  7^  From  this  example  it  is  clear  that 
all  of  the  variables  m  an  equation  are  important  when  considering  the 
probable  error  of  the  dependent  variable 


DATA  SYSTEMS 


Mic rometeorological  Stations 

Data  stations  are  a  necessary  part  of  the  toxic  test  facility.  Since  such  data 
stations  exist  and  are  well  within  the  state-of-the-art,  a  thorough  description 
of  these  stations  will  not  be  presented.  For  specific  test  programs  in  which 
the  particular  propellant  characteristics  are  known,  the  rate  of  heat  release 
vs.  the  buoyancy  rise  is  given  (Figure  13’.  The  heights  which  are  attainable 
in  advanced  development  programs  are  greater  than  those  covered  in  pre¬ 
vious  diffusion  studies.  Towers  and  other  equipment  are  required  to  deter¬ 
mine  the  prediction  capability  for  plume  rises  abo\’e  normal  inversion 
layers  and  in  excess  of  3000  feet  altitude. 

Vandenberg  Air  Force  Base  has  a  satisfactory  meteorological  station  net¬ 
work  for  diffusion  studies  (Reference  28)  and  most  of  their  hardware  can  be 
duplicated  for  such  stations  as  required  by  th, .  Rocket  Site. 

It  is  recommended  that  three  towers  150  to  200  feet  liigh  are  provided 
for  the  purposes  of  testing  and  monitoring  the  required  meteorological 
conditions  at  the  Rocket  Site.  These  towers  would  allow  the  obtaining  of 
data  pertaining  to  lapse  rates  and  wind  velocity  gradients,  which  were  not 
necessary  at  PMR  because  non-buoyant  releases  were  the  only  types 
covered.  The  use  of  towers  would  reduce  the  dependency  on  radiosonde  or 
wiresonde  equipment  w'hich  would  otherwise  be  required.  The  data  obtained 
from  towers  are  continuously  ai'ailable  and  will  reduce  pretest  complications 
by  permitting  an  uninterrupted  meteorological  forecast  to  be  made.  A  com¬ 
puter  is  the  optimum  means  for  achiexing  accurate  and  continuous  fore¬ 
casts  due  to  its  data  retention  and  short  time  readout  capabilities. 

These  towers  will  monitor  at  specified  heights  such  parameters  as  wind 
velocity,  wind  direction,  temperatures,  temperature  differences  between 
various  heights,  and  heat  radiation.  The  data  measured  at  these  towers 
should  be  transmitted  to  a  data  reduction  center  for  rapid  analysis.  A 
computer  is  recommended  for  data  correlation  at  this  point  so  that  infor¬ 
mation  can  be  transmitted  as  rapidly  as  possible  to  cognizant  personnel 
conducting  the  tests. 

A  helium  filled  radiosonde  or  wiresonde  capable  of  500  feet  altitude  with 
future  expansion  capabilities  to  2500  feet  is  also  desirable  in  order  to 
determine  lapse  rate.  The  information  transmitred  from  the  radiosonde 
should  also  be  transmitted  to  the  central  data  reduction  center  for  rapid 
data  reduction.  Radiosondes  and  wiresondes  have  inherent  economic  and 
accuracy  drawbacks  which  tend  to  reduce  their  use  to  a  necessary  minimum. 
High  low-level  winds  tend  to  reduce  their  usefulness  since  launching  is 
more  difficult  and  a  vertical  data  profile  cannot  be  obtained. 
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Cameras  are  required  in  addition  to  the  other  meteorological  instruments. 
The  cameras  will  be  used  to  record  height,  size,  and  velocities  of  plumes. 
The  recordings  will  be  made  of  nontoxic  plumes  for  meteorological  studies 
as  well  as  for  actual  test  conditions  when  the  plume  will  contain  toxic 
material.  The  use  of  smoke  plume  photography  is  a  relatively  recent  inno¬ 
vation,  but  is  particularly  appropriate  for  air  pollution  studies,  since 
direct  measurement  of  the  dispersion  caused  by  wind  turbulence  may  be 
obtained,  The  method  has  a  number  of  advantages  but  fails  when  the  visi¬ 
bility  is  low.  Analogous  techniques  to  the  photography  of  plumes  is  the 
use  of  infra-red  cinetheodolites ,  radar  tracking  by  use  of  radar  reflective 
additives,  or  the  addition  of  small  balloons  to  the  effluent  and  optically  or 
visually  following  them.  Cameras  have  been  used  by  NASA  at  the  Langley 
Research  Center,  Langley,  Virginia.  Also,  an  angle-of-attack  sensor 
attached  to  the  nose  of  a  Scout  vehicle  was  used  at  Marshall  Space  Flight 
Center. 

In  an  implementation  program,  it  is  recommended  that  the  three  towers 
and  the  low  altitude  measuring  devices  be  installed  first.  As  the  interest 
in  larger  toxic  rockets  increases,  the  high  altitude  and  camera  instrumen¬ 
tation  should  be  added. 


Sampling  System 


A  rocket  test  stand  configuration  employing  toxic  propellants  must  be 
designed  to  safely  handle  both  the  propellants  and  the  exhaust  products. 
These  chemicals  will  be  emitted  wholly  or  in  part  to  the  atmosphere  where 
they  will  be  dispersed  by  the  macroscopic  turbulent  available  energy.  Due 
to  the  possibility  of  a  large  variation  in  thrust  levels  between  tests  and  the 
change  in  quantity  or  type  of  toxic  e.xhaust  product,  no  single  device  could 
monitor  all  possible  contaminants.  Some  possible  toxic  materials  present 
during  rocket  testing  are:  (before  combustion)  hydrazine  and  its  deriva¬ 
tives,  nitrogen  tetro.xide ,  fluorine,  methyl  alcohol,  nitric  acid,  ammonia, 
pentaborane,  diborane,  bor ohydrides ,  ozone,  benzene,  and  others. 

Upon  combustion,  many  of  the  above  chemicals  become  relatively  harmless, 
but  many  exhaust  products  are  dangerous.  Poisonous  compounds  of  boron, 
fluorine,  chlorine,  as  well  as  carbon  monoxide,  are  sometimes  present 
and  must  be  considered. 

Because  of  the  extremely  low  exposure  limits  for  the  more  toxic  com¬ 
pounds,  the  concentration  must  be  monitored  to  determine  hazards  at  the 
test  stand  and  in  the  surrounding  area. 
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Monitoring  techniques  have  been  developed  for  use  with  toxic  liquid  pro¬ 
pellants,  such  as  hydrazine,  UDMH,  and  nitrogen  tetroxide.  These 
chemicals,  for  use  with  the  TITAN  II  and  other  weapons  systems,  must 
be  carefully  handled  and  every  precaution  must  be  exercised  to  prevent 
spills.  When  spillage  does  occur,  the  presence  and  spread  of  gaseous 
vapors  must  be  detected  quickly.  Examples  of  excellent  existing  systems 
for  this  purpose  are  found  in  those  developed  for  the  Atlantic  and  Pacific 
Missile  Ranges  (Reference  5),  Here,  a  controlled  emission  of  zinc  sulfide 
aerosol  was  allowed  to  disperse,  and  upwards  of  500  air  samplers  were 
used  to  detect  its  concentration  at  locations  downwind  of  the  source.  The 
samplers  were  positioned  in  two  and  in  one  case  three  arcs  in  excess  of 
90°  spread  in  the  prevailing  wind  directions  and  in  a  pattern  near  the 
aerosol  generator.  The  data  collected  by  this  system  were  used  to  derive 
highly  reliable  dispersion  coefficients  for  use  in  calculating  diffusion  rates 
for  spills,  In  addition,  a  series  of  equations  for  the  prediction  of  aerosol 
concentration  with  distance  was  developed  from  these  experiments. 

A  similar  program  was  performed  by  the  Stanford  Research  Institute 
(Reference  24),  at  Edwards  Air  Force  Base,  to  determine  dispersion 
coefficients.  The  coefficients  may  be  used  in  the  Sutton  diffusion  equations 
to  calculate  downwind  concentrations.  The  Stanford  Research  Institute 
study  was  performed  to  determine  the  hazards  of  a  hydrogen  and  fluorine 
test  stand  at  the  Edwards  Rocket  Site.  Ground  level  dispersion  was  the 
only  factor  investigated  in  this  study. 

Because  of  a  high  level  of  toxicity  possible  with  the  more  exotic  fuels,  a 
program  to  evaluate  dispersion  from  elevated  sources  is  required  and  would 
use  a  more  elaborate  data  gathering  system  than  used  heretofore  at  Edwards, 
Sampling  stations  are  required  both  at  the  Rocket  Site  and  outside  of  the 
boundaries  of  Edwards  Air  Force  Base.  Appreciable  concentrations  of 
toxics  could  exist  at  even  hundreds  of  miles  downwind  of  a  high  thrust  firing 
and  the  hazard  at  ground  level  must  be  known.  It  has  been  shown  that 
sampling  equipment  up  to  at  least  150  miles  downwind  of  the  Rocket  Site 
would  be  required  for  direct  release  into  the  atmosphere.  A  total  of  approxi¬ 
mately  75  units  would  be  needed.  Three  central  data  gathering  units  in  45° 
arcs  would  be  the  minimum  for  safety. 

Some  meteorological  systems  use  more  than  75  sampling  units.  At 
Vandenberg  Air  Force  Base  (PMR)  hundreds  of  zinc  sulfide  samplers  were 
used  in  the  development  of  the  "Ocean  Breeze"  formulae.  Also  at 
Vandenberg  samplers  are  required  for  momforing  propellant  spills  of 
hydrazine,  UDMH,  and  nitrogen  tetroxide. 
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A  small  toxic  test  stand  would  require  onl\  a  few  samplers  if  the  propel¬ 
lant  toxicity  is  of  a  low  order  whereas  the  larger  systems  require  greater 
sampling  complexity  for  safe  operation.  The  rate  of  toxic  emission  and 
the  toxic  properties  ot  the  propellant  wil'  determine  the  sampling  system. 

The  actual  equipment  needed  to  accou.j.  ;!(  •,<.  c t .'on  of  toxic  propel- 

ants  or  exhaust  products  may  vary  widely  in  form,  operation,  and  cost. 
For  the  volatile  chemicals  a  number  of  methods  are  possible;  gas 
chromatographs  or  sensitive  chemical  detector'^  being  but  two  ways  to 
detect  toxic  gases  Toxic  -olid  materials  require  different  methods  of 
analysis,  PTrst  panicle-  ol  finite  sue  have  to  be  collected  and  then  must 
be  analyzed  This  could  be  done  by  air  samplers  which  would  gather  the 
toxic  solids  close  to  the  level  at  which  they  would  be  inhaled.  The  sample 
may  be  collected  on  a  t.'.lter  paper  for  analyzing.  The  detection  of  radio¬ 
active  waste  material  w...ll  probably  become  necessary  in  relation  to 
propulsion  system-.  Therefore,  the  means  of  detecting  a  specific  toxic 
effluent  w>ll  be  chosen  as  rocket  test  schedules  require  since  new  propel¬ 
lant  combination'-  are  continually  being  di.scovered. 


Data  Handling  Sy.stem 

A  data  handling  system  is  an  important  au.xiliary  to  a  toxic  test  facility 
and  will  be  used  to  assist  in  the  prediction  of  meteorological  conditions 
and  indicate  whether  a  test  could  be  safely  initiated.  Input  into  the 
system  would  include  data  from  radiosonde  balloons,  local  weather 
stations,  general  weather  reports,  and  facility-based  instruments.  In 
addition,  the  type,  orientation,  and  thrust  level  of  the  toxic  propellant 
rocket  would  be  added  to  the  general  equation  of  the  computer.  The  safe 
firing  time  would  be  obtained  from  this  input.  In  addition,  a  plot  of  pre¬ 
dicted  dow'nwmd  concentrations  as  a  function  of  real  time,  distance,  and 
changing  weather  would  be  possible  to  derive 

Implementation  of  such  a  system  may  follow  a  number  of  paths.  The 
least  expensive  would  be  a  manually  operated  calculator  with  the 
calculations  done  only  as  required.  The  most  complex  system  would  be 
a  computer  providing  continuous  plots  of  test  conditions.  The  chosen 
system  for  the  requirements  presented  will  be  the  optimum  of  these  two 
extremes  It  would  give  the  go  -  no-go  prediction  based  on  short  time 
inputs  The  time  lag  between  the  taking  of  meteorological  and  other  data 
and  the  go  -  no-go  prediction  would  be  a  dominant  factor  in  the  choice  of 
data  handling  system  complexity. 


TYPICAL  applications  TO  ROCKET  TEST  SYSTEMS 


Design  of  a.  test  system  depends  on  the  development  of  design  criteria. 
Motor  sizes  and  characteristics  are  of  highest  importance  in  design 
considerations  since  every  effort  must  be  made  to  obtain  data  of  the  actual 
test  item  operating  conditions.  Where  highly  toxic  propellants  and  pro¬ 
pellant  exhaust  products  are  involved,  measures  must  be  taken  to  prevent 
overexposure  of  personnel  to  toxic  products.  The  test  system  must  be 
designed  to  provide  the  protection  needed.  Anticipated  test  frequency  for 
various  size  motors  must  be  established  to  prevent  excessive  set-up  time 
and  undue  scheduling  of  test  runs.  The  frequency  of  testing  is  important 
in  determining  allowable  release  of  toxic  products  both  from  the  stand¬ 
point  of  total  dosage  at  inhabited  areas  downwind  and  the  probability  of 
suitable  weather  conditions  for  safe  dispersion  of  toxic  products. 


Once  the  test  program  has  been  established,  meteorological  principles 
can  be  employed  to  estimate  the  toxicity  hazards  that  can  result  from  such 
a  test  program.  A  study  of  meteorological  phenomena  will  suggest  ways 
of  using  atmospheric  dispersion  to  the  greatest  advantage.  Emission  of 
effluent  from  one  large  stack  is  recommended  over  many  small  stacks. 

The  ideal  stack  for  maximum  dispersion  of  effluent  is  tall,  large  in  diameter, 
and  with  high  exit  velocity  and  buoyancy  of  effluent.  The  stack  should  be  high 
enough  to  prevent  downwash  due  to  turbulence  caused  by  structural  and  topo¬ 
graphical  features  near  the  stack.  A  stack  at  least  twice  as  high  as 
obstructions  should  be  empio\  ed. 

Current  meteorological  information  will  be  required  for  the  test  site  to 
determine  go-no  go  conditions  where  climatic  influences  can  result  in 
hazardous  conditions  from  a  test  or  malfunction.  Sufficient  meteorological 
measuring  instruments  must  be  provided  to  distinguish  between  safe  and 
unsafe  firing  conditions.  Data  reduction  equipment  with  the  capacity  to 
make  current  predictions  of  safe  conditions  must  be  provided.  Sophistication 
can  range  from  simple  pen  type  recorders  of  wind  velocity,  direction, 
variability  and  temperature  lapse  rate  to  an  elaborate  network  of  instruments 
feeding  information  into  a  computer  for  continuous  prediction  of  toxicity 
limitations.  The  more  complex  the  meteorological  station, the  greater  the 
initial  cost.  Maintenance  and  operating  costs  will  not  necessarily  increase 
with  the  complexity  of  the  meteorological  system.  For  highly  toxic  exhaust 
products  where  great  distances  may  be  involved  for  adequate  dispersion,  the 
weather  network  and  toxicity  monitoring  system  may  become  so  extensive  that 
greater  economy  can  be  realized  through  application  of  safety  measures  and 
devices  at  the  source.  Without  sufficient  monitoring  devices  for  weather  and 
for  toxicity  determinations  to  improve  confidence  levels  for  dispersion  de¬ 
terminations,  reliance  on  atmospheric  dispersion  is  dependent  upon  local 
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dispersion  determinations  and  effective  sta.ck  height  measurements. 
Sufficient  area  meteorological  information  must  be  available  to  determine 
rainfall  and  frontal  conditions  between  the  test  facility  and  downwind  areas. 

An  economic  balance  between  test  system,  meteorological  network,  and 
operating  costs  should  be  established.  Further  considerations  include  the 
value  of  meteorological  data  accumulated  for  use  on  other  programs  and 
estimated  savings  that  may  be  realized  in  future  facility  design,  The 
meteorological  station  will  be  useful  and  possibly  necessary  for  scheduling 
toxic  liquid  propellant  handling  operations  anticipated  in  the  near  future. 

A  highly  sophisticated  micrometeorological  network  is  needed  for  spill 
tests,  future  toxic  liquid  propellant  transfer  facilities,  small  toxic  solid 
motor  testing,  and  local  monitoring  of  large  toxic  solid  motor  test  system. 

One  method  of  simplifying  a  test  system  design  and  capital  investment 
would  be  the  addition  of  a  flare  to  the  rocket  or  stack  exhaust.  This  is 
well  within  the  state-of-the  art. 

The  graph  shown  on  Figure  15  is  a  plot  of  diffusion  from  a  point  source  with 
no  stack  height  for  ground  emission.  The  source  strength  is  doubled  to 
account  for  reflection  of  particulates  or  gas  reflection  from  the  ground.  Alon 
the  axis  of  an  elevated  plume,  the  concentration  will  be  half  these  values. 

The  ratio  of  X  /Q  plotted  with  respect  to  distance  to  show  dilution  ratios 
for  maximum  convenience  in  use.  The  dilution  ratios  required  for  selected 
motors  are  marked  for  concentrations  equal  to  one  microgram  per  cubic 
meter.  The  distance  required  to  these  dilutions  points  out  the  necessity  for 
elevated  release  and/or  buoyant  plumes  to  yield  adequate  dilution  at  ground 
level.  The  slightly  diverging  series  of  lines  represents  various  atmospheric 
stabilities  from  extremely  turbulent  to  strong  inversion  Dilution  ratios  at 
a  downwind  location  vary  a  thousand-fold  between  climatic  extremes.  How¬ 
ever,  from  an  elevated  point  source,  the  maximum  ground  level  concen- 
tation  remains  essentially  the  same,  but  the  distance  to  the  maximum  is 
increased  ten-fold  when  climatic  conditions  change  from  very  unstable  to 
very  stable. 

The  graph  shown  on  Figure  16  is  a  plot  of  the  empirical  formula  derived 
from  groundTevel  emission  under  Project  Ocean  Breeze  (Ref.  5).  The 
preliminary  report  under  this  project  presented  four  formulae,  one  for  each 
major  atmospheric  stability  category.  However,  it  was  found  that  changing 
from  one  formula  for  one  lapse  rate  to  the  formula,  for  the  adjoining  lapse 
rate  category,  discontinuities  resulted  of  greater  magnitude  than  the  devi¬ 
ation  resulting  from  one  all-inclusive  formula  including  a  lapse  term.  A 
high  confidence  level  has  been  established  for  this  new,  pre-release  formula. 
A  99%  probability  of  predicting  within  one  half  or  double  the  measured  con¬ 
centration  has  resulted  from  the  use  of  this  formula  for  ground  level  point 
source  emission.  For  ground  emissions,  such  as  toxic  liquid  spills,  the 
Ocean  Breeze  formula  is  very  promising 
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CONCLUSIONS  AND  RECOMMENDATIONS 


The  study  of  the  effect  that  meteorology  has  on  toxic  rocket  systems 
has  resulted  in  certain  basic  conclusions.  These  conclusions  must  be 
rather  general  until  a  specific  propellant  system  is  chosen,  after  which 
more  exact  calculations  may  be  made  from  definite  parameters  of  the 
propellant  toxicity  and  the,  associated  thrust  levels.  The  conclusions 
are: 

1.  Toxic  rockets  may  be  fired  directly  into  the  atmosphere 
with  the  following  qualifications. 

a.  The  test  stand  location  must  be  a  safe  distance  from 
inhabited  areas. 

b.  Any  limitations  on  thrust  levels,  and  toxic  emission 
rates,  must  be  maintained. 

c.  The  effluent  should  be  directed  in  a  vertical  direction 
for  maximum  natural  stack  height. 

d.  The  weather  conditions  must  be  suitable  or  artificial 
equipment  must  be  used.  This  could  include  flaring 

of  the  exhaust  gases  for  greater  buoyancy  of  the  plume. 

e.  The  weather  conditions  should  be  predicted  to  be 
constant  and  near  ideal  for  the  duration  of  the  test. 

f.  Air  sampling  must  be  conducted  to  insure  that  toxic 
limits  are  not  exceeded. 

Z.  Toxic  rocket  tests  of  large  size  may  be  conducted,  but  tests 
should  not  be  md.de  where  the  maximum  allowable  concentra¬ 
tions  (MAC)  might  be  exceeded  beyond  the  borders  of  a  test  site. 

3.  In  addition  to  the  atmosphere  diffusion  of  toxic  rocket  exhaust,  a 
filtration  or  scrubber  system  could  be  added  to  allow  greater 
toxic  emission  rates  prior  to  filtration.  Trade-offs  would  have 
to  be  made  between  the  two  extremes  of  complete  filtration  to 
complete  emission  of  the  toxic  components.  From  these  trade¬ 
offs  the  facility  could  be  optimized  for: 

a.  Maximum  personnel  safety  and  confidence  in  system  • 
performance. 
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b.  Optimum  cost  where  open  air  firing  would  require 
elaborate  and  costly  meteorological  networks  and 
experimentation,  and  where  complete  filtration  would 
require  elaborate  and  costly  blast  containment  and 
filtration  netccorks. 

c.  Maximum  facility  utilization  based  on  setup  and 
time  percentages. 

4,  A  health  and  safety  group  must  be  lormed  to  establish  and 
enforce  regulations,  establish  sampling  stations  for  back¬ 
ground  readings,  and  obtain  and  monitor  physical  examina¬ 
tions  ot  test  personnel.  This  group  s  activities  would 
increase  in  proportion  to  test  size  and  frequency. 

Based  on  these  conclusions  and  the  stud\  of  applicable  meteorological 
analytical  techniques,,  the  following  recommendations  have  been  developed 

1,  Atmospheric  diffusion  for  highly  toxic  rockets  of  small  size 
can  utilize  point  source  ground  ie\-el  emission  equations. 

Either  Sutton's  form  or  the  "Ocean  Breeze"  form  of  point 
source  ground  level  emission  formula  can  be  used.  Where 
the  downwind  distanc  e  to  the  control  area  limit  is  over  five 
miles,  these  equations  ha'  e  led  to  the  conclusion  that  small 
rockets  may  be  tested  during  turbulent  atmospheric  conditions. 

The  maximum  safe  size  would  ha\’e  to  be  calculated  for  the 
specific  propellant  combination, 

Z.  Large  toxic  rockets  mac  require  a  control  area  greater  than 
five  miles  downwind  to  pre\  ent  exceeding  the  weekly  allowable 
dosage.  It  is  important  to  strive  for  maximum  effective  stack 
heights  and  use  Sutton  s  elevated  instantaneous  point  source 
formula.  "Ocean  Breeze"  formulae  are  for  non-buoyant 
emissions  but  could  be  modified  for  elevated  plumes.  Dis¬ 
persion  coefficients  appropriate  for  the  effective  elevation  of 
the  source  should  be  determined  in  either  case. 

3.  Dispersion  coefficients  ha\  e  been  established  for  various 

atmospheric  condi'ions  and  the  variation  with  altitude  delineated. 
However,  experimentation  in  this  a,”ea  should  be  accomplished 
to  improve  reliability  for  long  distance  predictions  where  very- 
toxic  exhaust  products  are  found.  Reliable  dispersion  co¬ 
efficients  for  elevations  of  several  thousand  feet  will  be  required 
for  large  toxic  engines. 
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The  effective  stack  height  must  be  determined  from  the  relation 
specific  to  the  emission  type.  That  is,  the  recommended 
procedure  for  an  explosion  is  to  use  the  graphical  relation 
between  heat  release  and  cloud  height.  A  horizontal  release 
will  use  the  buoyancy  formula,  and  a  vertical  firing  will  use 
the  buoyancy  formula  plus  the  jet  penetration  equation;  the 
buoyancy  equation  was  developed  in  this  study.  The  effective 
stack  height  can  be  augmented  by  flaring  the  fuel  rich  exhaust 
gases  or  even  by  adding  additional  fuel  to  the  exhaust  prior  to 
flaring  to  provide  a  greater  buoyancy  to  the  gases. 

5.  Inversion  penetration  has  been  discussed  and  is  shown  to  be 
possible  with  large  engines.  Where  the  high  stability 
associated  with  an  iin-ersion  results  in  low  dispersion  and  high 
downwind  concentrations,,  piercing  the  inversion  can  result 
in  the  lowest  downwind  concentrations.  Experimental  work  with 
non  toxic  charges  or  fuel  flaring  should  prove  the  reliability  of 
inversion  penetration.  Release  and  flaring  of  ordinary  fuel 
could  be  used  prior  to  a  toxic  test  to  prove  that  piercing  can 
be  accomplished  for  the  particular  test  climate.  Large  motors 
should  be  capable  of  piercing  inversions,  and  smaller  motors 
may  be  augmented  to  do  so. 

6.  Air  sampling  stations  in  inhabited  areas  must  be  set  up  prior 
to  and  during  test  phases  to  monitor  background  concentrations 
and  extent  of  contamination.  A  few  of  the  more  obvious  areas 
where  these  stations  should  be  located  are  at  the  rocket  test 
stand,  the  support  facilities,  guard  houses,  nearby  public 
highways,  and  dowmw’ind  communities.  Other  air  sampling 
stations  both  in-plan*-  and  out-plant  will  be  evolved  from  the 
type  of  test  facility  system  used  and  at  the  direction  of 
cognizant  health  safety  personnel. 

Based  on  the  analysis  of  meteoroiogicai  phenomena,  a  conceptual  design 
may  be  prepared  for  specific  toxic  test  system.s. 
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ABSTRACTS 


The  following  abstracts  of  reports  and  books  are  presented  in  alphabetical 
order.  These  documents  cover  meteorology,  toxicity,  and  rocket  technology 
applicable  to  this  study, 

Barad,  M,  L,  ,  Project  Prairie  Grass,  A  Field  Program  in  Diffusion, 

GRD  Research  I^te  No,  59,  (July  1958)  Astia  No,  AD  152572, 

A  field  program  was  designed  to  pro\ide  experimental  data  on  the 
diffusion  of  a  tracer  gas  over  a  range  of  800  meters.  Seventy  tests 
were  conducted  in  which  the  gases  were  continuously  emitted  for 
ten-minute  intervals.  The  releases  were  made  over  flat  Nebraska 
prairie  under  a  variety  of  meteorological  conditions  during  July  and 
August  of  1956.  Tabulations  of  the  diffusion  data  and  the  meteorological 
data  collected  during  the  gas  releases  are  also  presented, 

Barad,  M,  L.  and  Haugen.  D,A,  ,  "A  Preliminary  Evaluation  of  Sutton's 
Hypothesis  for  Diffusion  from  a  Continuous  Point  Source,  "  Journal  of 
Meteorology,  X\  I,  No,  1,  (February  1959)  12-20, 

In  Sutton  s  equation  the  obseri'ed  concentration  distribution  is  pre¬ 
dicted  only  if  there  are  two  \  aiues  of  n,  a  ny,  and  a  n^.  Statistical 
tests  indicate  that  ny  and  n^  are  invariant  with  distance  between  100 
and  800  meters  of  the  source,  but  that  the  values  of  Uy  and  n^ 
appropriate  for  these  distances  exceed  the  values  within  100  meters 
of  the  source.  It  is  also  shown  that  neither  n^  nor  n^  can  be  speci¬ 
fied  by  n,^y,  the  value  of  n  found  from  a  power-law  fit  to  the  wind 
profile  in  the  lowest  8  meters. 

These  conclusions  were  obtained  from  data  from  Project  Prairie 
Grass,  The  v'alues  of  h  were  50  centimeters  and  1,  5  meters  over 
level  countryside  covered  with  grass  of  a  height  of  about  5  centimeters. 
These  conclusions  probably  do  not  present  valid  indications  for 
macroscopic  rocket  tests  where  h  is  much  greater  than  1,  5  meters 
and  where  ny  is  approximately  equal  to  n^. 

Beers,  Norman  R.  ,  "Stack  Meteorology  and  Atmospheric  Disposal  of 
Radioactive  AVaste,  "  Nucleonics,  IV,  No.  4,  (April  1949)  28-38. 

Discusses  the  problem  of  air  cooling  of  piles  since  argon  in  air 
becomes  radioactive  vvith  a  half-life  of  110  minutes.  Article  covers 
Sutton's  equations  for  a  continuous  point  source  and  gives  the  oil-fog 
smoke  test  method.  Concludes  that  for  the  high  reactor  stacks  at 
Brookhaven  there  is  no  danger  of  radioactivity  and  indications  are 
that  this  IS  true . 
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Best,  A,  C.  ,  "Maximum  Gas  Concentration  at  Ground  Level  from 
Industrial  Chimneys  "  Journal  of  the  Institute  of  Fuel.  XXX,  (June  1957)  197. 


The  maximum  ground  level  concentration  of  toxic  effluent  and  the 
distance  to  the  maximum  depend  upon  the  stack  and  meteorological 
characteristics,  Co'  ers  sei'eral  formmae  lor  computing  stack 
height,  giving  results  of  calculations  and  aMempts  to  show  that  for 
the  purpose  of  computing  the  maximum  ground- level  concentration 
the  three  are  in  reasonable  agreement. 

Bosanquet,  C.  H.  ,  "The  Rise  of  a  Hot  Waste  Gas  Plume,  "  Journal  of  the 
Institute  of  Fuel  XXX  June'  19571  322. 

it  IS  assumed  when  a  cloud  of  hot  gas  is  rising,  the  total  heat 
content  and  total  upward  momentum  are  unaffected  by  dilution  with 
air.  Also  assumed  that  upward  momentum  increases  at  a  rate 
proportional  to  the  heat  content.  A  set  ol  equa'ions  for  stack  height 
IS  dern-ed. 

Bosanquet,  Carey,  and  Hallon.  "Dust  Dc-position  from  Chimney  Stacks," 
Proc,  Inst,  of  Mech.  Engrs,  CLXIJ  .’1950,  155. 

A  general  description  of  'he  interaction  of  stacks  and  diffusion  is 
given.  Derives  formuJae  lor  piume  rise  and  diffusion  of  dust  with 
rates  of  deposition  of  dust.  N'umerous  nomographs  are  given  to 
assist  in  working  method, 

Lramer,  Record,  and  Vaughan,  'I  he  .Study  of  the  Diffusion  of  Gases  or 
Aerosols  in  the  Lower  Atmosphere,  Mass,  Ins',  of  Technology,  Dept,  of 
Meteorology  'January  1959,  Astia  Xo-  AD210482, 

Diffusion  measurements  were  taken  at  Round  Hill,  Mass.  ,  of  a 
similar  nature  to  Project  Prairie  Grass,  These  experiments  were 
for  the  determination  of  instantaneous  diffusion  of  coefficients  at 
ground  level.  Emission  times  were  of  30  second  and  three-minute 
intervals  and  are  compared  to  concentration  levels  obtained  at  10- 
minute  intervals.  A'-erage  concentrations  were  2.5  and  1,5  times 
larger  respectively  than  for  10-minute  periods. 

DeMarrais,  G.  A.  ,  Workbook  in  Atmospheric  Diffusion  Calculations, 

IJ,  S,  Weathe r  Bur eau,  Idaho  Fans,  Idaho,  iFebruary  1959), 

This  workbook  w'as  prepared  as  a  practical  supplement  to  Meteorology 
and  Atomic  Energy,  Included  are  the  behavior  of  stack  effluents,  the 
determination  of  n  and  the  coefficients  of  diffusion,  the  application 
of  the  basic  equations  and  nomographic  solutions  of  the  equations. 
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Estoque,  M.  A,  ,  Venting  of  Hot  Gases  Through  Temperature  Inversions, 
GRD  Research  Note  No.  3  December  Astia  No.  AD  160756. 

The  penetration  of  temperature  inversions  in  the  lower  atmosphere 
by  plumes  of  hot  air  is  investigated  as  part  of  a  safety  analysis  in 
the  nuclear  ai  rcraft  program.  This  program  is  studied  with  the  aid 
of  existing  theoretical  as  well  as  experimental  work.  A  nomogram 
showing  the  relationship  be'ween  the  maximum  height  attained  by  a 
hot  piume  of  a  given  heat  source  inter  sitv  and  the  temperature 
gradient  of  the  env;roi  mvnt  is  presented. 

Gifford,  F.  A.  Jr,  ,  The  Problem  of  Forecasting  Dispersion  in  the  Lower 
Atmosphere,  Wi  ather  Bureau  R(  s^-a  r  ch  St  at  ion  Oak  Ridge,  Tennessee 
(July  IQbFlT 

As  the  nuclear  po^er  industry  develops  v'eathe  r  forecasters  will 
be  called  on  mon.  and  more  to  make  estimates  of  dispersion  from 
isolated  sources  in  the  lower  atmosphere.  The  various  equations 
and  parameter  values  needed  for  this  purpose  are  given.  The 
generalized  Gaussian  plume  model  is  presented  in  some  detail, 
and  diffusion  parameters  based  or  recent  observational  and 
theoretical  studies  a-e  d’seussed. 

Glasstone  S,  "'ihe  Etfvct-,  ot  Nuciear  Weapons.  Published  by  USAEC, 
Washington,  D,  C.  ,  'April  1?6Z  , 

The  subiects  covered  are  the  ge.verai  principles  and  descriptions 
of  nuclear  explosions  blast  effects  the  effects  on  personnel,  and 
principles  of  protection.  Tabular  data  of  interest  included 
descriptions  of  ar’-ouriced  nuclear  detonations  in  which  the  name  of 
the  shot,  date  and  time  of  ig.nition.  lomatior.,  height  and  type  of 
burst,  yield  and  the,  altitude  of  the  cloud  top,  cloud  base,  and 
tropopause , 

Green,  H.  L.  and  Lane  W.  R.  ,  Particulate  Clouds;  Dusts,  Smokes,  and 
Mists,  D.  Van  Nostrand  London,  T957-. 

The  process  by  which  dust  is  passed  into  the  lungs  or  trapped  by 
the  cilia  ard  other  parts  of  the  respiratory  tract  is  covered.  A 
curve  for  deposition  in,  the  lung  of  particles  of  unit  density  is  presented. 
The  maximum  rate  of  deposition  is  for  one  micron  diameter  particles 
and  over  12  microns  the  rate  is  rii. 

Greenfield,  Harvey,  Theoretica,.  Derivations  ard  Nomographic  Methods  for 
Sutton's  Diffusion  Equations,  Research  Directorate,  Dugway  Proving  Ground, 
Dugwav,  Utah  .September  10,  .956,;  Astia  No.  AD  1  1  3307. 
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A  sur\  ey  is  made  of  diflus'.on  equaiions  including  hypotheses 
underlying  derivations  and  limits  to  applicability.  Nomograms  are 
derived,  and  their  use  is  demonstrated  for  Sutton's  diffusion 
equations.  The  use  of  nomograms  is  shown  to  simplify  and  expedite 
diffusion  calculations.  Six  computed  nomograms  are  included  with 
derivations  and  demonstrations  of  their  use. 

Haugen  D,  A.,  Barad,  M.  L.  and  Anianaitis,  H,  ,  "Values  of  Parameters 
Appearing  in  Suiton  s  Diffusion  .Modeis.  "  Journal  of  Meteorology,  XT’lIl, 

No,  i  ■'June  196  1.  i68-37Z. 

S'atisticai  summaries  arc  presented  ol  computed  values  of  para¬ 
meters  appearing  in  d  modilicaiion  of  Sutton  s  diffusion  equations. 
Experimental  da'a  used  for  ilie  computdlions  .are  those  obtained 
during  Project  Prairie  fjrass  at  OA'cn.  Nebraska  and  those  obtained 
at  RuLuid  Hill  .Massac  liLiset' s.  St  ratil  .cation  ol  the  parameters  ny 
and  n^  according  to  stabiiitv  c  lass  revea.s  systematic  variations  of 
Class  median  'aiues  of  n^  and  n^  vMth  stability  class  and  a  large  range 
of  iiy  and  n^  wnhui  any  go  eii  sttilxlitv  ciass.  No  difference  in  the 
results  was  found  between  the  two  sets  of  experiments  for  these  para¬ 
meters. 

The  parameter  has  a  mean  vabie  of  about  0,  4  for  the  Prairie 
Grass  experiments  and  .aliout  0.  Q  lor  the  Round  Hii)  experiments.  No 
relationship  between  Cv  and  stabiirv  was  found  for  the  Prairie  Grass 
experiments,  but  a  correlation  cuelfuient  of  0,  between  Cy  and 
stability  ratio  was  found  for  the  Round  Hill  experiments. 

The  parameter  C^,  lomputed  oniy  for  the  Prairie  Grass  experiments, 
has  a  mean  value  of  about  O.OV  and  is  independent  of  stability  for  those 
experiments  permitting  a  siat  i  si  i  cal  D  -  st  able  estimate  of  C^. 

Hawkins,  J.  E,  and  Nonhebel  .  (j,  "C  himnevs  and  the  Dispersal  of  Smoke,  " 
Journal  of  the  Rst.tute  of  Fuel,  XX ’v’.',.  November  1955)  178, 


The  authors  d..S(  uss  effects  of  disturbance  of  a  plume  by  passage  of 
wind  past  the  emitting  chimney.  The  smoke  is  drawn  down  to  ground 
lev'el  by  trailing  vortices.  A  plot  of  emission  v  elocity  vs.  wind  speed 
at  which  downwash  wnl  not  occur  is  given.  Contains  discussion  of 
Oak  Ridge  fhomus  and  Bosanquet  formulae. 

Hill  Thomas  and  Abersold,  "Dispersion  of  (jases  from  Tall  Stacks,  " 
industrial  and  Engineering  Chemist  ry  XLI  No.  11  November  1959) 
24^^ 

The  Thecreticdl  ec[uaIions  of  BoS'inquel  and  Pearson  and  Sutt'on  for 
the  dispers  on  of  sm<’>ke  [’•'"m  f'>ctc'rv  chtmnevs  hive  been  solved  in 
’erms  o.’  tlic  convent  ona!  uit's  of  the'  smelting  industrv  The 
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theoref.ical  curves  ard  corfirmi-^ig  data  itilustrate  forcefully 
the  be'^eficial  effects  of  the  use  of  tail  stacks  in  dispersing  air 
contamir  ants  from  factories  ir  creases  in  stack  temperature  increase 
effective  stack  height  ard  improve  dispersion. 

Launch  Siting  Criteria  tor  High- Thrust  Vehicles  Aeronutronic  Division, 
Ford  Motor  Company,  Te  chri  cai  Re  port  No.  U- 108:118,  (March  31,  1961). 

Launch  ha^ards  are  described  <vhich  include  acoustical,  explosion, 
and  toxicity  hazards  Rocket  e-.gi’"'e  ^oise  a’ d  its  propagation  and 
effects  at  various  sound  power  le  f  is  are  cove.'-ed.  Propellant 
explosio’"s  are  described  along  with  th*^  '•esulting  fragmentation  and 
overpressure  a’'d  a  tabic  of  dista*-crs  oased  o-  available  data,  is 
derived,  in  the  toxicit,'  hazard  section,  the  diffusion  of  poisonous 
propellants  ir  the  atmosphere  '.s  ii'i-^cnted  based  on  .Sutton's  equations. 

Magill  Holder,  a.ndAckle>,  A: :  Poljutio’'  Handbook.  McGraw-Hill,  New 
York  ;1956l. 

Presented  are  various  facets  of  pol.utio"  .  The  section  on  calculation 
of  concent I'atior:  of  pollula.-^ts  gi  as  the  ■ariation  in  the  generalized 
eddy  -  diffur  lo’"  coe  f!icie’'ts  :th  he  ig.ht  ard  .vith  lapse  rate,  The 
statemert  is  made  that  .Sn"0'^  s  -al  ;es  for  maximum  concentration 
from  a  sou’^ce  are  o-'d^'r  of  mag',  rude  correct. 

Meteorology  and  Atomic  E'er.;’,  AFCU  3066,  Published  by  USAEC, 
Washington  D,  C.  jul/  i955.. 

This  IS  cor.S'dered  to  be  ore  of  th'"  best  'oaimes  or  meteorology 
as  applied  to  poiiutior  by  air-borre  vastes,  especially  radioactive. 
Presented  are  outlines  oi  diffusio”,  theories.,  bsha’.uor  of  stack 
effluents  and  explosio'*  debris  ciouds  fall-out,  wash-out,  and  rain- 
out  from  air-borre  Ciouds  ard  i-adioacti’  cloud  dosage  calculations. 
Also  co'-ered  are  graphical  soiutio'ts  to  a*.mo5pher;.c  diffusion 
problems  reacto'  hazard  a''.alyse.->  a- d  a  selection  of  recommended 
equatior.s  pa.'-amete -s  a- d  conve rsi or  factors. 

Moses,  H,  and  Strom  G.  H.  "A  Comparisor  ot  Observed  Plume  Rises 
with  Values  Obtained  from  Weil-K’rowr  Formulas,  "  Journal  of  the  APCA, 

XI,  No.  iO,  October  i96i'  455-66, 

Actual  observatiors  were  made,  oi  over  one  hurd.“ed  different  smoke 
runs,  and  the  measured  v.aiues  of  stack  height  were  compared  with 
calculated  ores  from  six  sources:  1  Holland  'Oak  Ridge),  (2)  Bryant 
and  Davidson  (3'  Sutton,  '4  Sco  -e.-  5  Bosa-quet,  Carey,  and 

Haltoi"  and  6  Bosarquet.  Co'  eluded  that  'o  one  formula  is  out¬ 
standing  ir.  all  respects  a-rd  the  ra'.ure  of  the  problem  would  influence 
the  selection  of  the  o’'r  us^d, 
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Priestley,  C.  H.B.,  ''A  Workiiijj  theory  of  the  Bent-Over  Plume  of  Hot 
(las,  "  Quarterly  Journal  of  the  Ruy.il  iVleteorolo<>n  al  Society,  LXXXII, 

1  Q  5  6  .  -  7  6  ~  ' 

An  ascendiny  plume  of  hot  pa.->  is  init'a.:^  diffused  by  the  turbulence 
which  IS  induced  by  its  own  moiiun  ,ind  ,atcr,  but  before  all  upward 
iTiotion  and  buoyanc\  are  lost,  bv  the  natural  turbulence  of  the 
environmenl  A  woikinp  soiulion  is  jiresi'iited  for  the  first  phase  in 
terms  of  spread. np  cuefficieni  (.1  and  for  th"  second  in  terms  of  a 
mixMiL!  rate  K  ,  and  ,1  is  shown  that  the  transition  from  one  phase  to 
the  otl'.er  sIiolocI  be  quite  wel,  marked 

Robinson,  E  ,  M'.  c  rornet  eoroiop-c  ai  and  D.  ffusion  Study  for  High  Energy 
Test  Faciiitc,  S’antord  Research  nstitute  Proeject  \'o,  SE-1964, 

Contract  \o  AF  0-4  o  1  1  -  2  3 1  a  'August  30,  Pia?, 

Detailed  s'udies  uen  intide  of  thie  a;r  puiiuf.on  problem  posed  b^ 
a  planned  high,  cnerg\  ti-s'  f<'ic..itv  at  Edwards  Air  Force  Base. 

Both  theoretical  and  experimenta,  phtises  of  the  problem  were 
included  in  th.s  st  ud\  An  expcrimi  nttil  fieid  program  was  carried 
out  for  the  det  e  rmin.ati  cm  of  dillus  uii  loelfuients  and  stability 
parameters  at  E.duards  AF'B' 

Scorer.  R  .  S.  ,  "I^uimcs  from  I<..1  (diimnecs."  Weather,  X,  'April,  1955 

lOb-lOQ 

A  general  discussion  of  th.e  probU'm  !  he  (  onclusions  are  that  the 
gaseous  products  of  c  ombustion  are  omy  harmful  when  they  are  at 
the  ground  The  objecti'  e  should  be  'here  fore,  to  get  them  as 
high  into  the  air  as  poss'.b.e  bv  means  of  few  tall  wide  chimneys 
and  ha’  e  a  certain  amtauil  of  buovam  \  The  solid  particles  have  a 
harmful  effect  ;it  a.,  'C’-eis,  and  it  is  the  reirioval  of  these  that  should 
Claim  the  engineer's  attent’on 

Spurr,  G,  ,  "I  he  Pen  trutiun  of  A'mospheric  .ii' ersions  by  Hot  Plumes,  ” 
Journal  of  Meleoroiogy,  X-  February  1959  30-  37 

The  article  covers  a  derivation  of  the  magnitude  of  heat  sources  to 
penetrate  inv  ersion.-^,  Equat.on  for  height  is  a  function  of  lapse 
rate  and  hea'  reiease  rate  Preducts  that  the  worst  inversion  on 
record  can  be  penetrated  bv  a  power  plant  of  100  megawatts. 

Stern,  A.  C.  .  Air  Po.lution,  Academic  Press,  Hew  York  flObZ',, 

-  olume  1 . 

This  book  covers  air  pouution  and  its  dispersion,  the  effects 
of  air  poiiution,  measuring  and  monitoring  air  pollution.  A 
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discussion  of  stack  height  formulae  and  dispersion  formulae  is 
very  thoroughly  covered.  The  effects  of  air  pollution  on  humans 
and  animals  are  also  presented. 

Stokinger,  H.  E.  ,  "Air  Pollution  and  the  Particle  Size  -  Toxicity  Problem, 
Nucleonics  Magazine,  (December  1949) 

Included  is  a  general  discus sionof  natural  air  pollution  such  as 
dust  storms,  and  then  man-made  pollution,  which  is  divided  into 
two  types:  (1)  community  contamination,  such  as  smog,  and  (2) 
those  occurring  within  the  plant  itself.  The  toxic  substances 
covered  are:  sulfur  dioxide,  fluorine,  and  dusts  of  beryllium 
compoimds.  The  methods  of  investigation  for  finding  beryllium 
air  contamination  are  given.  The  effect  of  hydrogen  fluoride  on 
beryllium  poisoning  is  presented.  The  toxic  effects  of  beryllium 
sulfate  mist  were  found  to  have  been  doubled  by  the  addition  of 
concentrations  of  hydrogen  fluoride.  The  same  concentration  of 
hydrogen  fluoride  was  absolutely  without  effect  in  the  control  animals 
when  inhaled  alone. 

Sutton,  O.  G.  ,  "Dispersion  of  Hot  Gases  in  the  Atmosphere,  "  Journal  of 
Meteorology,  VII,  No.  5  (October  1950)  307-312. 

The  problem  of  the  disposition  of  a  stream  of  hot  gas  from  a 
point  is  considered.  It  is  shown  that  a  plausible  assumption 
concerning  the  mecl'ianism  of  entrainment  of  air  by  the  jet  leads 
to  simple  expressions  for  the  mean  temperature  and  mean  velocity 
of  a  jet  of  hot  air  rising  in  a  calm  atmosphere  of  uniform  potential 
temperature.  The  theoretical  expressions  are  shown  to  agree  with 
the  laboratory  measurements  of  Schmidt,  and  it  is  found  that  the 
coefficient  of  diffusion  for  these  conditions  is  of  the  same  order  of 
magnitude  as  that  which  is  derived  from  the  large-scale  spreading  of 
cold  smoke  in  the  atmosphere.  Finally,  an  approximate  solution  is 
given  for  the  shape  of  the  plume  from  a  hot  source  in  a  horizontal 
wind;  and  it  is  demonstrated  that  the  reduction  of  maximum  concen¬ 
tration  at  ground  level,  caused  by  adding  heat  to  the  effluent  from  a 
stack,  is  directly  proportional  to  the  strength  of  the  heat  source  and 
inversely  proportional  to  the  height  of  the  chimney  and  the  cube  of 
the  horizontal  wind  speed. 

Thring,  M.  W.  ,  Air  Pollution  ,  Butterworth ' s  Scientific  Publications, 
London  (1957) 

The  chapter  on  the  "Flow  of  Chimney  Gases"  is  by  C.  H.  B.  Bosanquet. 
He  derives  the  stack  height  formula  covered  in  the  text  above. 
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Toxic  Hazards  of  Rocket  Propellants,  Aeronulronic  Division,  Ford  Motor 
Company,  Technical  Report  Ho.  U  108:99  ^November  30,  I960)  Astia 
No.  AD253236. 

This  report  discusses  the  general  subject  of  toxicity  and  the 
method  for  obtaining  and  the  philosophy  for  establishing  the 
maximum,  allowable  concentration  fMACV,  Eighteen  liquid 
rocket  propellants  in  particular  are  discussed  along  with  their 
toxic  hazards  and  operating  procedures  now  used  to  minimize 
the  dangers  of  each  to  personnel.  Atmospheric  diffusion  was 
covered  with  Sutton  s  equation  and  modifications  thereof  used. 
Conditions  of  fallout  and  rainout  with  ground  deposition  and 
integrated  concentrations  are  included.  The  sources  are 
classified  as  continuous  point  line  and  \-oiume  sources.  The 
special  hazards  of  total  instantaneous  washout  as  caused  by  a 
rain  storm,  and  fumigation  as  caused  by  a  temperature  inversion 
are  included  with  typical  dispersion  problems. 

Turbulent  Diffusion  in  the  Atmosphere,  Technical  Note  No.  24,  World 
Meteorological  Organization,  WMO.  No.  77,  TP-31,  Geneva,  Switzerland, 
tl958) 

The  early  theories  of  Taylor,  Schmidt  and  Richardson  as  well  as 
Sutton  are  included  and  Sutton's  equations  are  given.  The  diffusion 
from  a  stack  source  with  fo’^mulae  for  stack  height  and  maximum 
ground  level  concentration  is  given. 

The  Use  of  Effective  Stack  Height  for  Dispersal  of  I  oxic  Rocket  Exhaust 
Gases,  Marquardt  Corporation,  Report  FE -223-2,  Contract  No. 

AF  04(6ll)-4304  (April  20,  1961. 

The  feasibility  and  safety  of  atmospheric  dispersal  of  rocket 
exhaust  gases  at  Edwards  AFB  are  studied.  Conclusions  are  that 
under  suitable  conditions  large  quantities  of  toxic  gases  can  be 
safely  discharged  directly  into  the  atmosphere,  "ncluded  in  the 
study  are  rockets  in  the  50,  000  to  800,  000  lb.  thrust  range,  using 
fluorine  and  chlorine  trifluoride  oxidizers  and  hydrogen  and 
hydrazine  fuels.  Stack  height  is  shown  to  be  the  most  important 
single  parameter  influencing  the  safe  dispersal  of  the  gases  in  the 
atmosphere. 
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APPENDICES 


The  meteorological  aspects  of  toxic  rocket  testing  have  been  investigated, 
and  the  areas  requiring  further  mathematical  analysis  have  been  covered. 
The  synthesizing  of  an  equation  relating  stack  height,  stability  constant  (n), 
wind  velocity,  and  diffusion  coefficient  has  been  done  to  facilitate  the  de¬ 
velopment  of  an  equation  for  buoyant  plume  rise.  They  will  be  found  in 
calculation  sets  A  and  B.  The  momentum  rise  of  a  stack  emission  was 
developed  empirically  by  wind  tunnel  testing  and  related  to  the  system  of 
units  used  in  this  report  in  calculation  set  C.  Dimensional  analysis  was 
used  to  develop  a  new  equation  for  the  momentum  and  buoyancy  rise  of  a 
plume,  and  is  presented  in  set  D. 

The  degree  to  which  errors  in  parameter  measurement  cause  the  calculated 
functions  of  error  theory  to  functions  of  several  variables  in  set  E.  A 
comparison  of  ground  level  concentration  versus  distance  for  gaseous  and 
particulate  emissions  is  done  in  set  F. 
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EQUATION  FOR.  OBTAINING  VALUES  OF  C 


Since  the  basic  derived  equation  for  plume  rise  due  to  buoyancy  includes 
values  of  "C",  an  equation  is  required  which  is  based  on  some  value  of  h 
to  solve  for  "C".  This  is  due  to  the  inciuded  mass  of  a  cloud  increasing 
because  of  diffusion  effects.  Using  the  extended  values  of  "C"  plotted  in 
the  "Meteorology  and  Atomic  Energy"  nomogram  as  functions  of  h,  n  and 
u,  an  equation  was  derived  which  included  these  values  plus  a  constant. 

An  alternate  solution  would  have  been  to  solve  the  plume  rise  equation 
for  various  heights  and  selecting  values  of  "C"  by  successive  approxi¬ 
mations  until  close  fit  values  are  obtained  for  each  calculation.  The 
equation  of  the  nomogram  was  solved  assuming  the  curves  were  plotted  on 
an  equation  of  the  type: 

y  =  CX*  X|  X3' 

where  y  =  "C"  Diffusion  coefficient 

C  -  Constant 


X, 


h  Plume  height 


X. 


n  =  Lapse 


X3  -  u  -  Mean  wind  velocity 


\'alues  for  "C"  were  obtained  from  the  nomogram  for  various  h,  n  and  u 
values  and  the  basic  equation  solved  algebraically  to  yield  "C"=  (3.  67  x 
1  0  ^  j  (h  ■  )  (n"-’’-!!  }  (u”'*'^^  )  With  arbitrarily  chosen  values  of  n 

and  u,  the  plume  rise  can  now  be  expressed  as  a  function  of  the  heat 
released.  A  check  of  the  equation  versus  nomogram  derived  values  of  "C" 
showed  very  good  agreement.  However,  due  caution  must  be  exercised  in 
the  use  of  this  equation;  it  is  doubtful  if  values  obtained  at  h<10  meters  or 
h  >  1000  meters  will  be  of  sufficient  accuracy  to  justify  the  use  of  the  equa¬ 
tion  under  these  conditions. 
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on  -LAt^o  -  Loi<^C  -  -o.9?Xh 

E^OANiAJi^  I'A^uPnoAJ  O’ ^  Fty^  h>(f  C 

0^  y^6-  t~  no.69Ai  h 

Su  8SrtT  u  TiAif-  D-^-f  />7V?  /T-T  /^A'O  .h 

-0.9^1-  -  d?.^  oC  ^0.6996  •’O.X9H  ~  o.doi.  h 
-0.1  (DZ-  0.09']  y 

h=  ~3.II 

Oo&'ST/TC/TAAId-  \//]L-U£':>  j9/i/TO  r/ £/P  /'"'J  ! c?  P/./iyO  d 

P'i  "l.x^d  ~  io&  C  -0.79IJ  '  0,9>d^-{-:5./l) 

yofc-  -'hx'jdfoo^i- i.Hir 

Loy  C  -  ~^>  ^935 
C  -  ^‘(’1 

^OE>PTA^A/  X-i  X}f 

/  —  7  ^  /■  '*^'  \  /•  \  J  \ 

XXi  )(^3  ) 

^tcpjJT 

Rc/iVIA^P  fppM  f]/^>rf<Nd^N^^ir7 


^  h 

n.  U 

0.09^  )00 

a2j'  jO 

A/rLA  T'lA.^Z-  T/i/rd  T}fL 

6zd 

(p.099  -(3.1,1  xto^ 

Yl^OO  )(00,%D'  ){10  ) 

DmEL.  MHH,  JOmOHj  MEHDENHALL 

JOB  NO. 

SHEET  No.  A  ~  Y 

OESIONCD  BY  U/*^  ^ 

DATE 

AMROVIO 

■f"  d 

Gzl  2  v^j'  -^.  /  r  y^  tJ 

^=  -'O.U^ 

^  =  (jyx  yxja  ^)[  ^  /in  ){^  / 

(Z^C  d  ^  T/c'a^ 

I,  J)~  rd  r^tfrt-yzf  ^ 

U’S'  c^  y/j '/-  f-d^’/^yY/.d  d9)  -d {'■'?. //-K-^,  </^ yy-oddY/adY) 
■z -di.'^ 3^- d>.  ^~S7?  i- /.‘V9C  ”  <::>< //r 
~  -l.'doi  . 

c  -  d9.dPV9Z.  /yV^y<S^/i/r^r 

X,  h~-  hL7ty<-^  ,  Id  -  '5 ne-nypc/Sizc^c^  /}  =dX,  2- i' 

UdC  ~  -^'‘i3:>-i’l^d/^lX'2.i-ifi~)-h  (-3Jiy-'C>,(>c?'i.')t  {~2P./i,^)Yp.Hl^r) 

'  o.  3^^ -hi.rix  -  dd’On 

'  -1.013 

(^  -  0>  o9y  ~  (O.dyyj 

1  h-  ^  li  -  /r'A/rr<i'/^j'/;/-'^^^-o^ 

ido di'-x.^js' ■/-y-d>,/^yyx.<^oi)'fC-^d/Y~o.  d9^)4C-o.)69Xi'/-)c 
~-x.H2r-o.9H  i-  x.  nr  -  o.iT) 

^-o.  rrx 

C  -  o.iZ I  /XAfdo/yy^di  ~  o,  13 

h-  Yldn'/if  ^  1/  '  X/1t-^iV^s^fiYdy•'2>  ^  y}~  O,  J3 

Ud-  d-  -2.  V  xr  i  C^d-f  V7Yuo2.)  i-Y-3'/0^0,in]  +  (-Or  ddoYo-  3di ) 

-  -X,9jr-  o.m  t  /.rao  -OrOr/ 

~  -h  36  9 

_ ^  -0.09X1  A/a/^nd/dA  '  O.  09 ^ _ 


111 


DANIEL,  MANN.  JOHNSON.^  MENDENHALL 


=■  e-i 


DCBIONED  BY  Hj  Ia/A 

DATE 

APPROVED  1 

1  .  . 

CONT/A'UOUS  FD/A/r  SQIJRCI^  PLUh^  RISE  DUB 
TO  3U(Jf/^NCY  BECAUSE  OF  CONTENT  /l3dUT 

T)A''/BJTA/T 

T/i^  W/D77/ OF  A  SFAO/CT  PIUME  AFOPl  COA/77A/UOUS 

PO/^T  SOUf^CE  /S  OET/A/ED  6y  77/^  POLLdWiA^Cr 

-  A.  /Ap\k  .  _  , 


avo  2  ■^) 


PH£AP  P  ^  %  OP  CENTEPUA/E  C0A/c£P7R/]r/0N 
yPT  BPOB  OF  CLOUP^  /07^  /5  FFEOUEnTlY 
U^ED  AND  APPPDX/MATB'^  OESEPUABLE  UM/TJ. 

/o  =2{2.-ipCy'^  ^3.02i-C7:.^ 

77/£  PRODUCT  OF  A  PAPAMETEP  COA/TPVNEP  VY/TP/P 
TP/E  CLOUT  UMIT  OF  /^-/<^%  W/LL  RE  LES^ 
yy/AA/  /D6p>.  OF  TEE  TOPIL  ,  PEAT  IN  7?/E  CLOUD 

A  TA/F  FARAMETEA  BE/N6  CONSIDERED, 

CONTOERINE  A  FORIZONTAL  CfL/NDRlCAL  SEE/'ICNT 

OF  TTF  CLOUD ^  THE  TOTAL  FEAT  coa/TA/Z/FP 

\NITF/T  TRLF  SECFJE/VT  of  l/NJT  LEFOT/  CAlF  pE 
FOUNT  .  FT  ///TFCAA'u'FE  TFE  FOLLOwzNO  EXPRES?/LN 


(2  =  y  ^rrXyC  Cz) 

'  ft) 

TRC  u  A- 

(see  text  FCR  DEFNIT/ON  of  tepms() 
W/HBflE  Ji  ^  c  /y"^ 

©  , "  fcvH 


X 


rrCu.<^'^ 


BUT  C  E 


/of 


PPC^N  0) 


X  PPPE] 

HBice  Q„=  L-UAxT- T  Yo-  / 

J  ITd  X’ 


MoPibE 

IZ. 
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DAHIEL.  MAHN,  JOHHSOHj  MEHDEHHALL 


JOB  No. 


DEBIONED  BY 


APPROVED 


iv  U/ 


SHEET  No.,??'  "2. 


DATE 


Yo  /S  D£P/AJ£i?  r?/£  po/^r  V^H£P£.  77/^ 

CQA/CBA/mAT70A/  0/2CaPP£p  /6%  o/- 

Th/^.  C0A/r^R^UA/P  COA/CBVrM770A/  CAR 


/)NP 


OR  THE  /d7o  CONCEMTRAVbN  EP6E  CQNTAJA/S  96^'  Or 
TPE  TOTAL  PARAPETS  A  (  /-/BAJT). 


/ATMOS PRESS  ORE  PS//f/7y  PCA.  STAA/PAPD 


CONDfTTDNS 

VAR/ES 

APPROX IMATELY  ACCOPD/A/G 

THE  POLL  OVU/ATE 

B/PPEST/ONS 

Pi 

=  ro- 

■•  s“  /  1  S'/  2  b  4 

2, 1ESXM  hj  . 

C7) 

R 

“  S }) 

1.  -Z  ES//0  hJ 

&)  ' 

THE  PLUME  /SSUES  HORl-ZONTP  LLY  DUE  TO  TfE 
W/ND  W/TP  TEE/yiO  QUOYANCY  A/EGLECT/P& 

EjdT  VELOCITY^  /an  E/PPESS/6aJ  pipAiP 

P/SE  PAS  /A  p'cyA/CTVOA/  COE  POp/zoz/TAL  OETAACB 
/S  PEQU/RE  POP  ^  RUOYAaJT  PLUME,  ' 

h^f  C'dpxs) 

PRO/]  EXPERIMENTS  PERPCRMEo  BY  3RYAXT  PAW 
COW£>REY  (rep  a) 


As  BEST  APTRLXIMAVAJ  OF  LAJA, 


TEE  OAmo&E  formula  (reP  i)  ms  DE/irUEU  for  //eat 
TELEA\SES  OP  ASAAUT  7E0O  BTV/sEC  AND  700  hETPPS  LOVJhVAiND, 
TA/E  OAKRJO&E  FORMULA  FOR  EFFECT/E  STACK 
//E/GPT  .  Buoy  A  A/ or  term  only  /s  — 


DANIEL.  MANN.  JOHNSON, i  MENDENHALL 

r~ 

JOB  No. 

SHEET  Ho.  S  Y 

DCeiONED  BY  (/J  LA/ A 

{ 

DATE  I 

APPROVED 

Ak-  (0  A-U  ~  BUdY^NC/  RISi^  FEBT 

0  Q-  H£Ar  RciBASE  CAl/s£r^^ 

U  -W/NO  VELOCITY  METERjysCC. 

coiweRT/N(:>  uAiir:  Qt,  Cp^tv/sec); 

h  -  0^03  h  =  pLUnE  PsluT  Due  77; 

U  6(;6YANCY  at  Bi^Of^.ETDU 


h  =  pLUnE  PsluT  Due  77; 

6(;6YANCY  AT  3^0  MEYT/iS 
(9i^^heat  helease.  erj/sEc 
0  -  MEAM  ^IND  vELOC.cn' 


VY HE/IE  iSi^-  E>EOO  BW/SEC.  AND  U  =  pj  MElERS/iEc 


h  -  OJ6^  y  35^0  -  S.6C  METEPA 

f.  Y7 

BcEAVY^  PBPIVEO  F/iCYi  M  EASUREMENIE  r71KEN  AT 
BOO  METBFE  pAc^O  FAdh  A  FAUX  FFIFASF  OF 
HFAT^  77//J-  BQUAFAnj  /D  r?//F  A  ay  7'  AVAILABLE 
EDO  ACE  OF  JNF6RAA7?6C/  /-EVe  PFuF'YM//N/V6 
r//E  FROPOA-fyOWAL/py  OOA/syoi/’r  FDR  FLUFF 

tfa  JFctofy,  'H/b  oak  and  cb  formula  is 

FFK}ub:p  for  feutfal.  stability  vvm/  a 

FACJd/l  F6F.  STAFF  A  FA  LAPSE,  TNBRBFiFB  //=  ZZTALTCF 

A  =JlYtA  ^  ZIlAFlL  =.  2.4'S 

Q^KO^J  3S0o(s6E)‘O^ 


TY/BREFOPB  /)  z  "2,  SI  A/0  ^ 

FNO  ^  ~  JlI^^ - -  ^/'^  U 


.  h2jL _ 

(2.53  Aid  o.Ao 

l,~7X  /  S'T /~\L7S  /j 


li7S  I  Y 

^  w  //Q  (E) 

FROM]  -p-IELSE  EQUAF/ON^  ty/E  CLOUD  'yCLu/uE 
PUE  TO  PIEFUS/OM  AFR  AVPVAS'^TC  ESFP/J/'Yf/DA 
F-/B  CLOUO  WF/Ocyr  A1/VYP  F/SAFCEO  AYR  CA/Y  BE 
PBTBAF/YA/BP, 


l  l  A 


DAHIEL,  JOHHSOtl^  MEt^OEHHALL 


JOB  NO. 

SHEET  No.  B  ~~ 

DEBIQNED  BY  l>J  {yj  A- 

DATE 

APPROVED 


VOLUME  OF  PLL/ME.  DjE  TO  PJFFU5ION  CONSlDELINO 
ONE  lECC^NO  £AVSS/6^/  FROM  CON>T/NUQU^  PO/NT  SOUPSF 

^  fTr^u  BUT  r ^  /.sn  C Qo) 

TF£  P/SPBKS/QN  COFFF/C/Fajt  j  W/LL  PFCTEASF 
VU/TF  ALJ77V,0/F  ACCOFOING  TD  /F/£  FCLLdVjJA'O 
BAPKEOyO/V' 

BY  SUESVTUF6/Y 

^  r  -i/rO,/97Y 

-Tro  [Af/?x3A7  x/o  A0>1S  Ao  J/f 

V7HBPE  X)  -  0>'^E  Fbfl  FTAN[)APJ>  uap^B 

,  ^  2  2  -<5. 27y  //7r 

\/,  ^  F  U  X/.S/7  x3,^7  X  k  U  ^ 

(O^XF)  //o  ^ 

y  7  <-  ,  -  7.  r  /A-  >’^ 

but  7^ '  =  ^  XI0_ 

1,  5-4 


2(V/2  7  2  7/ 

EE  IX  X  /.s/7  X3.i7  X  /O  h 


7/206 


{ii)/  2 


^^/2  2 


X  2.54  ACaP/P'-^ 


/?  2.  <//2  /  5-  20C. 


VOLUME  OF  cloud  from  /  SECOND  EPlSS/ON  /A/CLUP///& 
0/fFUS/ON  and  BXPFNS/Oa/ 


V,.=  V4  +f 


V.  ^  fP)'^  3XT  F=  F  aEA) 

V  [pj 


5",20'i» 


ANO  ^  r  — A - r  // •  Z.  2Si'7/&‘' b) 

/.  /  -  Tj\~  j) 

muL  -  -  j.7P'/A2.2r-r //i  /p-i.ios xhj  hj  Bi, 


-IN  07 


I  / 

£/i^  r  ^  ^0^J(/-2.F so  /'If  'AAJ  " 

pit 
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DAHIEL.  MANH.  JOHNSOH^  MENDEHHALL 


JOB  NO. 

SHEET  NO.  A.'  ;3' 

- 1 

DEBIONED  BY  1>J  v'J  A 

j  DATE 

/I  v  '-//Z  )  j"  ^  ^ 

/S^r  K  =  ^  1^^/- 2/ //v'jr/o  ^  ^ 

/IW-?  Vc  'h  f  ^  Jli 

^  j  a ‘\. 


TA/d'/c^Fry/:^  VC^LVHIE  6I~  CLbUO  I3'y  0/FFl/S/OAA  AND  EX  PANS!  0/^  IS- 

/z  l,L//-i  j  3,ZCiC, 

14  =  2,  //V^/0  _ 4 


/z  i,L/n ,  3,xciC 
14  =  2,  //V^/0  4 

2-V/2  /. 

2,//^/X^^  /„  7 


^  ? '  7S~^/J(^  2.2b" yx  jj(l  -2,2  Ef  )!JCi^}j 


\NEjy^T  OF  CLOUO  yPOPl  /  SECtND  EMiy/ON  F/i6M  C.O/'rr/rpJO^y  Pi/'Ti 
SOU/ZC&  /Aj  AMa//£AfT  A/P  /S  D/AFliPyO  ////O  77/e  CT^auO 


0^  -  7- 'i^'T  XIO 


\0  r  2. 


A  ,  2.'//Z 


I’^OL 


j  ,  IZ  2.V'2  /  jj 

etX,  T  3,200  x-L,!!/  KIES  u  - _ 4 


iNHQZe  F  ^  PE/D  try  Or  /'/>'< 

'  /-T  HBKEDT  h 

i  VOLUNIE  Or  /LO/P 

pue  772  /?//^/7:7//7yy 

c  DENS  try  or  Atti 
AT  h~o  (MrrBTj 
h  r  Z/F/^///'  ADe!/E. 

e/EojNo  C^T/r]  ix/b) 

r\erBRE 

U  ■:=  \r>lNO  vELOCiT/ 

hBveps  /second 


Te'/BPEFOiZE  \//a/6//7'‘  OF  CJ-oOO  /j 

/)  ^  /i  ^ ^  ^  y,zsX  ,  z.zor:  . 

u/,  :=  /  ><  (>fz7  m  u  J ( i^zFsrxto  Aj  qA  Qe) 


VVOAK  00E{B  ABAINST  GRAU/T}'  f3Y  CLOUD  AB  /r  F/fFS 

AJsp/77A/e  s dunce  /j  ujep  for  Acc:oA///ijs^//A/(c  ^ 

WCaRK  of  077/EFU//JF.  FOf/  ^  73  UOYAf  T  PLUMB.  W(0XK=-W. 

//  /  ■ 

+  Ia/^^  c/Ij 

|-bUT  ^  p,  ^  l-z>  ZfFX/O  hj 

J^or^  ^TANDAf^O  ATKl^W^HefAE  TO  3^000  X7]  / 


1  1  L 


DANIEL.  MANN,  JOHNSON,  &  MENDENHALL  I 

JOB  NO. 

SHEET  NO.  /S"/^ 

delbioned  by  U;  /\ 

I  -  - - — 

DATE 

APPROVED 

7!^ ERE FO HE  71 /E  B)(PHf£^S/6N  RO/l  WO/ZK  -QECOIAES 
U^'^'  /(“V'cPd+  Q^) 


W 


/f.  /I  2i  4/^  I  3,^04  7^  /  -S'  >,^,2  51^  ,. 

rfpJ/VX/O  U  _ h _ {[,[}- 2, HB7XI0  hj  Jh 

/L  "  J  .j, 

2  1//2/  ’’S\f/C/  ■  -cy  N- y<  7^*^ /?' 2^353 

-  U  '  ^■^ncL^j.^.-LSS^ia^y-l,  2SSX10  hJ  h 


+ 

o 


I'l' e,  ^  '"'^‘‘dL^Jl,  Cio) 


=  Zil  Lli/MS77/Z^  /^r  S5A  i-euei^  1/V  /;w  L5-PIE1S^S> 

IZUT  ay  c^aa/)  aoi/fvoA/v/  oee/a^/76a/  90^^  c>/=^ 

TOTTiL  RE  AT  TElEA!>ED  W/LL  BE  Cc>A/rA/AED  //K  CLOUP 

Afyp  Z17  LE-MErars  =  j  aru,  -///^^aroae  ia/'Z/3iS/., 

T//£BETO/za  emANCf  EaUAV(^^  m  /SaCCNO  RELB/^SE  of  plumb  IS^. 

Z/zQl-JP _  -2.//V/i'/0  X  l.(^T  /-l,ZBSx/o  hJ  //  R// 

Qu ‘ 


li  h 


-  z 


•  // y  xio ’bns‘/j{^- 2, zs7//a¥/jj^i  -I'l^ sx/o  ¥  d h 


X  ^1-2,255-  Jh 

h  A 


0  < 
(? 


-ZiLUi.-ny X 10 'Yni T'''‘‘'‘ c/ljcM^ 

Oil) 
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OAHIEL.  MAHH,  JOHHSOf^^  MEHOEHtiALL 


JOB  NO. 

- 55 - 

SHECT  NO./)  /  1 

DESIONED  BY  UJ  UJ A 

DATE 

APPROVED 

1 

- 1* 

/3/A/0M:'ALS  70  TBR-l^'lS  d.  r /5c;d /?//’/«- 


-  i.2b%  ^2.2irr^Y/u  Jh 


+  25’C  h.  ZirL)(2  .IS'S'X/'^  ^ 


)  y 


-  V/2 


^2ir4Y 2  2  Si"//a  J  /j 


S-\  >  /  (>,Z6P~J 


3/264  .  . 

-(-V,7^Vy/2  2i5/^6 


_5\,  V/264  j 


Vs^j-5jsyi7, Z‘SSXil^y  k 


2  /  STi  20t> 


^-A-7r‘;Y-^.7r'/y; 
6  , 


-  ^-V/7IvY~6”-/ ip^Y-L7shY ? ,2s.yA7^] f/'i, 26r7.>v67, 


v-2i‘yj 
zj.r.ojil  Up 


2,6V  X  (cpiy xio 


4-  %JShilLzi^yrLS7xEoj\ 


■  2.2-OU  /  -7)  :?'26C^ 

h  -  V/ 2^4/2 -255 //o  //i 


'SUl  V.204 


/  A/  N-  -s']?  / 

-  4^/  2  3, 2r(Y 2  issx/PJ  h 

iNresfwriNS  th^  mjr  77/7 /r 


^  ^  3  Qj,  __  _  ^5,  5-^/  XI 0 

,  ,  2,  VV2 


V/  204 


I  '(•  ,  „,  /  -S-\  /  5/264, 

_i2 _ _  .  v/2-4  (z^i'i.zss Xio  jL 

V/  2o  6  5"/  70  <6 

/  6, 26  G  /  \/  \/  26  (i 


•f  7,  2  S'6 /  j ,'2 6%/ 2 /  2. 55  ?,'/0  yy  ^  ^  V'/,72“6  (^£3)Yi£^(2.25j/{’/0'']  // 
2X6/20  4  \  ^  6X7/'2.oC',  y 

-  7  ,  /  /  V  X  /  6 ''  3-  75 vY  2/  75 5X /  0  ^ J  /  _  (~  V  7  5  'Y 1  /  7  55  X  /  j  ''  ^ 

y/ p/.ioL  ppic 


-t-  (^-  y. 75^/- 6. 757 Yz / :/  >5  X/oy  "  /Y'  75^}-  TS’^^I.ZSS  p/jl 

_  2.  X  <o  •  £o2S  _  W - 2*__ - - -  /  j 

]/  _  *ir|\f/4i>6;  I  ,i/f7  3,206  /  :  y/206  I 

[(/-Z-276A/0  /)J  J  '1-7.cv/^/777X/6  /  IJlI  .  ^ilP  ( n.rs^ApY  j 

/  \  %  /  '  7. 1  t*  i.''^^  ‘'  '.  ‘ 

/,  ^7,0  ^~^V4  ,/„-5Z/ 5/20 ;  ,  /v  _  _  7,1:5 


_  7,7  0/. 

P  h 


2.  X  6/  _  A  /  5  X6;6 


.  c  '7  !'■<'■■■'  • 

ildUJU  "UP, 


\  DANIEL.  MANN.  JOHNSON,  s.  MENDENHALL 


DESIGNED  BY  [fj  VJ  />,  \  DAT 


APPROVED 


INTEGRATINi^  THE  Tli'']E 

m  a2,Al.s.se,x,o'^l, 


2.  ^  ^ t  2  OC 


4(25^^  f7. 2r^. 2^4/2.7,55^/^ 


-5^’/3 


4>  X  7^  2  o  4 


-  2,11 ‘fX id V-  3.  Zi"^ 2 52"  2VO  _ (  y,  2^4/2  .  ? 


'S^j  I  (> .IOC. 


.,  .  f  V  -•S’)'  /  7/i06 

2  5<^-^  ( i<  ^  ii.  7  2, ;  ^,sx>o  1  L 


2  V  H,TOCi  X  7. 204 


V/254  4  2  5^)4  ■1£^22.o.^ 

4  ^  X  i’/TOiOi  h 


-S\lje.70!= 


t  f>7Si/2  -1E2  A 

yioc.  A(  6..  2.0G 


'jh  ^  ISH  (‘i/'Z  zdY 1 .  7  SEX!  0  ^ 


6'  I0L  /7.  24  4 


V-  7/ 2Ty  /V,  2  5^4-  2^-442 . 2 4T//B'^) h  5  y,  7SYfv.2zc.Yz idl^. 

^  -206  /T  i  /  <P,  2(36  'sTX^^X^  9,  2f<i 


^/d:  2  sC(9, 7syY^'  72h)(7 . 2x5  a  4  V 


(?.  3^6 


7  K  L  20^  /723f 


2.^  6.  204  A”  (T.  204 


j9.2S(^12sCU7Ei)(E. 7Eh):7.  ^SSXao^) H 

■2- AT  4,  20G  y2.X  '9/  20l 

^  ("y.  2  sc^y.  7  E^2.  Ts  d)d,  jp/Y E  ysyh,  tsta d  '^' 
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